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ABSTRACT  We  report  here  the  application  of  a 
previously  described  method  to  directly  determine  the 
C02  permeability  (PC02)  °f  the  cell  membranes  of 
normal  human  red  blood  cells  (RBCs)  vs.  those  defi¬ 
cient  in  aquaporin  1  (AQPl),  as  well  as  AQP1 -express¬ 
ing  Xenopus  laevis  oocytes.  This  method  measures  the 
exchange  of  180  between  C02,  HC03“,  and  H20  in 
cell  suspensions.  In  addition,  we  measure  the  alkaline 
surface  pH  (pHs)  transients  caused  by  the  dominant 
effect  of  entry'  of  C02  vs.  HC03-  into  oocytes  exposed 
to  step  increases  in  [C02].  We  report  that  1)  AQP1 
constitutes  the  major  pathway  for  molecular  C02  in 
human  RBCs;  lack  of  AQP1  reduces  Pc(>2  from  the 
normal  value  of  0.15  ±  0.08  (so;  n= 85)  cm/s  by  60%  to 
0.06  cm/s.  Expression  of  AQP1  in  oocytes  increases 
PC()2  2-fold  and  doubles  the  alkaline  pHs  gradient.  2) 
pCMBS,  an  inhibitor  of  the  AQP1  water  channel, 
reduces  PCo2  °f  RBCs  solely  by  action  on  AQPl  as  it 
has  no  effect  in  AQP  1-deficient  RBCs.  3)  PC02  deter¬ 
minations  of  RBCs  and  pHs  measurements  of  oocytes 
indicate  that  DIDS  inhibits  the  C02  pathway  of  AQP1 
by  half.  4)  RBCs  have  at  least  one  other  DIDS-sensitive 
pathway  for  C02.  We  conclude  that  AQPl  is  responsi¬ 
ble  for  60%  of  the  high  PC02  °f  red  cells  and  that 
another,  so  far  unidentified,  C02  pathway  is  present  in 
this  membrane  that  may  account  for  at  least  30%  of 
total  P co2* — Endeward,  V.,  Musa-Aziz,  R.,  Cooper, 
G.  J.,  Chen,  L.,  Pelletier,  M.  F.,  Virkki,  L.  V.,  Supuran, 
C.  T.,  King,  L.  S.,  Boron,  W.  F.,  Gros,  G.  Evidence  that 
aquaporin  1  is  a  major  pathway  for  C02  transport 
across  the  human  erythrocyte  membrane.  FASEB  J.  20, 
1974—1 98 1  (2006) 

Key  Words:  human  red  cell  membrane  •  C02  permeability 
•  pCMBS  •  DIDS 

One  of  the  most  widely  held  assumptions  in  biology 
has  been  that  virtually  all  gases  rapidly  cross  all  cell 
membranes,  and  do  so  simply  by  dissolving  in  mem¬ 
brane  lipids.  The  first  major  challenge  to  this  dogma 
was  the  observation  that  the  apical  membranes  of 
gastric  gland  cells  have  no  demonstrable  permeabil¬ 
ity  to  NHS  or  C02  (1).  The  second  was  the  observa¬ 


tion  that  the  water  channel  aquaporin-1  (AQPl) 
appears  to  serve  as  a  conduit  for  C02  (2-  1),  at  least 
in  artificial  systems.  The  third  was  the  finding  that 
DIDS  markedly  reduces  COL,  permeability  of  the  red 
blood  cell  (RBC)  membrane  (5,  6).  The  first  evi¬ 
dence  for  the  physiological  relevance  of  gas  channels 
was  the  demonstration  that  AQPl  plays  a  key  role  in 
C02  irptake  during  photosynthesis  by  tobacco  plants  (7). 

To  determine  the  C02  permeability  of  red  cells  and 
oocytes,  we  used  a  previously  developed  method  of 
quantitating  the  C02  permeability  of  cells  in  suspen¬ 
sion  by  the  mass  specirometric  lsO  exchange  technique 
(8-10),  which  is  especially  suitable  to  detect  rather 
high  CO«2  permeabilities  that  arc  difficult  to  measure  by 
stopped-flow  techniques.  We  also  developed  a  new 
technique  to  visualize  increases  in  C02  permeation 
across  membranes  by  measuring  alkaline  surface  pH 
transients  associated  with  augmented  transfei  of  molec¬ 
ular  C02  from  the  extra-  to  the  intracellular  space  of 
Xenopus  laevis  oocytes.  From  experiments  in  which  we 
compare  red  blood  cells  from  normal  and  Colton  null 
(AQPl  )  humans,  we  report  that  AQPl  appears  to  he 
responsible  for  60%  of  the  C02  permeability,  a  finding 
confirmed  by  studies  of  COL,  permeability  and  alkaline 
surface  pH  transients  in  AQPl -expressing  oocytes. 


MATERIALS  AND  METHODS 

Red  cells 

Human  red  cells  were  taken  from  several  members  of  our 
laboratory  and  used  for  mass  sped romc trie  experiments  on 
the  same  and  the  following  day.  AQPl -deficient  (Rh -positive) 
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red  tells  were  taken  from  two  individuals  with  the  Colton  null 
phenotype  (of  only  seven  kindreds  identified  in  the  world) 
(11,  12).  All  red  cells  were  washed  three  times  in  physiological 
saline  before  being  used  in  the  mass  spectromctric  experi¬ 
ments  and  were  controlled  for  hemolysis.  Hematocrit  and 
blood  cell  count  were  taken  to  determine  mean  corpuscular 
volumes.  Informed  consent  was  sought  and  given  in  accor¬ 
dance  with  the  Declaration  of  Helsinki,  blood  samples  were 
shipped  chilled  and  red  cells  were  used  within  2-4  days  after 
the  samples  were  taken  It  was  ascertained  by  control  exper¬ 
iment  that  this  prolonged  time  interval  between  blood  with¬ 
drawal  and  determination  of  permeability  does  not  affect 
CO.j  permeability  values  (P(  ()2):  red  cells  from  normal  blood 
investigated  on  the  day  of  withdrawal  showed  P(  (>2  of  0.16  ± 
0.06  cm/s  (w  =  35),  while  another  group  of  normal  red  cell 
samples,  after  having  traveled  for  several  days,  yielded  Pi:o2  of 
0.16  ±  0.07  cm/s  (/t=20). 

Inhibitors 

Most  mass  spectromctric  experiments  with  intact  red  cells 
were  carried  out  in  the  presence  of  the  extracellular  carbonic 
anhydrase  (CA)  sulfonamide  inhibitor  l-[5-sulfamoyl-l,3,4- 
thiadiazol-2-yl-(aminosulfonyl-4-phenyl)  ]-2, 4,6-t  rimethyl-pyri- 
diniinn  perchlorate,  STAPTPP,  described  by  Casey  et  ol.  (13), 
at  a  final  concentration  of  1  *  10  1  M.  The  inhibitor  was 
added  in  order  to  ensure  inhibition  of  small  amounts  of  CA 
that  may  be  set  free  b\  hemolysis  during  the  mass  spectro- 
metric  experiment  due  to  the  continuous  stirring  of  the  red 
cell  suspension.  We  demonstrated  that  the  intracellular  CA 
activity  of  red  cells  was  not  affected  by  30  min  preincubation 
of  red  cells  with  the  same  concentration  of  this  inhibitor. 
Para-(chloi omercurij-bcnzenesulfonate,  or  pCMBS,  is  a 
known  inhibitor  ol  the  water  as  well  as  CO«2  permeability  of 
AQ1M  (3,  12  14).  4,4'-Diisothiocyanato-2,2'-siilbencdisul- 
fonate  (DII)S)  was  used  as  an  established  inhibitor  of  the 
red  cell  Cl  HCO:i  exchanger  AE1  (band  3)  and  other 
ion  transporters,  which  we  have  previously  shown  to  also 
inhibit  red  cell  C02  permeability  (5). 

AQP1  expression  in  Xenopus  laevis  oocytes 

Oocyte  isolation  was  performed  according  to  established 
procedures  (15).  Brielly,  female  Xenopus  laevis  were  anesthe¬ 
tized  using  0.2%  MS-222  (ethyl  3-aminobenzoatc  methanesul- 
fonate;  Sigma,  St.  Louis,  MO,  USA).  Ovarian  lobes  were 
removed  and  placed  in  0  Ca  solution  (08  rnM  NaCI,  2  mM 
KCl,  1  mM  MgCl2,  5  mM  HEPES,  pH  7.5).  Enzymatic  disso¬ 
ciation  and  dcfolliculation  of  oocytes  were  perf  ormed  using  2 
mg/ml  type  1A  collagenase  (Sigma)  in  0  Ca  solution  Stage 
Y-Yrl  oocytes  were  selected  and  kept  at  18°C  in  sterile-filtered 
OR3  medium  containing  (per  2  L)  one  pack  of  powdered 
Leibovitz  L-15  media  with  t -glutamine  (GIBCO-Life  Technol¬ 
ogies,  Inc.,  Gaithersburg,  MD,  USA),  100  ml  of  10,000  U 
penicillin,  10,000  U  streptomycin  solution  (Sigma),  and  5 
mM  HEPES,  pH  7.5. 

Capped  cRNA  was  transcribed  from  a  plasmid  containing 
human  AQP1  (kind  gift  from  Peter  Agie,  Johns  Hopkins 
University)  or  mouse  CA  11  (kind  gift  from  William  S.  Sly, 
Saint  Louis  University)  using  Message  Machine  kit  (Ambion, 
Austin,  TX,  USA).  Oocytes  were  injected  with  50  nl  0.02 
pg/pl  cRNA  encoding  CAH  or  with  a  mixture  containing 
cRNAs  encoding  CAII  (0.02  pg/pl)  and  AQP  I  (0.005  pg/pl). 
Experiments  were  performed  3—6  days  after  injection.  Ex¬ 
pression  of  AQP1  in  oocytes  was  substantiated  by  determining 
the  time  from  their  immersion  in  deionized  water  to  lysis  (2). 
This  time  was  >20  min  in  control  oocytes  and  <1  min  in 
AQP  l -expressing  oocytes.  Similarly,  the  expression  of  CA  was 
verified  by  determining  the  oocytes’  CA  activity  by  mass 


spectrometry,  yielding  an  activity  of  —  10  U  for  control  oocytes 
and  of  150-200  U  for  CAU-exprcssing  oocytes. 

Principle  of  mass  speetrometrie  measurement 

We  adapted  the  technique  of  ltada  and  Forster  (16)  to  study 
the  exchange  of  IH0-labelcd  C02  and  HCO-U  between  an 
extracellular  solution  and  the  intracellular  space  of  cells 
containing  CA  activity,  in  the  present  study  either  red  cells  in 
suspension  or  oocytes  expressing  CA  also  suspended  in  fluid, 
ltada  and  Foster  (16)  showed  that  by  observing  the  time 
course  of  the  species  C,HOlhO  (mass/cliarge,  m/z,  46  peak 
height)  in  the  suspension  with  a  mass  spectrometer,  it  is 
possible  to  determine  the  intracellular  CA  activity  A,  and  the 
bicarbonate  permeability  Pnuo.v-  of  the  cell  membranes.  The 
new  feature  of  the  present  approach  is  the  theoretical  treat¬ 
ment  of  the  process  of  ,hO  exchange,  which  docs  not  assume, 
as  ltada  and  Forster  had,  that  C02  is  infinitely  permeable  and 
cannot  be  a  limiting  step  in  the  overall  exchange  process  (8,  10). 

In  normal  human  red  cells,  C02  permeation  across  the  cell 
membrane  has  a  detectable  effect  on  the  kinetics  of  ihO 
exchange  and  thus  it  is  possible  to  determine  the  membrane 
C02  permeability  from  measurements  of  ,mO  exchange  kinet¬ 
ics  (8,  5)  The  set  of  equations  describing  the  entire  process 
of  membrane  transport  of  ,H0-labeled  C02,  HCO:i  ,  and 
1 120  as  well  as  the  chemical  reactions  between  them  has  been 
reported  (8,  10),  and  is  used  here  to  estimate  Pt;o2  from  the 
kinetics  of  the  extracellular  concentration  of  ClH0H,0  in 
suspensions  of  cells.  These  equations  were  solved  numerically 
and  a  fitting  procedure  was  used  to  find  best-lit  values  of 
membrane  PCQ2  and  P|,(:orU .  The  value  of  the  intracellular 
CA  activity  A,,  which  is  required  for  the  calculations,  was 
determined  by  measuring  by  mass  spectrometry  the  CA 
activity  of  red  cell  lysates  under  conditions  of  pH  and  [Cl  ]  as 
they  prevail  inside  cells. 

Surfaee  pH  and  Pf  measurements  in  oocy  tes 

Others  (17,  18)  had  monitored  transient  changes  in  extracel¬ 
lular  pH  caused  by  the  addition  of  C02/I1C0  .*  to  the 
extracellular  fluid.  We  monitored  the  surface  pi  1  (plls)  of 
Xenopus  oocytes — injected  either  with  50  nl  of  1 12()  or  50  nl  of 
a  solution  containing  RNA  (0.5  ng/nl)  encoding  human 
AQP!  —  using  a  glass  microelectrode  (tip  diameter  of  —  15 
pm),  containing  a  pli-sensitive  ionophore  (Hydrogen  lono- 
phore  l-Cocktail  B  #95293,  Fluka).  When  plls  was  relatively 
stable,  wc  periodically  calibrated  the  pH  electrode  bv  using  a 
model  M PC-200  motorized  inicroinanipulator  (Sutter  Instru¬ 
ment  Co.,  Novato,  CA,  USA)  to  move  the  tip  —300  pm  away 
from  the  oocyte  surface.  At  other  times,  the  electrode  tip 
pushed  up  against  the  oocyte  surface  —60  pm.  The  pi  1  signal 
was  amplified  by  a  model  FD223  high-impedance  electrome¬ 
ter  (World  Precision  Instruments,  Inc.,  Sarasota,  FL,  USA). 
The  external  reference  electrode  was  a  calomel  half  cell  with 
a  3-M  KCl  bridge.  The  C02/HCO>free  extracellular  ND96 
solution  contained  (in  mM):  96  NaCI,  2  KCl,  1  MgCl.,,  1.8 
CaCl2,  5  HEPES,  was  titrated  to  pH  7.50,  and  had  an 
osmolality  of  —200  mosmol/kg  H20.  The  CO., /I  ICO s  solu¬ 
tion  contained  (in  mM):  66  NaCI,  2  KCl,  1  MgCl2,  1.8  ("a Cl.,, 
5  HEPES,  33  NaHCOv 

Wc  measured  water  permeabilities,  Pf,  according  to  the 
method  of  Preston  et  ol.  (19).  /\s  described  by  Virkki  et  ol. 
(20),  wc  used  a  video  camera  to  monitor  the  projection  area 
of  an  oocyte  while  exposing  the  cell  to  a  solution  similar  to 
C02/ H CO-.rfiee  solution,  except  that  we  reduced  [NaCI]  to 
lower  the  osmolality  to  80  niOsm/kg.  A  small  metal  sphere 
next  to  the  oocyte  scivcd  as  a  size  reference.  All  experiments 
were  performed  at  room  temperature  (  — 22°C). 


AQUAPORIN  1  IS  A  PATHWAY  FOR  CO 


1975 


RESULTS 

C02  permeability  of  the  human  red  cell  membrane 

Using  a  large  number  of  measurements  with  red  cells  of 
healthy  donors,  we  determined  a  P(02  of  0.15  ±  0.08 
(so)  cm/s  for  the  human  red  cell  membrane  (?t=85). 
The  median  in  this  group  of  data  is  0.14  cm/s.  We  note 
that  the  maximal  PCQ2  that  can  be  determined  with  the 
present  technique  for  the  conditions  of  the  human  red 
cell  is  —  I  cm/s.  The  same  group  of  85  measurements 
yielded  a  Pnc.o:*-  of  1.3  •  10  3  ±  0.3  *  10  3  (so)  cm/s, 
with  a  median  of  1.3  *  10  3  cm/s. 

Effect  of  pCMBS  on  C02  permeability 

Figure  1  illustrates  an  experimental  series  studying  the 
effects  of  various  concentrations  of  pCMBS  on  PHcos- 
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Figure  1.  Semilogarithmic  plot  of  a  titration  of  P| it. os-  a,)d 
P(;02  in  human  red  cells  by  pCMBS.  Intracellular  CA  activity 
in  all  experiments  was  taken  to  be  identical  to  that  of  red  cells 
in  the  absence  of  pCMBS.  Number  of  experiments  10  for 
control  (P(  02  =  0.12  cm/s,  not  shown  in  the  graph),  3  for  0. 01 
mM  pCMBS,  1 1  for  0.05  mM,  9  for  0.1  mM,  3  for  0.2  mM,  14 
for  1  mM,  5  for  2  mM,  and  5  for  3  mM.  Bars  =  so.  The  IW)  of 
the  pCMBS  effect  on  P<X)2  can  be  estimated  to  be  —0.5  mM. 
The  Ir>0  of  the  pCMBS  effect  on  Pn<  orv_  is  known  from  the 
literature  to  be  2  mM.  ^Significantly  different  from  control 
value  (/*<0.05). 


and  P ( ;o2*  pCMBS  is  an  efficient  inhibitor  of  the  water 
permeability  of  AQPI  (14),  but  is  also  known  to  inhibit 
band  3  protein  to  some  extent  (1D50  for  pCMBS 
binding  to  band  3  was  reported  to  be  2.0  mM  by  Zhang 
and  Solomon,  ref.  (21).  If  PC02  of  Colton  null  red  cells, 
0.06  cm/s  (Fig.  2),  is  taken  as  the  value  to  be  expected 
with  full  inhibition  of  the  AQPI  C02  pathway,  then  the 
data  of  Fig.  1  lead  to  an  estimate  of  lr>0  for  inhibition  of 
this  pathway  of  —0.5  mM  pCMBS.  This  is  somewhat 
higher  than  t  he  Kj  estimated  for  inhibition  of  the  AQPI 
water  pathway  by  pCMBS  of  0.13  mM  (14),  but  is 
markedly  lower  than  the  150  of  pCMBS  for  band  3.  In 
conclusion,  pCMBS  inhibits  the  C02  permeability  of 
human  red  cells  at  concentrations  of  —4-lbld  lower 
than  needed  for  inhibition  of  band  3  and  it  will  be 
shown  below  that  this  occurs  by  the  action  of  pCMBS 
on  AQPI. 

C02  permeability  of  aquaporin  1-deficient  red  cells 

When  Colton  null  (AQP1  )  red  cells  are  compared 
to  normal  human  red  cells,  Phcos-  is  identical  in  both 
types  of  red  cells  (Fig.  2 A)  hut  P(Q2  is  significantly  (by 
—60%)  lower  in  Colton  null  cells  than  in  control  red 
cells  (Fig.  2/i),  indicating  that  AQPI  contributes  signif¬ 
icantly  to  CO,2  permeation  across  the  red  cell  mem¬ 
brane.  Consistent  with  this,  the  AQPI  inhibitor  pCMBS 
at  concentrations  of  1  and  2  mM  has  no  significant 
effects  on  PC(J2  of  AQPI -deficient  red  cells,  but  it 
significantly  reduces  Pco2  of  control  cells  (Fig.  2 Zi). 
Even  at  the  higher  pCMBS  concentration  of  3  mM, 
where  the  inhibitory  effect  of  pCMBS  on  band  3 
becomes  apparent  by  a  significant  reduction  of  PIlco.s_ 
to  —1/2  both  in  control  and  Colton  null  red  cells  (Fig. 
2 A),  this  drug  still  produces  no  effect  on  P(:02  in  Colton 
null  red  cells  (Fig.  2/i).  We  conclude  that  pCMBS  at  all 
concentrations  studied  is  without  effect  on  P(;c)2  of 
AQPI -deficient  red  cells,  confirming  that  pCMBS  re¬ 
duces  PC02  solely  by  its  action  on  AQPI. 

We  have  shown  that  the  AF1  inhibitor  D1DS  mark¬ 
edly  reduces  PC02  as  well  as  PHCOH_  in  normal  red  cells 
(5).  D1DS  reduces  Pj-jco.s-  lo  a  similar  extent  in  Colton 
null  and  in  control  cells  (Fig.  2A),  as  expected.  On  the 
other  hand,  PC02»  which  is  reduced  in  Colton  null  cells 
even  in  the  absence  of  D1DS,  is  further  reduced  by 
1)1  DS  to  —0.012  cm/s,  —  1  /5th  of  the  figure  obtained  in 
the  absence  of  DIDS,  The  lowest  Pc:Q2  value  seen  in 
Fig.  2/i,  the  value  for  Colton  null  cells  treated  with 
DIDS  (0.012  cm/s),  amounts  to  only  1  / 1 2t h  of  the 
P(:o2  found  in  normal  red  cells.  In  conclusion,  P(:o2  of 
AQPI -deficient  red  cells  is  reduced  compared  to  con¬ 
trol  cells;  it  is  further  reduced  by  DIDS  but  is  not 
affected  by  pCMBS.  The  inhibitory  effect  of  DIDS 
lowers  1\ 02  lu  normal  red  cells  by  0.10-0.1 1  cm/s,  but 
lowers  Pct)2  of  Colton  null  red  cells  by  only  0.04-0.03 
cm/s.  These  data  suggest  that  part  of  the  inhibitory 
effect  of  DIDS  occurs  by  an  action  on  AQPI,  whereas 
the  other  part  results  from  its  action  on  a  different 
DIDS-sensitive  pathway. 
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Figure  2.  Comparison  of  the  effects  of  D1DS  and  pCMBS  on 
P, ,(  ( rv_  and  PC()2  of  normal  human  red  ceils  and  on  Colton 
null  human  red  cells.  Number  of  blood  samples  investigated 
was  >  5  for  each  condition;  number  of  donors  of  Colton  null 
blood  was  2.  It  is  seen  that  P(:02  is  significantly  reduced  when 
aquaporin-I  is  absent.  DIDS  reduces  P(:C)2  of  both  normal  and 
Colton  null  red  cells.  pCMBS  reduces  Pc:c)2  of  normal  but  not 
of  Colton  null  red  cells.  *Significant  differences  to  the 
leftmost  black  column  (/*<0.05),  "Significant  differences  to 
the  leftmost  gray  column  (/*< 0.05).  Error  bars  =  sn. 

Surface  pH  transients  of  AQP1  -expressing  oocytes 
and  effect  of  DIDS 

To  test  the  hypothesis  that  DIDS  reduces  the  C02 
permeability  of  AQPl ,  we  examined  the  effect  of  DIDS 
in  oocytes  either  expressing  human  AQPl  or  injected 
with  H20.  Preliminary  experiments  (not  shown)  indi¬ 
cated  that  100  jxM  DIDS  reduces  the  rate  at  which  the 
addition  of  C02  causes  intracellular  pH  to  fall  in 
oocytes  with  vitelline  membrane  removed.  Figure  3A,  B 
illustrates  more  sensitive  experiments  on  oocytes  (with 


vitelline  membrane  intact),  in  which  we  monitored  the 
transient  rise  in  extracellular  surface  pH  caused  by 
switching  to  an  extracellular  solution  containing  5% 
C02/33  niM  HCOf  (pH  7.50).  Adding  C02/HC03" 
causes  a  rapid  rise  in  pHs,  followed  by  an  approxi¬ 
mately  exponential  decay.  This  pHs  trajectory  reflects 
the  effect  of  C02  influx  (less  any  effect  of  HCO.,~ 
influx)  on  the  equilibrium  C02  +  H20  HCO.f  + 
H  near  the  cell  surface.  The  spike  height,  an  index  of 
maximal  C02  influx,  was  substantially  less  in  the  H2Q- 
injected  oocyte  (Fig.  3 A)  than  in  the  AQPl -expressing 
oocyte  (Fig.  3/J),  confirming  that  AQPl  increases  C02 
permeability.  After  removal  of  C02/HCO;f,  we  applied 
100  |x\l  DIDS  and — in  the  continuous  presence  of 
DIDS — exposed  each  oocyte  to  C02/HCO;t~  a  second 
time.  DIDS  had  little  effect  on  spike  height  in  the 
H20-injected  oocyte  but  reduced  spike  height  in  the 
AQPl -expressing  cell 

Figure  3C  summarizes  similar  paired  data  from  a 
total  of  six  experiments  each  on  H20-injected  and 
AQPl -expressing  oocytes,  and  indicates  that  expression 
of  AQPl  increases  spike  height  by  96%  and  that  DIDS 
reduces  the  AQP  1-dependent  spike  height  by  49%. 
Figure  3 1)  summarizes  experiments  in  which  we  found 
that  DIDS  has  no  effect  on  Pf.  Thus,  DIDS  substantially 
reduces  the  C02  permeability  of  AQPl  without  affect¬ 
ing  11,0  permeability. 

Effects  of  aquaporin  1  expression  on  C02 
permeability  of  Xenopus  laevis  oocytes 

Figure  4  gives  the  results  of  Pc:o2  determinations  ob¬ 
tained  with  oocytes  in  suspension.  In  oocytes  expressing 
CA  11  only,  we  measured  a  PC02  of  0.057  cm/s,  while 
coexpressing  both  CA  and  AQPl  raised  P(:()2  signifi¬ 
cantly  to  0.1  1  cm/s.  This  agrees  with  the  finding  of  Fig. 
3  of  a  marked  increase  in  the  alkaline  surface  pH 
transient  due  to  the  expression  of  AQPl  in  oocytes.  The 
increase  of  Pc02  in  AQPl -expressing  oocytes  comple¬ 
ments  the  finding  of  a  reduced  PC02  in  aquaporin-1- 
deficient  red  cells. 


DISCUSSION 

The  CO 2  permeability  of  normal  human  red  cells 

The  value  of  0.15  cm/s  for  P(;o2  reported  here  may 
somewhat  underestimate  the  true  C02  permeability  of 
the  red  cell  membrane  for  two  possible  reasons.  First, 
the  theoretical  treatment  of  the  tH0  exchange  process 
assumes  perfect  mixing  in  the  intracellular  space,  and 
second,  it  assumes  that  unstirred  layers  around  the  red 
cells  have  no  effect  on  the  calculated  PC02.  The  latter 
assumption  is  supported  by  our  previous  conclusion 
(10)  that  the  present  technique  is  highly  insensitive 
toward  unstirred  layers,  likely  because  there  is  a  store  of 
bicarbonate  in  the  external  solution,  including  the 
membrane  surface,  which  can  replenish  C02  that  has 
permeated  the  membrane  into  the  cell  interior.  As  the 
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Figure  3.  Surface  pH  and  Pf  measurements  in  Xenopus  Uievis  oocytes  exposed  to  a  solution  containing  5%  C02/33  mM  HCO“;i 
at  pi  I  7.50.  A)  Time  course  of  pi  Is  in  a  H20-injectcd  oocyte.  During  the  second  pulse,  100  pM  DIDS  was  present.  The  step-graph 
below  the  plls  record  indicates  when  the  electrode  was  on  the  oocyte  surface,  and  when  it  was  in  the  bulk  extracellular  (pi  I 
7.50)  solution.  B)  Time  course  of  plls  in  an  AQ.P1  -expressing  oocyte.  C)  Summary  of  paired  plls  data.  P  values  indicate  the 
results  of  paired,  two-tailed  t  tests.  D)  Summary  of  Pf  data.  The  effect  of  D1I)S  was  not  significant  in  either  H^O-injected  or 
AQP 1 -expressing  oocytes.  Error  bars  =  si). 


Figure  4.  Effect  of  the  expression  of  carbonic  anhvdrase  II 
alone  (CA  II)  and  of  CA  II  and  aquaporin-l  together  (CA  II 
+  AQP  I)  on  the  C02  permeability  of  Xenojms  laevis  oocytes. 
n  =  33,  error  bar  =  so,  *P  <  0.04. 


former  assumption  is  likely  not  to  be  entirely  fulfilled, 
the  apparent  PCX)2  will  tend  to  be  an  underestimate.  It 
may  be  noted  that  from  steady-state  measurements  of 
CO*2  fluxes  across  layers  of  red  cells,  we  earlier  postu¬ 
lated  a  Pc02  of  1-2  cm/s  (22,  23). 

The  lipid  phase  of  the  membrane  is  likeh  to  impart 
a  “basal”  C02  permeability  to  the  membrane,  which 
may  be  represented  by  the  value  of  0.012  cm/s  ob¬ 
served  in  AQP  1 -deficient  red  cells  exposed  to  DIDS 
(Fig.  2/?).  This  latter  value  represents  the  lowest  C02 
permeability  observed  in  this  study  and  would  be 
30-fold  lower  than  the  C02  permeability  of  0.3  cm/s 
observed  by  Gutknccht  et  ai  (24)  in  artificial  lipid 
bilayers,  whose  composition,  however,  was  not  repre¬ 
sentative  of  the  lipids  in  the  red  cell  membrane  and 
which  were  devoid  of  membrane  proteins.  Our  data 
indicate  that  AQP1  and  an  additional  D1  OS-sensitive 
membrane  protein  markedly  enhance  red  cell  P(  Q2  to 
the  point  that  the  resistance  of  the  membrane  for  CO., 
diffusion  is  at  least  two  orders  of  magnitude  lower  than 
that  for  HCO.,~. 

Bicaibonatc  permeabilities  reported  here  are  1.5-  to 
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3-fold  higher  than  values  reported  in  the  literature,  and 
tlieii  inhibition  by  D1DS  appears  less  than  has  been 
established  by  other  methods.  This  is  due  to  some  lysis 
that  occurs  by  the  stirring  during  the  mass  spectromet- 
ric  experiment.  It  can  be  shown  theoretically  that  the 
CA  set  free  by  this  lysis  will,  before  the  inhibition  by 
STAPTPP  is  complete,  effect  an  apparent  increase  in 
PH(  of  up  to  4-fold,  but  will  cause  an  overestimation 
of  P(Q2  by  no  more  than  a  few  percent. 

Aquaporin-l-mediated  permeation  of  C02 
across  membranes 

Affuaponn-l  constitutes  a  pathway  for  C02 

The  absence  of  AQPl  from  the  red  cell  membrane 
reduces  C02  permeability  from  0.15  cm/s  to  0.06  cm/s 
(i.o.,  by  60%).  The  abundance  of  band  3,  Glutl,  the 
urea  transporter,  and  spectrin  are  similar  in  AQP1- 
deficient  and  normal  human  red  cells  (12),  supporting 
our  conclusion  that  the  decrease  observed  in  PCQ2  must 
be  attributed  to  aquaporin-1.  This  conclusion  is  con¬ 
firmed  here  by  two  other  findings  that  were  obtained  by 
entirely  different  experimental  approaches:  1)  the 
2-fold  increase  in  P(:02  of  the  oocyte  membrane  on 
AQPl  expression  as  determined  by  mass  spectrometry, 
and  2)  the  marked  increase  of  the  alkaline  surface  pH 
transient  on  AQP1 -expressing  oocytes.  Our  findings  for 
P(:()2  in  oocytes  are  consistent  with  earlier  studies 
examining  C02  influx  into  oocytes  from  the  time 
course  of  intracellular  pH  (2,  3). 

There  has  been  some  controversy  as  to  whether 
aquaporin-1  can  act  as  a  C02  pore.  Yang  et  al.  (25) 
observed  the  intraeiythrocytic  acidification  on  mixing 
of  red  cells  with  a  CO<2  solution  in  a  rapid  reaction, 
stopped-flow  spectrophotometer  and  obtained  similar 
Pf:o2  values  for  red  cells  from  wild-type  (YVT)  mice 
(0.012  cm/s)  and  for  red  cells  from  AQP1  -deficient 
mice  (0.01 1  cm/s).  Similarly,  they  found  identical  P(:02 
values  of  10  *  cm/s  in  liposomes  with  and  without 
reconstitution  of  AQPl  using  stopped-flow.  They  con¬ 
cluded  that  aquaporin  1  does  not  act  as  a  C02  channel. 
Of  note,  P(:02  of  10  s  cm/s  for  liposomes  is  100-fold 
lower  that  the  value  reported  by  Gutknecht  et  al.  (24), 
and  the  P(;{)2  value  observed  for  red  cells  is  >  10-fold 
lower  than  the  figures  r  eported  here.  These  differences 
could  be  due  to  considerable  unstirred  layers  around 
the  cells  or,  more  likely,  to  inadequate  mixing  in  the 
stopped-flow  chamber  that  caused  an  apparently  low 
C02  permeability;  therefore,  the  contribution  of  AQPl 
was  not  evident.  Our  findings  arc  in  line  with  those  of 
Prasad  et  al.  (4),  who  used  stopped-flow  spectrophotom¬ 
etry’  to  demonstrate  a  C02  permeability  in  liposomes  of 
0.5  cm/s  and  in  AQPl  reconstituted  protcoliposomcs 
of  1.9  cm/s,  a  nearly  4-fold  increase  that  was  inhibited 
by  the  aquaporin  inhibitor  HgCl2. 

pCMBS  inhibits  the  C02  permeability  via  AQPl  only 

The  present  data  confirm  previous  reports  that  pCMBS 
inhibits  a  pathway  for  molecular  C02  across  aquaporin 


1 -expressing  oocyte  membranes  (3)  and  across  normal 
human  red  cell  membranes  (6).  The  data  of  the 
present  paper  show  that  in  red  cells  this  pathway  is 
constituted  by  AQPl:  pCMBS,  an  established  inhibitor 
of  the  water  permeability  of  AQPl  (14),  at  concentra¬ 
tions  between  2  and  3  mM  reduces  Pro2  of  normal  red 
cells  by  50-70%  (Figs.  1 , 2)  but  has  no  effect  on  Pro2  of 
Colton  null  red  cells. 

DIDS  inhibits  partially  the  C02  pathway  via  AQPl 

In  normal  human  red  cells  DIDS  reduces  Pro2  from 
0.15  cm/s  to  0.05-0.04  cm/s,  producing  a  AP(  ()2  of 
~0. 1  cm/s,  whereas  in  Colton  null  red  blood  cells  it 
reduces  Pc;02  from  0.06  cm/s  to  0.012  cm/s,  causing  a 
APro2  of  ~0.05  cm/s  (Fig.  2 B).  This  is  most  simply 
interpreted  to  indicate  that  one-half  of  the  DIDS  effect 
is  due  to  inhibition  of  AQPl  by  DIDS,  causing  a  AP(()2 
of  0.05  cm/s,  and  that  the  other  half  of  the  DIDS  effect 
is  due  to  an  effect  other  than  that  on  AQPl.  This  is 
most  clearly  substantiated  by  the  inhibitory’  effect  of 
DIDS  on  the  surface  pH  transient  observed  in  aqua- 
porin-expressing  oocytes  (Fig.  3). 

DIDS  inhibits  in  addition  a  C02  pathway  that  is  not 
linked  to  AQPl 

The  effect  of  DIDS  on  P(:Q2  seen  in  AQPl -deficient  red 
cells  (Fig.  2 B)  shows  clearly  there  is  another  pathway 
for  C02  in  the  red  cell  membrane  that  l)  is  inhibitable 
by  DIDS,  2)  is  independent  of  AQPl ,  and  3)  probably  is 
not  inhibited  by  pCMBS.  Although  it  cannot  be  ex¬ 
cluded  that  this  pathway  is  constituted  by  the  lipid 
bilayer  and  DIDS  acts  in  an  unknown  manner  on  its 
lipid  phase,  as  speculated  by  Forster  et  al.  (5),  it  should 
be  noted  that  DIDS  has  no  effect  in  oocytes  that  do  not 
express  AQPl  (Fig.  3),  and  thus  it  appears  more  likely 
that  this  CO<2  pathway  is  provided  by  another  1)1  DS- 
inhibitable  membrane  protein  of  the  red  cell,  such  as 
AE1  (band  3),  the  Cl-HCOs  '  exchanger.  The  identity 
of  this  protein  has  yet  to  be  determined.  From  the 
numbers  given  above,  its  contribution  may  make  up 
32%  of  normal  C02  permeability,  while  the  basal 
permeability  due  to  the  membrane  lipids  may  be  re¬ 
sponsible  for  8%. 

Pathway  of  C02  across  the  AQPl  molecule 

What  is  the  pathway  in  the  aquaporin  1  honioletramcr 
through  which  the  C02  molecule  passes?  It  is  well 
established  that  within  the  AQPl  tetramer  each  mono¬ 
mer  possesses  a  water  channel.  Some  have  proposed 
that,  in  addition,  AQPl  exhibits  a  nonsclcctivc  cation 
conductance  mediated  by  the  single  central  pore  of  the 
homotetramer  (26).  Each  AQPl  monomer  possesses  a 
curvilinear  water  channel,  which  exhibits  wide  open¬ 
ings  on  the  extracellular  and  the  cytoplasmic  sides,  and 
within  the  membrane  narrows  down  to  a  size-selective 
pore  in  the  center  of  the  bilayer.  This  channel  is 
~2.8-4  A  at  its  minimal  diameter  and  has  a  length  of 
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18-20  A  (27-29).  Half  of  the  channel  surface  is  hydro¬ 
philic  and  the  other  half  hydrophobic.  Water  with  its 
size  of  2.8  A  is  thought  to  strip  its  waters  of  hydration  via 
several  water  binding  regions  along  the  pore,  thus 
allowing  passage  of  a  single  water  molecule  through  the 
pore.  Ions  are  excluded  on  the  basis  of  their  charge 
and,  in  the  hydrated  form,  their  size  (29).  The  channel 
does  not  conduct  protons.  The  region  around  the 
4-fold  symmetry  axis  of  the  tetramer  has  been  proposed 
to  constitute  the  cation  pore  of  aquaporin  (26).  It 
exhibits  the  narrowest  constriction  close  to  the  external 
side  of  the  pore;  with  an  ~3  A  diameter,  its  interior  is 
largely  hydrophobic.  It  possesses  a  large  central  cavity 
near  the  membrane  midpoint  and  has  many  similarities 
with  K+  channels  (26-28).  The  cation  channel  appears 
to  be  cGMP-gated.  Carbon  dioxide  with  a  MW  of  44  is 
larger  than  water  with  its  MW  of  18,  but  has  an  effective 
molecular  diameter  between  2.8  and  3.4  A,  which  is 
only  slightly  greater  than  the  diameter  of  water  of  2.8  A 
(30,  31). 

because  C02  is  a  linear  molecule  with  a  smaller 
transversal  than  longitudinal  diameter,  it  may  be  ex¬ 
pected  to  be  able  to  permeate  either  of  the  two  types  of 
pores,  each  of  which  has  a  constriction  site  of  ~3  A. 
COL>  is  a  quadrapole,  that  is,  each  of  the  two  C-O  bonds 
has  a  small  charge  separation.  This  may  explain  the 
observation  that  C02,  which  is  highly  soluble  in  water 
and  even  more  soluble  in  nonpolar  solvents  such  as 
benzene,  toluene,  and  heptane,  exhibits  its  highest 
solubility  in  solvents  such  as  acetone  or  ethanol,  which 
contain  both  hydrophobic  and  strongly  polar  groups 
(32).  In  view  of  these  properties,  it  is  reasonable  to 
propose  that  C02  can  transit  the  water  channel  of 
AQP1  with  its  mixed  hydrophobic-hydrophilic  surface. 
This  conclusion  would  be  in  line  with  the  observation 
by  Cooper  and  Boron  (3),  who  reported  that  pCMBS, 
though  significantly  reducing  the  rate  of  C02  uptake 
and  intracellular  acidification  rate  in  AQP1  expressing 
oocytes,  does  not  affect  CO<2  permeation  in  the  AQP1 
mutant  with  serine  replacing  the  mercury-sen  si  tive  cys¬ 
teine  (C189S)  (33).  Since  Cl 89  is  close  to  the  inner 
surface  of  the  water  channel  of  the  AQP1  monomer 
(28)  and  inhibits  the  water  channel  when  mercury 
binds  (33),  this  observation  suggests  that  C02  and  H20 
pass  through  the  same  channel.  On  the  other  hand,  the 
data  in  Fig.  3 — showing  that  D1DS  blocks  half  of  the 
AQPl-dcpcndent,  C02-induced  pHs  transient  in  oo¬ 
cytes  but  has  no  effect  on  water  permeability — suggests 
that  half  of  the  C02  passes  through  AQPl’s  central 
pore  and  half  passes  through  the  four  water  pores.  The 
absence  of  an  effect  of  DIDS  on  water  permeability  in 
the  present  study  is  consistent  with  the  earlier  observa¬ 
tion  that  DIDS  does  not  inhibit  the  water  permeability 
of  red  cells  (34,  35). 

The  apparent  incongruity  of  the  pCMBS  and  DIDS 
data  in  oocytes  could  have  two  explanations  that  are 
not  mutually  exclusive.  First,  because  the  earlier  intra¬ 
cellular  pH  measurements  (3)  are  less  sensitive  than 
the  current  pHs  measurements,  it  is  possible  that  the 
earlier  work  could  not  distinguish  a  50%  inhibition  by 


pCMBS  of  C02  permeability  from  full  blockade.  Sec¬ 
ond,  it  is  possible  that  the  interaction  of  the  large 
pCMBS  molecule  with  Cl 89  in  the  mouth  of  the  water 
pore  causes  a  conformational  change  that  lowers  the 
C02  permeability  of  the  central  pore.  We  add  that, 
from  the  molecular  dimensions  of  02  (effective  diam¬ 
eter  2. 3-2. 9  A;  ref.  30) ,  it  is  conceivable  that  this  gas  can 
also  pass  through  the  same  channel (s),  but  so  far  this 
has  not  been  tested. 

Physiological  role  of  the  AQP1  pathway  for  C02 

Recently,  Uehlcin  et  at.  (7)  have  shown  that  aquaporin 
of  the  tobacco  plant  acts  as  a  C02  pore  and  that  this 
pore  plays  a  physiological  role  in  the  plant's  uptake  of 
C02  and  growth.  The  present  work  is  the  first  demon¬ 
stration  of  the  physiological  role  of  AQP1  as  a  channel 
for  a  dissolved  gas  in  an  intact  animal  cell.  For  normal 
individuals  at  rest,  the  02  and  C02  in  pulmonary 
capillaiy  blood  is  believed  to  equilibrate  with  alveolar 
air  within  the  first  200  ms  of  the  capillary  transit  time  of 
700  ms.  To  estimate  the  role  of  the  gas  diffusion 
resistance  of  the  red  cell  membrane,  we  have  done  a 
simple  model  calculation  in  which  the  intracellular 
diffusion  resistance  is  estimated  by  treating  the  red  cell 
as  a  plane  sheet  of  1.6  jxm  thickness  and  using  the 
known  intracellular  C02  diffusion  coefficient  (22,  36); 
to  this  the  inverse  of  the  membrane  permeability  is 
added  as  a  second  diffusion  resistance  for  C02.  The 
chemical  reaction  inside  the  red  cell  is  assumed  to  be 
infinitely  fast,  which  allows  one  to  describe  the  diffu¬ 
sion  and  reaction  process  by  replacing  the  solubility  a 
in  the  diffusion  equation  by  p,  a  value  defined  as  the 
change  in  total  intracellular  C02  concentration  (C02 
4-  HCO/  )  per  C02  partial  pressure  difference,  i  e.,  the 
slope  of  the  C02  binding  curve  of  red  cells  (37).  The 
time  required  by  Cl  -HC03~  exchange  is  neglected  in 
this  model.  We  have  calculated  in  this  manner  the  time 
necessary  for  C02  release  by  a  red  cell  to  reach  95% 
completion  (t<>5%),  using  the  CQ2  permeabilities  re¬ 
ported  here.  Because  PCQ2  values  may  be  somewhat 
underestimated  due  to  the  assumption  of  perfect  intra¬ 
cellular  mixing  in  the  calculations  of  P(:02,  the  95% 
times  given  in  the  following  may  be  somewhat  overes¬ 
timated,  but  this  does  not  severely  affect  the  compari¬ 
sons  made.  With  a  normal  PC02  of  0.15  cm/s,  t<)r>%  is 
110  ms,  and  is  increased  to  200  ms  when  P(X)2  assumes 
a  value  of  0.06  cm/s  in  the  absence  of  AQP1.  In  view  of 
the  capillar)7  transit  time  of  700  ms  under  resting 
conditions,  it  is  not  expected  that  AQP1  deficiency 
causes  an  impairment  of  C02  exchange  in  the  capillaiy, 
not  even  with  the  reduced  transit  time  that  prevails 
during  exercise.  However,  when  the  AQP1 -indepen¬ 
dent,  DIDS-sensitivc  pathway  for  C02  is  also  inhibited 
and  P co2  falls  to  0.01  cm/s,  one  obtains  U)T>%  —  990  ms. 
Thus,  if  both  pathways  for  molecular  C02  are  absent 
and  PCo2  assumes  its  “basal”  value,  possibly  represent¬ 
ing  C02  diffusion  across  the  lipids  of  the  membrane, 
C02  release  in  resting  individuals  may  be  less  than 
complete  by  the  end  of  the  pulmonary  capillaiy.  Under 
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conditions  of  exercise,  when  capillary  transit  tune  in 
the  lung  is  reduced,  C02  release  may  be  subject  to  a 
marked  diffusion  limitation.  [FJ] 
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Abstract  Others  have  shown  that  exposing  oocytes  to  high 
le\els  of  NH3/NH4 f  ( 10-20  111M)  causes  a  paradoxical  fall 
in  intracellular  pH  (pH,),  whereas  low  levels  (e.g.,  0.5  mM) 
cause  little  pH,  change.  Here  we  monitored  pHj  and  extra¬ 
cellular  surface  pH  (pHs)  while  exposing  oocytes  to  5  or 
0.5  mM  NH3/NH4+.  We  confirm  that  5  mM  NH3/NH4+ 
causes  a  paradoxical  pH,  fall  (— ApH,  =  0.2),  but  also 
observe  an  abrupt  pHs  fall  (— ApHs  =  0.2) — indicative  of 
NH4  influx — followed  by  a  slow  decay  Reducing  [NH3/ 
NH4+]  to  0.5  mM  minimizes  pH,  changes  but  maintains  pHs 
changes  at  a  reduced  magnitude.  Expressing  AmtB  (bacterial 
Rh  homologue)  exaggerates  — ApHs  at  both  NH3/NH4+ 
levels.  During  removal  of0.5  or  5  mM  NH4/NH4+,  failure  of 
pHs  to  markedly  overshoot  bulk  extracellular  pH  implies 
little  NH-*  efflux  and,  thus,  little  free  cytosolic  NH4/NH4‘‘".  A 
new  analysis  of  the  effects  of  NH3  vs.  NH4+  fluxes  on  pH$ 
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and  pH,  indicates  that  (a)  NH^  rather  than  NH4+  fluxes 
dominate  pHj  and  pHs  changes  and  (b)  oocytes  dispose  of 
most  incoming  NHV  NMR  studies  of  oocytes  exposed  to 
,5N-1abeled  NH3/NH4+  show  no  significant  formation  of 
glutamine  but  substantial  NH3/NH4  +  accumulation  in  what 
is  likely  an  acid  intracellular  compartment.  In  conclusion, 
parallel  measurements  of  pH,  and  pHs  demonstrate  that  NH3 
flows  across  the  plasma  membrane  and  provide  new  insights 
into  how  a  protein  molecule  in  the  plasma  membrane— 
AmtB — enhances  the  flux  of  a  gas  across  a  biological 
membrane. 

Keywords  NH3  permeability  •  Surface  pH  measurement 
Xenopus  oocytes  •  AmtB 

The  movement  of  NH3  across  cell  membranes  is  important 
for  several  physiological  processes,  including  nitrogen 
metabolism  by  the  liver  and  acid-base  transport  by  the 
kidney.  In  1897,  Overton — monitoring  the  precipitation  of 
tannins  in  the  algae  Spirogyra — demonstrated  that  it  is  NH^ 
rather  than  NH4+  that  readily  crosses  the  cell  membrane 
Later,  Warburg  (1922),  Harvey  (191 1),  and  Jacobs  (1922) 
reached  similar  conclusions  working  on  a  variety  of  prep¬ 
arations  and  using  different  approaches  to  show  that  the 
influx  of  NH4  produces  a  rise  in  internal  pH  (pH,).  More 
recently,  work  with  pH-sensilive  microelectrodes  by  Boron 
and  De  Weer  (1976b)  on  squid  giant  axons  demonstrated 
that,  although  the  influx  of  the  weak  base  NH4  causes  a 
large  and  rapid  increase  in  pH,,  the  lower  influx  of  the 
weak  acid  NH4+  produces  a  slow  fall  in  pH,  during  the 
“plateau  phase”  of  the  NH3/NH4^  exposure.  Moreover, 
this  influx  of  NH4+  during  the  exposure  to  extracellular 
NH3/NH4+  leads  to  large  undershoot  of  the  original  pH,, 
once  the  NH3/NH44  is  removed  from  the  extracellular 


*?)  Springer 


16 


R.  Musu-Aziz  el  al.:  NH3  Pcnncability  of  oocytes 


solution.  This  effect  is  the  basis  of  the  widely  used 
“ammonium  prepulse”  technique  that  they  introduced 
(Boron  and  De  Weer  1976a,  1976b). 

Aickin  and  Thomas  (1977)  subsequently  showed  that  the 
uptake  of  NH4f  can  be  so  powerful  in  mammalian  skeletal 
muscle  that  the  alkalinizing  effect  of  NH^  entry  is  barely 
detectable.  Kikeri  et  al.  (1989)  made  similar  observations 
when  introducing  NH3/NH44  into  the  lumen  of  the  renal 
thick  ascending  limb.  Waisbren  et  al.  (1994)  were  the  first  to 
delinitively  identify  a  gas-impermeable  membrane,  show¬ 
ing  that  neither  NH3  nor  NH4+  (nor  C02  nor  HCO3”)  could 
cross  the  apical  membranes  of  gastric  gland  cells. 

In  small-diameter  A 'e nopus  oocytes  (Keicher  and  Meech 
1994),  an  exposure  to  20  mM  extracellular  NH3/NH44 
produces  the  classic  biphasic  rise  in  pH],  followed  by  a  fall. 
However,  in  large-diameter  oocytes  (Burckhardt  and 
Fromter  1992;  Keicher  and  Meech  1994),  an  exposure  to 
20  mM  NH3/NH4+causes  a  paradoxical  fall  in  pH,  and  a 
strong  positive  shift  in  membrane  potential  (Vm).  To 
explain  the  above  effect,  Burckhardt  and  Fromter  hypoth¬ 
esized  that  NH4^  enters  via  a  nonselective  cation  channel, 
dissociates  to  form  NH3  I  H+  possibly  followed  by 
sequestration  of  NH3  in  lipid  stores.  This  model  calls  for 
negligible  NH3  permeability  and  accounts  for  the  then- 
available  data.  Keicher  and  Meech  found  that  the  NH3/ 
NH44 -induced  fall  in  pH,  requires  extracellular  Cl“  and 
proposed  that  the  fall  in  pH,  is  due  in  part  to  an  Na/K/Cl 
cotransporter  that  carries  NH4^  in  place  of  K\  However, 
the  acidification  was  only  slightly  inhibited  by  Na+  removal 
or  by  bumetanide.  Moreover,  Na/NH4/Cl  cotransport  would 
not  account  for  the  positive  shift  in  Vm,  which  was  not 
abolished  by  Cl  removal.  Note  that,  in  the  Keicher-Meech 
model,  the  pH,  decrease  would  require  that  the  NH4+  influx 
be  accompanied  by  an  NH3  efflux,  as  outlined  originally  by 
Boron  and  De  Weer,  and  lead  to  a  buildup  of  cytosolic 
NH4+.  However,  if  the  oocytes  could  mediate  NH^  efflux, 
then  NHi/NH44  removal  would  lead  to  a  rapid  decrease  in 
PH„  which  is  not  observed.  Thus,  although  the  Keicher- 
Meech  model  has  interesting  features,  it  cannot  account  for 
all  of  the  then-available  observations. 

Later,  Bakouh  et  al.  (2006)  confirmed  that  Xenopus 
oocytes  exposed  to  relatively  high  levels  of  NH3/NH4+  (i.e., 
10  mM)  exhibit  the  fall  in  pH,  noted  above,  but  also  found 
that  oocytes  exposed  to  only  0.5  mM  NH3/NH4+  exhibited 
no  change  in  pH,  and  no  substantial  induced  inward  current 
(related  to  net  entry  of  positive  charge,  namely,  NH4  +  ) 
(Bakouh  et  al.  2004,  2006).  The  above  results  are  consistent 
with  the  hypothesis  that  the  NH3/NH4  ^-induced  fall  in  pH, 
and  positive  shift  in  Vm  are  due  to  a  low-affinity  NH4  +  - 
uptake  mechanism  that  is,  to  some  extent,  Cl-  dependent  and 
that  is  largely  inoperative  at  0.5  mM  NH3/NH44 

More  recently,  it  has  become  clear  that  the  AmtB/Rh 
family  of  membrane  proteins  can  serve  as  a  conduit  for 


NH3/NH44.  The  crystal  structures  of  several  prokaryotic 
family  members  are  consistent  with  the  idea  that  NH^ 
passes  through  pores  in  each  of  the  three  monomers  of  the 
homotrimer  (Fabiny  et  al.  1991;  Soupene  et  al.  2002; 
Khademi  et  al.  2004;  Zheng  et  al.  2004;  Andrade  et  al 
2005;  Khademi  and  Stroud  2006;  Conroy  et  al.  2007).  In 
Xenopus  oocytes  expressing  RhCG,  Bakouh  et  al.  (2004, 
2006)  found  that  an  exposure  to  0.5  mM  NH3/NH4  + 
produced  a  rapid,  small,  and  short-lived  alkalinization  that 
was  followed  by  a  slow,  small,  and  sustained  acidification. 
These  authors  proposed  that  the  transient  pH,  increase 
reflects  the  influx  of  NH3,  whereas  the  subsequent  pH, 
decrease  reflects  NH44  influx  in  concert  with  NH*  efflux — 
all  fluxes  mediated  by  RhCG.  However,  as  noted  in  con¬ 
nection  with  the  Keicher-Meech  data,  if  RhCG  could 
mediate  NH^  efflux — and  if  appreciable  NH44  accumu¬ 
lated  inside  the  oocyte  during  the  preceding  NH3/NH4  + 
exposure — then  NH3/NH44  removal  should  have  led  to  a 
large  and  rapid  pH,  decrease,  which  was  not  observed. 

The  purpose  of  the  present  work  was  to  elucidate  the 
movements  of  NH3  vs.  NH44  across  the  plasma  mem¬ 
brane — and  the  potential  role  of  Amt  in  mediating  these 
fluxes — using  the  oocyte  as  a  model  system.  An  ancillary 
goal,  necessary  to  verify  our  interpretation  of  the  data,  was 
to  explore  the  unusual  handling  of  NH3/NH4f  by  the 
oocyte.  Our  approach  was  to  extend  to  NH^  fluxes  a 
technique  that  we  introduced  earlier  to  assess  C02  fluxes 
across  the  oocyte  membrane  (Endeward  et  al  2006;  Musa- 
Aziz  et  al.  2009).  In  the  earlier  work,  we  pushed  a  blunt  pH 
microelectrode  against  the  oocyte  membrane  and  found 
that  an  exposure  to  C02/HC03~  causes  a  transient  rise  in 
surface  pH  (pHs).  First  observed  in  1984  by  De  Hemptinne 
and  Huguenin  (1984),  who  worked  on  rat  soleus  muscle, 
this  rise  in  pH  occurs  as  C02  uptake  causes  the  depletion  of 
C02  at  the  cell  surface,  leading  to  the  reaction 
HCOr+H4  -  C02  +  H20.  In  1986,  Chesler  (1986) 
found  that  exposing  lamprey  neurons  to  the  weak  base  NH3 
causes  a  transient  decrease  in  extracellular  pH.  We 
hypothesized  that  a  hitherto  unseen  flux  of  NH*  into  the 
oocyte  would  deplete  NH3  at  the  extracellular  surface.  This 
depletion  would  lead  to  the  reaction  NH4+  — >  NH3  -F  H4, 
which  would  lower  pHs  and  also  partially  replenish  surface 
NH3  (Fig.  1).  We  found  that  exposures  to  both  5  and 
0.5  mM  elicited  transient  pHs  decreases,  and  that  these 
decreases  were  augmented  by  the  expression  of  AmtB.  The 
exposure  to  5  mM  NH3/NH4+  also  caused  a  paradoxical 
fall  in  pHj;  thus,  in  this  case,  the  pH  fell  on  both  sides  of 
the  membrane.  We  have  developed  a  model  that  allows  us 
to  predict  the  direction  in  which  fluxes  of  NH^  and  NH4  + 
would  affect  pH  on  both  sides  of  the  membrane.  Guided  by 
this  model  and  our  new  pHs  data,  we  conclude  that  (a) 
earlier  models  of  how  20  mM  NH3/NH4f  affects  oocyte 
pH,  are  incomplete;  (b)  in  oocytes  exposed  to  5  or  0.5  mM 
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NH3/NH4\  the  influx  of  NH3  produces  the  dominant 
effects  on  pHs  (though  we  cannot  rule  out  substantial 
NH44  (luxes);  (c)  oocytes  metabolize  or  sequester  incom¬ 
ing  NH3;  and  (d)  AmtB  enhances  permeability  to  NH3 
more  than  that  to  NH44 

Methods 

Expression  in  Xenopus  Oocytes 
cRNA  Synthesis 

As  described  elsewhere  (Musa-Aziz  et  al.  200c>),  we  used 
AmtB  that  we  had  tagged  at  the  C  terminus  with  EGFP 
(enhanced  green  fluorescent  protein).  We  ihen  used  Not  I  to 
linearize  AmtB  cDNA  constructs  in  pGH19,  purified  line¬ 
arized  cDNA  using  the  QIAquick  PCR  purification  kit 
(Qiagen  Inc.,  Valencia,  CA),  transcribed  capped  cRNA 
using  the  T7  niMessage  mMachine  kit  (Ambion,  Austin, 
TX),  and,  finally,  purified  and  concentrated  cRNA  using 
the  RNeasy  MinElute  RNA  Cleanup  Kit  (QIAGEN).  We 
determined  the  RxNA  concentration  using  ultraviolet 
absorbance  and  assessed  quality  using  gel  electrophoresis. 

Xenopus  Oocyte  isolation 

We  surgically  removed  ovaries  from  anesthetized  frogs, 
separated  the  oocytes  using  a  collagenase  treatment 
(Romero  et  al.  1998;  Toye  et  al.  2006),  selected  Stage  V-VI 
oocytes,  and  stored  them  until  use  at  18°C  in  OR3  medium 
supplemented  with  500  U  of  penicillin  and  500  U  of 
streptomycin 


Fig.  1  Model  of  an  oocyle  exposed  lo  NH^/NRC.  The  large  purple 
arrows  indicate  the  direclon  of  bulk  solution  flow.  USL  unslirred 
layer,  Bulk  ECF  bulk  exiracellular  fluid 


Microinjection  of  cRNAs 

One  day  after  isolation,  we  injected  oocytes  with  either 
25  ng  of  cRNA  encoding  AmtB-EGFP  cRNA  (50  nl  of  a 
0.5  ng/nl  cRNA  solution)  or  50  nl  of  sterile  water 
(Ambion,  Austin,  TX)  in  the  case  of  control  (”H20”) 
oocytes.  We  stored  the  oocytes  4-6  days  after  injection  for 
use  in  experiments.  We  verified  delivery  of  EGFP-tagged 
AmtB  to  a  region  near  the  plasma  membrane  using  a  96-well 
plate  reader  (BMG  Labtechnologies,  Inc.,  Durham,  NC)  to 
assess  whole-oocyte  fluorescence  (Toye  et  al.  2006). 

Solutions 

The  nominally  C02/HC03--free  ND96  solution  contained 
(mM):  96  NaCl,  2  KC1,  1  MgCl2,  1 .8  CaCl2,  and  5  HEPES. 
We  titrated  the  solution  to  pH  7.50  using  NaOH  The 
osmolality  of  all  solutions  was  200  mosmol/kg  H20.  We 
made  the  solution  containing  5  mM  NH3/NH4  f  in  ND96 
by  replacing  5  mM  NaCl  with  5  mM  NH4C1.  We  made 
0.5  mM  NH3/NH4  f  in  ND96  by  diluting  the  5  mM 
NH3/NH44  solution  1.10  with  ND96  solution. 

Electrophysiological  Measurements 

Chamber 

Oocytes  were  placed  in  plastic  perfusion  chamber  with  a 
channel  3  mm  wide  x  30  mm  long,  and  constantly  su- 
perfused  at  a  flow  of  3  mi/min.  Perfusing  solutions  were 
delivered  from  plastic  syringes  and  Tygon  tubing  using 
syringe  pumps  (Harvard  Apparatus,  South  Natick,  MA). 
Switching  between  solutions  was  performed  by  pneumati¬ 
cally  operated  valves  (Clippard  Instrument  Laboratory, 
Cincinnati,  OH).  All  experiments  were  performed  at  room 
temperature  (^22°C). 

Measurement  of  Intracellular  pH  (pHj) 

Our  approach  was  similar  to  that  detailed  previously  (Toye 
et  al.  2006;  Lu  et  al.  2006;  Parker  et  al.  2008).  Briefly,  we 
impaled  the  oocyte  with  two  microelectrodes,  one  for 
measuring  membrane  potential  (connected  to  a  model  725 
two-electrode  oocyte  voltage-clamp  amplifier,  Warner 
Instruments  Corp  ,  Hamden,  CT)  and  the  other  for  mea¬ 
suring  pH,  (connected  to  a  model  FD223  high-impedance 
electrometer;  World  Precision  Instruments,  Sarasota,  FL). 
Both  electrodes  had  tip  diameters  of  ~  1  pm.  The  pH 
microelectrode  was  of  the  liquid  membrane  style  (proton 
cocktail  no.  95293;  Fluka  Chemical  Corp.,  Ronkonkoma, 
NY).  The  analog  subtraction  of  the  Vm-electrode  signal 
from  the  pH-electrode  signal  produced  the  voltage  due  to 
pHj.  We  acquired  and  analyzed  data  by  computer,  using 
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software  written  in-house.  We  calibrated  the  pH,/Vm  in  the 
chamber  with  pH  standards  at  pH  6.0  and  8.0  (to  determine 
the  slope),  followed  by  a  single-point  calibration  with  the 
standard  ND96  solution  (pH  7.50)  in  the  chamber  with  the 
oocyte  present,  just  before  impalement.  The  mean  sponta¬ 
neous  initial  Vm  in  this  study  was  —39  ±  2  mV  (n  =  30). 

Measurement  of  Surface  pH  ( pHs ) 

The  microelectrode  for  measuring  pHs  (connected  to  a 
FD223  electrometer;  World  Precision  Instruments)  had  a 
tip  diameter  of  15  pm  and  was  of  the  same  liquid-mem¬ 
brane  design  as  the  pH,  microelectrode.  The  external 
reference  electrode  for  the  Vm  and  pHs  measurements  was 
a  calomel  half-cell  (connected  to  a  model  750  electrometer; 
World  Precision  Instruments)  contacting  a  3  M  KCl-hlled 
micropipette,  which  in  turn  contacted  the  fluid  in  the 
chamber.  The  analog  subtraction  of  the  calomel-electrode 
signal  from  the  pHs-electrode  signal  produced  the  signal 
due  to  pHs.  The  analog  subtraction  of  the  calomel-elec¬ 
trode  signal  from  the  Vm-electrode  signal  produced  the 
signal  due  to  Vni.  The  voltage-clamp  amplifier  generated 
the  virtual  ground  via  a  Ag/AgCl  half  cell  (connected  to  the 
1  Sense  input)  contacting  a  second  3  M  KCl-hlled  micro¬ 
electrode,  the  tip  of  which  we  positioned  close  to  the 
oocyte.  We  used  an  ultrafine  micromanipulator  (model 
MPC-200  system;  Sutter  Instrument  Co.,  Novato,  CA)  to 
position  the  pHs-electrode  tip  at  the  surface  of  the  oocyte, 
and  then  to  advance  it  ^40  pm  further,  whereupon  we 
observed  a  slight  dimple  in  the  membrane.  Periodically,  we 
withdrew  the  electrode  300  pm  from  the  surface  of  the 
oocyte  for  recalibration  in  the  bulk  extracellular  fluid  (pH 
7.50).  The  fluid  bathing  the  oocyte  flowed  at  3  ml/min. 
Relative  to  the  flowing  solution,  the  tip  of  the  pHs 
microelectrode  was  just  in  the  “shadow"  of  the  oocyte 
(Fig.  1). 

Analysis  of  pH  Data 
Initial  clpH/dt 

We  computed  the  initial  rate  of  change  of  pH,  (dpH,/dt) — 
produced  by  introducing  5  m\l  NH3/NH4+ — by  deter¬ 
mining  the  line  of  best  lit  to  the  steepest  part  of  the  pH,  vs. 
time  record. 

Maximum  pHs  Spike  Heights 

We  used  the  following  approach  to  compute  the  maximum 
magnitude  (i.e.,  “spike  height"  or  ApHs)  of  the  pHs  tran¬ 
sient  elicited  by  applying  5  or  0.5  m\l  extracellular  NH3/ 
NH4+.  We  determined  the  initial  pHs — that  is,  before 
application  of  NH3/NH4+ — by  comparing  the  pHs-electrode 


voltage  signal  when  the  electrode  tip  was  at  the  oocyte 
surface  with  the  voltage  signal  obtained  when  the  tip  was  in 
the  bulk  extracellular  fluid  (assumed  pH  =:  7.50).  We 
determined  the  maximum  pHs  during  the  NH3/NH4+ 
exposure  by  comparing  the  voltage  signal  (at  a  time  corre¬ 
sponding  to  the  extreme  pHs  value)  when  the  electrode  tip 
was  at  the  oocyte  surface  with  the  voltage  signal  obtained  a 
few  minutes  later,  when  the  tip  was  in  the  bulk  extracellular 
fluid  (assumed  pH  =  7.50).  The  ApHs  was  the  algebraic 
difference  between  the  extreme  and  the  initial  pH  values, 
one  based  on  a  calibration  in  an  ordinary  ND96  solution  and 
the  other  based  on  a  separate  calibration  in  the  NH3/NH4+- 
containing  solution.  We  performed  separate  calibrations  of 
the  pHs  electrode  in  the  plain  and  NH3/NH4+-containing 
ND96  solutions  to  minimize  potential  errors  in  the  event  that 
NH3/NH4+  affects  the  proton  cocktail. 

NMR  Spectroscopy 

Measurements  of  Metabolites  with  1 H-/1  *C/NMR 

Twenty  oocytes  were  incubated  with  0.5  mM  1SNH/ sNH4+, 
washed  extensively,  and  homogenized  with  ethanol  to  extract 
water-soluble  metabolites.  A  small  quantity  of  [2- nC]glycine 
(30  nmol)  was  added  as  an  internal  concentration  reference 
and  to  control  for  potential  sample  losses  during  the  extraction 
procedure.  After  centrifugation  of  the  homogenate,  the 
supernatant  was  lyophilized  and  resuspended  in  10%  DiO 
buffer  containing  a  small  amount  of  3-(triinethylsilyl)- 
[2,2,3,3-cLi  1-  propionatefNa-4")  (TSP)  as  a  chemical-shift  ref¬ 
erence.  lH-[l  *C]NMR  spectra  were  acquired  fully  relaxed  at 
11.7  T  (Patel  et  al.  2(X)5),  using  a  high-resolution  Bruker 
A vance  NMR  spectrometer  operating  at  500. 13  MHz  ( 1 H) 
and  125.7  MHz  (nC).  Spectral  data  were  acquired  fully 
relaxed  with  a  sweep  width  of  6009  Hz,  8192  data  points, 
interscan  delay  of  20  s,  and  128  scans.  Free-induction  decays 
were  zero  filled,  exponential  filtered  (LB  =  0.5  Hz),  and 
Fourier  transformed.  Metabolite  peak  intensities  were  mea¬ 
sured  by  integration  and  referenced  to  1  C-labeled  glycine 
with  correction  for  differences  in  the  number  of  hydrogen 
atoms.  The  glycine  signal  was  determined  to  be  1(X)%  enri¬ 
ched  with  13C,  indicating  no  significant  overlap  with 
endogenous  signals. 

Measurements  of,5NH/5NH4+  with  HSQC 

Eighty  oocytes  were  incubated  with  I5NH3/15NH4+  for 
30  nun,  followed  by  extensive  washing  with  ISN-free/NH3/ 
NH4^-free  buffer.  The  washed  oocytes  were  cooled  on  ice 
and  homogenized,  yielding  ^0.4  ml  of  supernatant  after 
centrifugation.  The  medium  was  then  acidified  with  50  pi 
of  2.5  M  HC1,  50  pi  of  D20  was  added  for  field-frequency 
lock,  and  samples  were  loaded  into  5-mni  NMR  tubes. 
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NMR  spectra  were  obtained  at  5°C  to  minimize  the 
exchange  of  hydrogen  atoms  between  ammonia  and  water. 
NMR  experiments  were  performed  at  1 1.7  T  using  a  high- 
resolution  Bruker  A vance  NMR  spectrometer  operating  at 
500.13  MHz  (lH)  and  50.683  MHz  (t5N).  Measurements 
of  total  l5NHy15NH4+  concentrations  in  oocytes  were 
determined  with  the  Heteronuclear  (!H-15N)  Single  Quan¬ 
tum  Coherence  (HSQC)  NMR  spectroscopy  technique 
(Grzesiek  and  Bax  1993;  Kanamori  et  al.  1995).  HSQC  is 
well  suited  for  detection  of  nuclei  (e.g.,  l<iN)  with  a  low 
gyromagnetic  ratio  (y).  Because  y  for  ,SN  is  —  10  times 
lower  than  for  !H,  the  sensitivity  for  ammonia-nitrogen 
detection  is  enhanced  by  a  factor  of  (yn/Visr)5  2  ~  306. 
Water  suppression  was  achieved  using  the  Water  Sup¬ 
pression  by  Gradient-Tailored  Excitation  (WATERGATE) 
technique  (Piotto  et  al.  1992).  15 N  decoupling  was  obtained 
using  Waltz  16  (Shaka  et  al.  1983).  Free-induction  decays 
were  acquired  with  a  spectral  width  of  6009  Hz,  8192  data 
points,  interscan  interval  of  10  s,  and  128  scans.  All  spectra 
were  zero  tilled  to  65,536  points.  The  total  ammonium 
(NH3  +  NH4+)  concentrations  in  oocytes  were  determined 
by  comparing  the  integrated  intensities  of  the  ammonia  1 H 


resonance  to  the  incubation  medium  containing  0.5  inM 
isNH3/,?NH44,  using  the  same  pulse  acquisition  condi¬ 
tions,  and  assuming  that  oocytes  are  spheres  of  1.2-mm 
diameter  and  contain  40%  water  by  weight.  lH-‘\N  HSQC 
spectra  of  oocytes  and  their  incubation  medium  were 
acquired  in  a  paired  manner. 

Statistics 

We  present  data  as  mean  ±  SE.  To  compare  the  difference 
between  two  means,  we  performed  a  Student's  /-test  (two 
tails).  To  compare  more  than  two  means,  we  performed  a 
one-way  ANOVA  and  a  Student-Newman-Keuls  (SNK) 
multiple  comparison,  using  KaleidaGraph  (Version  4,  Syn¬ 
ergy  Software).  We  considered  p  <  0.05  to  be  significant. 

Results 

Figure  I  summarizes  how  the  influx  of  NH4  would  affect 
pHs.  As  NH3  enters  the  cell,  [NH4]  near  the  external  sur¬ 
face  of  the  membrane — [NH3|S — falls.  The  replenishment 
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Fig.  2  Intracellular  pH  (pH,)  and  surface  (pHs)  records  during 
application  and  removal  of  5  mM  NH3/NH4+.  a  IHO-injected  oocyte, 
b  AintB  EGFP-tagged  expressing  oocylcs.  At  ihe  indicated  times,  we 
switched  the  extracellular  solution  from  ND%  to  5  mM  NH3/NH4+ 
and  then  back  again,  b  shows  two  representative  experiments:  on  the 
left,  pHs  at  first  decays  rapidly  and  then  more  slowly,  whereas  on  the 
right,  pHs  decays  slowly  throughout  the  NH3/NH4+  exposure.  In  both 
a  and  b  the  purple  record  is  pHs,  the  blue  record  is  pH,,  and  the 
brown  record  is  Vm.  The  vertical  gray  bands  represent  periods  during 


which  the  pHs  electrode  was  moved  to  the  bulk  ECF  for  calibration 
(fixed  pH  of  7.50).  In  H20  oocytes,  mean  pH,  was  7. 10  ±  0.04  before 
and  6.84  ±  0.06  during  exposure  to  5  mM  NH3/NH4+  (n  —  6; 
p  =  0.005),  the  comparable  mean  Vm  values  were  —43  ±  2  and 
—25  ±  2  mV  (n  =  6;  p  —  0.0002).  In  AmtB  oocytes,  the  mean  pH, 
was  7.14  ±  0.04  before  and  6.98  ±  0.06  during  exposure  to  5  mM 
NH3/NH4+  (n  =  8;  p  =  0.037);  the  comparable  mean  Vm  values 
were  -38  ±  4  and  -26  ±  2  </i  =  8;  p  =  0.010) 
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of  the  lost  NH3  near  membrane  surface  occurs  by  two  routes. 
First,  NH3  diffuses  from  the  bulk  extracellular  fluid  to 
approach  the  membrane,  an  effect  that  we  cannot  detect. 
Second,  NH4  +  near  the  membrane  dissociates  to  form  both 
NH3  as  well  as  the  H  f  that  we  detect  as  a  fall  in  pHs.  As  the 
concentration  of  NH3  inside  the  cell  ([NH3],)  rises,  the  influx 
of  NH3  slows,  and  pHs  relaxes  toward  the  pH  of  the  bulk 
extracellular  fluid  (pHBuik)  In  Fig.  1,  we  assume  that  only 
NH3  enters  the  cell.  If  only  NH4+  enters,  the  reverse  set  of 
events  would  occur  and  pHs  would  rise.  In  the  Discussion, 
we  consider  the  situation  in  which  both  NH3  andNH4+  enter 
We  explored  the  model  in  Fig.  1  by  exposing  oocytes  to 
5  or  0.5  111M  NH3/NH4+.  Figures  2  and  3  show  represen¬ 
tative  experiments  at  these  two  levels  of  NH3/NH4+,  and 
Fig.  4  summarizes  the  mean  values  for  these  experiments. 

Effects  of  5  mM  NH3/NH4+ 


Figure  2  shows  the  results  of  representative  experiments  on 
oocytes  exposed  to  a  relatively  high  level  of  NH3/NH4+, 
5  mM. 
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Fig.  3  Intracellular  pH  (pH,)  and  surface  (pHs)  records  during 
application  and  removal  of  0.5  mM  NH3/NH44'.  a  AmlB  EGFP- 
tagged  expressing  oocyte,  b  H20-injecled  oocyle.  At  the  indiealed 
times,  we  switched  the  extracellular  solution  from  ND96  to  0.5  mM 
NH3/NH4+  and  then  back  again.  In  both  a  and  b.  the  purple  record  is 
pHs.  the  blue  record  is  pH,,  and  the  brown  record  is  Vni.  The  vertical 
gray  bands  represent  periods  during  which  the  pHs  electrode  was 
moved  to  lhe  bulk  ECF  for  calibration  (fixed  pH  of  7.50).  In  H20 
oocytes,  the  mean  pH,  was  7.21  ±  0.05  before  and  7.17  db  0.05 
during  exposure  to  0.5  mM  NH3/NH4+  (m  =  8;  p  =  0.61);  the 
comparable  mean  Vm  values  were  — 44  ±  2  and  —40  ±  2  mV 
(n  =  8;  />  =  0,14).  In  AmlB  oocytes,  the  mean  pH,  was 
7  24  ±  0.05  before  and  7.21  ±  0.05  during  exposure  to  0.5  mM 
NH3/NH4+  (n  =  8;  p  =  0.52);  the  comparable  mean  Vm  values  were 
-31  ±4  and  -30  ±  4  (n  =  8;  p  =  0.87) 


HjO -Injected  Control  oocyte 

Figure  2a  refers  to  a  control,  H20-injected  oocyte.  Here — 
as  previously  reported  by  others  for  20  mM  NH3/NH44 
(Burckhardt  and  Fromter  1002;  Keicher  and  Meech 
1094) — the  addition  of  5  mM  NH3/NH4+  to  the  extracel¬ 
lular  fluid  causes  a  moderately  slow  but  sustained  fall  in 
pH,.  vSimultaneously,  pHs  abruptly  falls  by  ^0.2  and  then 
slowly  decays  (i.e.,  rises  toward  pHBuik)-  The  initial  fall  in 
pHs  indicates  that  the  influx  of  NH3  (as  opposed  to  NH4+) 
is  dominant  with  respect  to  effects  on  pHs.  Note  that,  at  a 
time  when  pH;  has  reached  a  stable  value,  pHs  is  still 
substantially  lower  than  pHBulk,  indicating  that  NH3  is  still 
entering  the  cell.  Applying  5  mM  NH3/NH4+also  causes  a 
very  slowly  developing  depolarization  that  reaches  a  peak 
after  nearly  15  min — not  the  instantaneous  depolarization 
expected  of  a  pre-existing  NH4  +  channel — that  slowly 
decays  by  nearly  9  mV. 

The  removal  of  the  NH3/NH4  +  leads  to  a  slow  recovery 
of  pHj  that  is  incomplete  over  the  time  frame  of  this 
experiment  The  removal  of  NH3/NH4  f  causes  pHs  to 
overshoot  pHBuik  by  a  small  amount  (0.048  ±  0.004; 
n  =  5),  and  then  to  decay  slowly.  The  presence  of  an 
overshoot  implies  that — regarding  effects  on  pHs — a  small 
efflux  of  NH3  dominates  over  any  efflux  of  NH4+.  This 
exiting  NH3  presumably  arises  from  a  small  amount  of 
intracellular  NH4+  that  dissociated  to  form  H  (which 
would  lower  pHj)  and  NH3  (which  would  exit  the  cell).  The 
short  duration  of  the  overshoot  implies  that  the  efflux  was 
of  short  duration  (i.e.,  that  the  pool  of  intracellular  NH4^ 
was  small)  If  intracellular  NH4+  had  accumulated  to  a 
substantial  degree  during  the  NH3/NH ^exposure,  then  the 
removal  of  NH3/NH4+  should  have  led  to  an  abrupt  fall  in 
pH,  (which  we  did  not  observe)  and  a  pHs  overshoot  whose 
magnitude  should  have  matched  that  of  the  initial  NH3/ 
NH4+-induced  fall  in  pH$  (i.e.,  ~0.2).  Removing 
NH3/NH4+  also  causes  Vm  to  undershoot  its  initial  level 
by  ~  1  mV  and  then  to  relax  to  the  pre-NH3/NH4+  value. 
Note  that  the  speed  of  the  repolarization  during  NH3/NH4  + 
removal  is  much  faster  than  the  depolarization  during 
NH3 /NH4+  application,  but  still  slower  than  the  speed  of  the 
pHs  increase  during  NH3/NH4+  removal 

Oocyte  Expressing  A/ntB  Tagged  with  EGFP 

Figure  2b  shows  two  experiments  on  oocytes  expressing 
AmtB-EGFP.  At  the  left,  the  switch  to  5  mM  NH3/NH4 
causes  a  fall  in  pHj,  though  one  that  is  slower  than  for  the 
H20  oocyte.  On  the  other  hand,  the  maximal  fall  in  surface 
pH  (which  we  refer  to  as  — ApHs)  is  substantially  larger  for 
the  AmtB-EGFP  oocyte  than  for  the  H20-injected  oocyte, 
implying  that  the  relative  uptake  of  NH3  over  NH4+  is 
higher  in  the  AmtB  oocyte.  In  the  oocyte  at  the  left  in 
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Fig.  2b,  this  shift  in  pHs  at  first  decays  rapidly  and  then 
more  slowly.  In  the  oocyte  at  the  right  in  Fig.  2b,  the  fall  in 
pHs  decays  very  slowly  throughout  the  NHj/NFU*  expo¬ 
sure.  In  both  AmtB  oocytes — as  was  true  for  the  H2O 
oocyte — pH.s  fails  to  decay  fully  to  pHBuik  during  the 
NH1/NH44  exposure.  The  Vm  changes  induced  by  appli¬ 
cation  and  removal  of  NHi/NHa4  are  similar  to,  but 
smaller  than  those  for  H20  oocytes. 

The  removal  of  NH3/NH44  produces  about  the  same 
effects  in  the  AmtB  oocyte  as  in  the  H20  oocyte  (mean 
overshoot  in  AmtB  oocytes,  0.037  ±  0.004;  n  =  6). 

Summary  of  AmtB- EOF P  Versus  H:Q 

As  summarized  on  the  left  side  in  Fig.  4  for  a  larger  number 
of  experiments,  the  effect  of  AmtB  vs.  H20  is  statistically 
significant  for  the  initial  rate  of  pH,  decrease  (dpH/dt) 
elicited  by  NH3/NH44  (Fig.  4a).  The  slower  fall  of  pH,  in  the 
case  of  AmtB  oocytes  is  consistent  with  the  hypothesis  that, 
regarding  effects  on  pH,  (a)  AmtB  acts  by  predominantly 
enhancing  NH-*  vs.  NH4+  influx  (b)  the  enhanced  influx  of 
NH-*  to  some  extent  promotes  the  reaction  NH3  +  H4  -* 
NH4^  in  the  cytosol,  and  (c)  the  consumption  of  cytosolic  H + 
slows  the  NHj/NH4+- induced  fall  in  pH,. 

The  effect  of  expressing  AmtB  did  not  reach  statistical 
significance  regarding  the  NHj/NH^-induced  ApH, 


(Fig.  4b).  However,  AmtB  oocytes  had  a  greater  — ApHs 
(Fig.  4c),  indicating  a  greater  relative  influx  of  NH3  and 
thus  supporting  the  hypothesis  in  the  previous  paragraph. 
Finally,  the  AmtB  oocytes  had  a  smaller  AV,n  (Fig.  4d). 

Effects  of  0.5  mM  NH3/NH4+ 

Figure  3  shows  representative  data  from  a  pair  of  oocytes 
exposed  to  a  relatively  low  level  of  NH3/NH44,  0.5  mM. 

H2O -Injected  Control  Oocyte 

As  shown  in  Fig.  3a  for  a  H20-injected  oocyte,  the  switch 
to  0.5  mM  NH3/NH4  +  causes  a  slow  and  very  slight  fall  in 
pH,.  Although  pHs  abruptly  falls,  — ApHs  is  substantially 
lower  than  for  H20  oocytes  exposed  to  the  much  higher 
level  of  NH3/NH44  (i.e.,  5  mM;  see  Fig.  2a).  Neverthe¬ 
less,  regarding  effects  on  pH$,  the  effect  of  NH3  influx 
(over  NH4+  influx)  is  clearly  dominant.  The  Vn,  changes 
caused  by  application  and  removal  of  0.5  mM  NH3/NH44 
are  much  smaller  than  those  in  the  experiments  with  5  mM 
NH3/NH4+. 

The  removal  of  NH3/NH4+  evoked  virtually  no  recovery 
of  pH,  in  Fig.  3a,  although  some  oocytes  exhibited  a  slight 
pHj  recovery  In  Fig.  3a,  pHs  rose  to  a  value  —0.02  higher 
than  pHBulk.  The  mean  pHs  overshoot  was  0.006  ±  0.003 


Fig.  4  Summary  of  the  mean 
initial  rale  of  pH,  decrease 
(dplt./di)  and  magnitude  of  ihe 
NHr induced  change  in  steady- 
slate  pH,  (ApH,),  ApHs,  and  Vm 
for  larger  groups  of  experiments 
like  those  in  Figs.  2  and  3 
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(//  =  8;  0.006  vs.  0;  p  =  0.28).  These  data  indicate  that — 
regarding  effects  on  pHs — the  efflux  of  Nhh  is  negligible 
relative  to  that  of  NH4^. 

Oocyte  Expressing  AmtB  Tagged  with  EGFP 

Figure  3b  shows  that,  with  an  oocyte  expressing  AmtB- 
EGFP,  the  switch  to  0.5  mM  NH3/NH4+  causes  little 
change  in  pH,.  The  maximal  fall  in  pHs  is  much  larger  for 
the  AmtB-EGFP  oocyte  than  for  the  H20  oocyte,  though 
both  values  are  much  smaller  than  for  their  counterparts  in 
exposures  to  5  mM  NH3/NH4+.  The  Vm  change  is 
minimal. 

The  removal  of  NH3/NH4  +  from  AmtB-EGFP  oocytes 
(Fig.  3b)  evoked  virtually  the  same  response  as  in  H20 
oocytes.  The  mean  pHs  overshoot  was  0.014  ±  0.006 
(//  =  8;  0.014  vs.  0;  p  =  0.30),  indicating  that — regarding 
effects  on  pHs — the  efflux  of  NH^  is  negligible  relative  to 
that  of  NH4+. 

Summary  of  AmtB-EGFP  Versus  H2O 

As  summarized  on  the  right  side  in  Fig.  4,  the  effect  of 
AmtB  vs.  H20  for  an  exposure  to  0.5  mM  NHj/NH^ 
follows  essentially  the  same  pattern  as  for  an  exposure  to 
5  mM,  except  that  the  difference  for  AmtB  vs.  H20  is  not 
statistically  significant  for  dpH./dt  (Fig.  4e).  The  dpH,/dt 
for  H20  oocytes  is  signilicantly  different  from  zero 
ip  =  0.05)  but  not  for  AmtB  oocytes  (/;  =  0.26).  As 
summarized  in  Fig.  4f,  the  ApH;  was  not  different  from 
zero  for  H20  oocytes  (0.038  ±  0.0 10;  //  =  8;  p  =  0.15) 
but  was  significantly — but  not  substantially  different — for 
AmtB  oocytes  (0.040  ±  0.005;  //  =  8;  p  =  0.01).  As  was 
the  case  for  the  exposures  to  5  mM  NH3/NH4+,  the  AmtB 
oocytes  also  had  a  higher  —  ApHs  (Fig.  4g),  indicating  a 

A)  Sham  incubation 


ppm 


Fig.  5  H-NMR  spectra  of  ethanol  extracts  of  oocytes,  a  Sham 
incubation  (i.e.,  oocyics  nol  exposed  lo  NH^/NH4+).  b  Oocytes 
incubated  wiih  0.5  mM  LSNHyISNH4+  for  10  min.  Each  spectrum 
represents  lhe  exiract  of  20  oocytes.  TSP,  3-trimethylsilyl-[2,2,3,3- 
d4J-propionate(Na+))  was  added  as  a  chemical-shift  reference 
(0.0  ppm).  *Glycinc  was  added  during  oocyte  extraction  to  serve  as 
an  inlemal  concentration  slandard  and  lo  control  for  poleniial 
metabolite  losses.  ** Residual  ethanol  contaminant  remaining  after 


greater  relative  NH3  influx.  As  also  was  the  case  for  the 
5  mM  exposures,  the  AVm  was  lower  for  the  AmtB  oocytes 
(Fig.  4h). 

Comparing  ApHs  data  for  5  mM  vs.  0.5  mM  (Fig.  4c 
vs.  g),  we  see  that  the  difference  (always  higher  lor  the 
experiments  with  5  mM  NH}/NH4+)  is  statistically  sig¬ 
nificant  both  for  H20  oocytes  and  for  AmtB  oocytes. 

NMR  Studies  of  Oocytes  Exposed  to  ISN-Labeled 
0.5  mM  NH3/NH4+ 

A  paradox  in  our  results  is  that  AmtB  increases  —ApHs 
(Fig.  4g) — indicating  that  AmtB  enhances  the  influx  of 
NH3 — and  yet  has  no  significant  effect  on  the  initial  dpH,/ 
dt  (Fig.  4e)  or  ApH[  (Fig.  40-  Moreover,  removing  extra¬ 
cellular  0.5  mM  NH}/NH4+  produces  little  pHs  overshoot 
beyond  pHBuik  (Pig-  3a,  b),  indicating  little  efflux  of  NH3. 
The  most  straightforward  explanation  for  these  data  is  that 
the  oocyte  somehow  disposes  of  most  incoming  NH3  by 
either  metabolizing  it  to  a  neutral  product  or  sequestering 
it.  We  test  this  hypothesis  in  the  next  two  sections.  Note 
that,  a  priori,  we  can  rule  out  the  possibility  that  the  oocyte 
metabolizes  or  sequesters  NH4+  (as  opposed  to  NH3).  If  the 
oocyte  first  converted  the  entering  NH^  to  NH44 
(NH3  4-  H  f  -*■  NH4^")  and  then  metabolized  that  NH4+  to 
a  neutral  product(s),  the  result  would  be  an  increase  in  pHj, 
which  we  would  have  easily  observed  in  experiments  with 
0.5  mM  NH3/NH4+. 

Assessment  of  Glutamine  Formation 

The  leading  candidate  lor  the  metabolism  of  NH^  to  a 
neutral  product  would  be  the  conversion  of  glutamate  to 
glutamine  via  glutamine  synthetase:  Glu-  4-  MgATP=  4- 
NH3  -  Gin  4-  MgADP  4-  H2PQ4  . 


B)  0.5  mM  NH3/NH/  incubation 


ppm 


extraction.  Estimated  concentrations  of  metabolites  in  control  nonin- 
jecied  oocytes  not  exposed  lo  NH3/NH4+.  (A)  Glutamate  (Glu): 
6.48  mM.  Estimated  metabolite  concentrations  after  incubation  with 
l,NHyl,NH4+  for  10  min:  Glu:  H20,  1.69  mM.  Glutamine  (Gin): 
1.1  mM  Glu  and  Gin  methylene  protons  (H4  and  H4)  and  lactate 
methyl  protons  (H3)  are  depicted.  Glu  H4  and  H3  were  observed  in 
control  noninjected  oocytes  incubated  with  1<iNH3/l5NH4* 
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To  test  this  hypothesis,  we  incubated  noninjected 
oocytes  for  10  min  either  in  ND96  alone  (for  sham  incu¬ 
bation)  or  in  ND96  plus  0.5  mM  15NHy,sNH4+  We  then 
homogenized  the  oocytes  with  ethanol  to  extract  water- 
soluble  metabolites  and  subjected  the  extract  to  !H-[nC] 
NMR  spectroscopy  (Fig.  5).  The  extract  from  oocytes 
incubated  in  15NH3/15NH4+  exhibited  a  very  low  Gin  sig¬ 
nal  ([Gin]  ^  1.1  mM)  vs.  the  shams  (undetectable  Gin). 
RhAG,  a  component  of  the  erythrocyte  Rh  complex,  has  a 
higher  NH3/CO2  permeability  ratio  than  AmtB  (Endeward 
et  al.  2006;  Musa-Aziz  et  al.  2009),  and  RhAG-expressing 
oocytes  (like  AmtB  oocytes)  exhibit  high  ApHs  values  and 
very  small  pHj  changes  when  exposed  to  0  5  mM  NH3/ 
NH4  +  .  However,  extracts  from  RhAG  oocytes  exposed  to 
0.5  mM  l5NH3/l5NH44  have  undetectable  Gin  (not  shown). 

'H-^N-HSQC  spectra  (acquired  as  in  Fig.  6  from  a 
separate  batch  of  RhAG  oocytes  incubated  with  0.5  mM 
i<'NH4C1)  showed  no  significant  1  N-enrichment  of  amides 
in  the  metabolites,  although  we  only  examined  Gin  and  a 
few  other  likely  metabolites.  Although  we  cannot  rule  out 
the  possibility  that  oocytes  sequester  NH3  via  unknown 
metabolic  pathways,  our  !FT[I3C]NMR  and  Mf-^N-HSQC 

A)  Control  non-injected 

(NH3+NH4*) 

6.5  5.5  4.5  3.5  2. 

ppm 

B)  AmtB 


Fig.  6  Representative  HSQC  spectra  of  Xe nopus  oocytes 

incubated  wilh  i*,NHj/I5NH4+  a  Control  noninjcclcd  oocylcs.  b 
AmtB  oocytes  Total  ammonia  appears  as  a  single  resonance  at 
^6.5  ppm  (relative  to  water,  assigned  to  4.7  ppm)  Spectral  inten¬ 
sities  were  scaled  relative  to  the  intensity  of  0.5  mM  |SN-ammonia  in 
the  incubation  medium  of  each  sample  acquired  with  the  same  pulse 
parameter  settings.  The  residual  water  signal  artifact  appears  al 


data  rule  out  substantial  conversion  of  NH4  to  Gin  and 
provide  no  evidence  for  other  pathways. 

NH3/NH4+  Accumulation  in  Oocyte  Water 

To  determine  if  oocytes  accumulate  NH3/NH4  per  se,  we 
incubated  noninjected  control  oocytes  or  AmtB  oocytes  in 
ND96  containing  0.5  mM  I5NH3/15NH4+  for  30  min  After 
the  incubation,  we  washed  the  oocytes  with  a  NH^-free 
ND96  solution  and,  following  homogenization  and  centri¬ 
fugation,  measured  the  total  1SNH3  +  15NH44  of  the 
supernatant  using  *H-I5N  HSQC  (Fig.  6a,  b).  We  found  that 
15NH3/15NH4"t’ — which  reflects  only  NH3/NH44  taken  up 
during  the  incubation — was  ~  30%  higher  in  AmtB  oocytes 
(4.15  mM  computed  for  total  oocyte  H20)  than  in  control 
noninjected  oocytes  (3.33  mM).  Because  pH,  was  —  7.40  in 
these  experiments — and  pHw  was  7.50  and  [NH3/NH44  ]() 
was  0.5  mM — we  can  conclude  that  if  NH4  full)  equili¬ 
brated  across  the  cell  membrane,  cytosolic  [NHj/NH4+] 
(i.e.,  the  concentration  in  the  water  in  direct  contact  with  the 
inner  surface  of  the  plasma  membrane)  could  only  have  been 
0.5  mM  x  10(7  5O_74O)  =  0.6  mM.  Note  that  even  this 


C)  Summary 

P  =  0.38 


Control  AmtB 
N  on-injected 


4.7  ppm  c  Total  ammonia  (NH3  4-  NH4+)  concentration  in  control 
noninjected  and  in  AmtB  oocytes  Values  arc  mean  ±  SB,  with 
numbers  of  oocyte  samples  for  each  group  in  parentheses.  The 
statistical  comparisons  were  made  using  unpaired  two-tailed  /-tests. 
The  experimental  error  was  large  due  to  fasl  evaporation  of  Ihe 
ammonia  in  both  the  standard  solution  and  ihe  exlracts 
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level  of  cytosolic  NH3/NH4+  would  have  led — upon  the 
removal  of  extracellular  NH^/NKj"1" — to  a  substantial  fall  in 
pHj  and  a  large  transient  overshoot  of  pHs,  neither  of  which 
we  observed  Thus,  we  conclude  that  cytosolic  [NH^/NH^] 
must  have  been  substantially  lower  than  0.6  mM.  The  most 
straightforward  explanation  for  the  extremely  high  level  of 
intracellular  NH^/NH4"  that  we  actually  observed  in  the 
NMR  experiments  (i.e.,  for  the  34  -4  +  mM)  is  that  the  vast 
majority  of  intracellular  NH3/NH.44  was,  in  fact,  confined  to 
an  acidic  intracellular  compartment. 

Discussion 

Overview 

The  pH,  data  in  the  present  paper  confirm  the  earlier  obser¬ 
vation  that  an  exposure  of  a  control  Xenopus  oocyte  (i.e.,  not 
heterologously  expressing  other  proteins)  to  a  relatively  high 
level  of  NH3/NH4+  leads  to  a  paradoxical  intracellular 
acidilication  (Burckhardt  and  Fromter  1992;  Keicher  and 
Meech  1994).  Earlier  investigators  had  also  found  that  an 
exposure  to  a  relatively  low  level  of  NH3/NH44  leads  to  a 
negligible  pH,  change  (Bakouh  et  al.  2004).  We  now  confirm 
both  observations  in  H20  oocytes  (Fig.  4b  and  f,  respec¬ 
tively).  In  addition,  we  make  the  novel  observation  that  the 
expression  of  AmtB  slows  the  pH,  fall  elicited  by  an  exposure 
to  5  mM  NHi/NH4+  (Fig.  4a),  consistent  with  the  hypoth¬ 
esis  that  AmtB  predominantly  promotes  the  influx  of  NH3 
over  NH4V  (regarding  effects  on  pH,).  However,  we  believe 
that  the  most  important  contributions  of  the  present  study  are 
(a)  the  simultaneous  application  of  pHs  and  pH,  measure¬ 
ments.  which  provide  new  insights  into  how  Xenopus  oocytes 
handle  exposures  to  NH3/NH/,  and  (b)  novel  pHs  data  that 
show  that  AmtB  enhances  NH^  fluxes  during  exposures  to 
both  5  mM  and  0.5  mM  NHVNH4+. 

Evidence  for  Unusual  NH3/NH44  Handling 
by  Xenopus  Oocytes 

Regardless  of  whether  the  concentration  is  5  or  0.5  mM, 
exposing  an  oocyte  to  NHj/NH4+  causes  an  initial  pHs 
decrease,  the  explanation  lor  which  is  straightforward  the 
influx  of  NH3  causes  a  depletion  of  NH3  at  the  extracellular 
surface  of  the  oocyte,  leading  to  the  reaction  NH4+  — ► 
NHi  4-  H  1 .  However,  five  other  observations  indicate  that 
the  handling  of  NH^  by  Xenopus  oocytes  is  unusual  com¬ 
pared  to  that  by  other  cells  studied  thus  far. 

First,  when  one  exposes  an  oocyte  to  a  relatively  high 
level  of  NHVNH4+,  pHj  paradoxically  falls.  Almost  all 
other  cells — the  exceptions  being  membranes  that  mediate 
either  a  massive  influx  of  NH44  (Aickin  and  Thomas  1977; 
Kikeri  et  al  19S9)  or  no  apparent  flux  of  either  NH3  or 


NH4+  (Waisbren  et  al.  1994;  Singh  et  al  1995) — exhibit  a 
rapid  pH,  increase,  due  to  the  influx  and  protonation  of 
NH^  (Roos  and  Boron  1981).  Working  on  oocytes,  previ 
ous  investigators  had  observed  NH3/NH4 ^-induced  pH, 
decreases  when  introducing  10-20  mM  NHj/NH^ 
(Burckhardt  and  Fromter  1992;  Keicher  and  Meech  1994; 
Bakouh  et  al.  2006).  In  the  present  paper,  we  observed 
similar,  though  smaller  pH,  decreases — with  both  H20- 
injected  controls  and  AmtB  oocytes — when  introducing 
5  mM  NH}/NH4+  (Fig.  2).  This  paradoxical  fall  in  pH, 
implies  that  the  rate  of  intracellular  acid  loading  is  higher 
than  that  of  acid  extrusion.  We  examine  potential  expla¬ 
nations  below,  under  “Models  of  NH^/NH4+  Handling  by 
Xenopus  Oocytes:  ‘High’  [NH3/NH4+|0." 

Second,  when  one  exposes  an  oocyte  to  a  relatively  low 
level  (e.g.,  0.5  mM)  of  NHi/NHa4,  pH,  either  falls  very 
slowly,  and  by  a  small  amount  (Figs.  3,  4e,  f;  H20  bars),  or 
does  not  change  significantly  (Fig.  4e,  f;  AmtB  bars). 
These  observations  on  H20  oocytes  confirm  earlier  work 
(Bakouh  et  al.  2004,  2006). 

Third,  after  an  exposure  to  relativ  ely  high  levels  of  NH^/ 
NH4+ — 10  to  20  mM  in  the  case  of  others  (Burckhardt  and 
Fromter  1992;  Keicher  and  Meech  1994;  Bakouh  ct  al.  2004, 
2006)  or  5  mM  in  the  present  work  (Fig.  2) — the  removal  of 
extracellular  NH3/NH44  causes  pH,  to  increase  slowl), 
rather  than  to  decrease  rapidly  (due  to  NH^  efflux  followed 
by  the  intracellular  reaction  NH44  — >  NH3  -f  H  1 ),  as  has 
been  observed  by  numerous  investigators  working  on  most 
other  cell  types. 

Fourth,  although  introducing  NH3/NH.44  produces  an 
abrupt  fall  in  pHs — indicative  of  NH3  influx — the  sub¬ 
sequent  decay  in  pHs  is  generally  very  slow,  regardless  of 
whether  the  oocytes  are  exposed  to  5  or  0.5  mM  NH3/NH4+. 
In  H20-injected  oocytes  (Figs.  2a,  3a),  the  decay  is  mono¬ 
tonic  and  extremely  slow  throughout.  In  AmtB  oocytes,  the 
decay  in  pHs  is  sometimes  fast  at  first,  but  then  extremely 
slow  (see  Fig.  2b,  left,  and  Fig.  3b).  In  some  other  experi¬ 
ments  on  AmtB  oocytes,  even  this  rapid  phase  of  pHs  decay 
is  absent  (Fig.  2b,  right).  In  neither  case  does  pHs  return  to 
the  pHBu!k  value  of  7.50  during  the  course  of  an  NH3/NH4 1 
exposure  lasting  ~~  15  min  or  more.  Thus,  the  decay  in  pHs 
in  these  experiments  with  NH3  is  decidedly  slower  than  the 
rapid  decay  in  pHs  in  oocytes  exposed  to  C02  (Endeward 
et  al.  2006) — a  decay  that  is  complete  in  —  5  min  These 
observations  are  consistent  with  the  hypothesis  that — unlike 
C02  (which  equilibrates  across  the  membrane  in  a  few 
minutes) — NH3  is  still  far  from  equilibrated  across  the 
plasma  membrane  and  continues  to  enter  the  oocyte  for  a 
prolonged  period  (Figs.  2,  3).  A  corollary  is  that  the  oocyte 
somehow  maintains  a  relatively  low  level  of  NH*  near  the 
inner  surface  of  the  cell  membrane. 

Fifth,  although  one  would  expect  that  removing  extra¬ 
cellular  NH}/NH44  would  cause  pHs  to  increase  abruptly 


Springer 


R.  Musa-Aziz  cl  al. :  NH3  Permeability  of  oocytes 


25 


and  overshoot  the  initial  pHs  (i.e.,  near  the  pHBulk  of 
7.50) — by  an  amount  that  approximates  the  magnitude  of 
the  pHs  decay  during  the  preceding  NH3/NH4+  expo¬ 
sure — in  fact  we  observe  little  or  no  pHs  overshoot 
(Figs.  2,  3).  Moreover,  the  expression  of  AmtB  (vs.  the 
injection  of  HFO)  increases  the  magnitudes  of  the  rapid 
pHs  changes  but  does  not  produce  a  larger  overshoot 
These  results  imply  that  the  NH3  efflux  from  the  oocyte  is 
small,  consistent  with  the  hypothesis  that  (1)  the  pathways 
for  NH3  movement  across  the  oocyte  membrane  are  highly 
rectified  (i.e.,  favoring  influx  over  efflux),  or  (2)  very  little 
NH44  is  present  in  the  cytosol  near  the  plasma  membrane 
at  the  end  of  an  NH3/NH4+  exposure.  Thus,  upon  removal 
of  extracellular  NH3/NH4+,  very  little  free,  cytosolic  NH44 
is  available  for  the  reaction  NH4+  -►  H+  +  NH3.  This 
reaction  is  necessary  not  only  for  the  pH|  to  fall,  but  also  to 
generate  the  NH3  that  subsequently  exits  and  produces  the 
pHs  overshoot.  Moreover,  a  limited  cytosolic  would 

minimize  the  outwardly  directed  diffusion  potential  for 
NH4+,  and  minimize  a  rebound  hyperpolarization  upon 
removal  of  extracellular  NH3/NH4+.  Because  we  know  of 
no  examples  of  a  rectifying  of  a  gas  flux  through  either 
lipid  or  a  protein  channel,  option  2  is  the  simplest  expla¬ 
nation,  and  would  account  for  the  observations  that 
removal  of  extracellular  NH3/NH4+  fails  to  cause  (a)  a  fall 
in  oocyte  pH*  (b)  a  large  pHs  overshoot,  and  (c)  a  rebound 
hyperpolarization.  Instead,  NH3/NH4+  removal  merely 
eliminates  the  processes  that  cause  pHs  to  fall  and  Vm  to 
become  more  positive,  causing  these  parameters  simply  to 
return  to  their  pre-NH3/NH4f  values. 


Models  of  NH3/NH4+  Handling  by  Xenopus  Oocytes: 
“Low”  [NH3/NH4+]0 

In  the  section  following  this  one,  we  consider  potential 
explanations  for  the  pH,  and  pHs  data  obtained  with  the 
exposure  to  and  withdrawal  of  relatively  high  NH3/NH4  + 
levels  (e.g.,  5  rnM),  with  their  attendant  large  and  paradox¬ 
ical  pH;  changes.  In  this  section,  we  examine  experiments 
with  relatively  low  NH3/NH4  +  levels  (e.g.,  0.5  mM),  where 
the  paradoxical  pH,  changes — mediated  by  a  process  with  a 
low  affinity  for  NH3  and/or  NH4+ — are  minimal.  We  con¬ 
sider  five  models  of  NH3  handling  by  Xenopus  oocytes 
(Fig.  7)  that  could  account  for  why  an  exposure  to 
NH3/NH4+  would  cause  virtually  no  change  in  pH,. 


Model  A;  Is  the  pHj  Effect  of  NHj  Influx  Perfectly  Matched 
by  the  pHj  Effect  of  NH 4 4  Influx  (Fig.  7a)? 


In  order  for  the  influxes  of  NH3  and  NH44  to  cause  no 
change  in  pH,,  the  influx  of  NH3  (7nh,)  and  the  influx  of 
NH4+  (7nu4  )  would  have  to  be  in  the  ratio  I0(pll‘  pKa), 
where  pH,  is  the  initial  pH,  and  pKa  refers  to  the  intra¬ 
cellular  equilibrium  NH4+  «=*  NH3+H+.  To  make  this 
point  more  clearly,  we  imagine  that  the  oocyte’s  cytosol 
contains  a  tiny  amount  of  NH3/NH4+  that  is  equilibrated  at 
the  initial  pH,.  If  we  assume  that  the  pKa  of  this  reaction  is 
9.2,  and,  for  the  sake  of  simplicity,  that  the  pH,  is  7.2,  then 


iop».  pk  =  io7-2-9-2 

[nh4+], 


100 


(0 


Fig.  7  Models  of  NH>/NH4+ 
handling  by  Xenopus  oocytes 
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Thus,  if  7nh,/7n h4-  also  were  1/100,  then  the  parallel 
influxes  of  NH3  and  NH4+  would  not  disturb  the  ratio 
|NH3]i/[NH4+]i  and  hence  would  not  disturb  the 
NH3/NH4  +  equilibrium  or  pH;.  (In  Fig.  7a,  we  assume  an 
NH3  influx  of  100  arbitrary  units.)  We  now  define  the 
virtual  ratio  that  would  produce  no  change  in 

the  actual  pH, 


JN'H>  A  [NH.l],  _  .QpH.-pK. 

Nun  [NHj  *  ]j 


(2) 


In  our  example  (7nh./^Hj*)nuii  *s  1/100.  To  be  more 
general,  we  can  dehne  (7mi*/7NH4-)nuIi  in  terms  of  the  pH 
on  either  side  of  the  membrane: 


_Jnir_\  _  [NHjji  _  | QPi i  pK„ 
JNH.  /Null  [NH4  +  ],  ' 


(3) 


Because  solutions  on  opposite  sides  of  the  membrane 
will  generally  have  different  pH  values,  they  will  also  have 
different  requirements  for  (7nh,/7nh4* )nuu»  which — as  we 
see  below — can  have  interesting  consequences. 

We  can  also  define  the  virtual  value  of  pH  that — given  the 
actual  fluxes  of  NH3  and  NH4+ — would  produce  no  pH 
change  in  the  compartment  under  consideration.  Starting  with 
an  expression  analogous  to  Eq.  3  and  solving  for  pH,  we  have 

pHNun  =  pKa  +  log-^  (4) 

7nH4~ 


This  equation  has  the  form  of  the  familiar  Henderson- 
Hasselbalch  equation,  except  that  here  we  replace 
concentrations  with  fluxes.  If  the  Mir  ratio  were  1/ 
100  and  pKa  were  9.2,  then  pHNun  would  be 


PHnuh  =  9.2  +  logjj^Q  =  7.2  (5) 

If  the  actual  pH,  were  7.2  (i.e.,  the  same  as  pHNun  in  this 
example),  the  NH3/NH4f  fluxes  would  not  alter  pH,. 
However,  because  the  pH  of  the  bulk  extracellular  fluid  is 
7.5,  this  same  ratio  °f  1/100  would  necessarily 

alter  pHs.  To  determine  the  direction  of  the  effect  on  pHs, 
we  return  to  Eq.  3.  For  the  NH3/NH44  fluxes  to  produce 
no  change  in  pHs  (where  the  initial  pHs  =  pHBulk  =  7.50), 


I  Qfr^Bulk 


pK.)  _  1  q(7.5  —9.2)  _  J_ 

50 


To  produce  the  observed  decrease  in  pHs,  as  in  Fig.  1 ,  the 
fluxes  would  have  to  be  in  a  ratio  >1/50  (or  2/100). 
However,  if  7mi,/7mv  were  only  1/100 — the  value 
necessary  to  stabilize  pH) — the  influx  of  NH3  would  be  too 
low  (or  the  influx  of  NH4+  would  be  too  high)  to  stabilize 
pHs,  and  the  following  reaction  on  the  extracellular  surface 
of  the  oocyte  would  replenish  some  of  the  lost  NH4+: 
NH,  -f  H+  ->  NH4+  (Fig.  7a).  Thus,  if  the  7Ni h/7Nn4- 
ratio  were  positioned  to  stabilize  pH,,  pHs  would  rise. 


rather  than  fall  as  actually  observed.  Obviously, ./ nii,/7nh4* 
cannot  simultaneously  be  1/100  (to  explain  the  stability  of 
pH,)  and  >2/100  (to  explain  the  fall  in  pHs).  Therefore,  no 
combination  of  NH3  influx  and  NH4+  influx — regardless  of 
mechanism  (e  g.,  influx  through  a  Na/K/Cl  cotransporter  or 
a  channel) — can  simultaneously  account,  by  itself,  for  the 
pH,  and  pHs  data.  Thus,  we  can  rule  out  model  A  for 
experiments  at  0.5  mM  NH3/NH4  +  Below,  we  see  that  a 
similar  line  of  reasoning  rules  out  NH4+  influx  as  a 
potential  mechanism  for  the  paradoxical  fall  of  pH,  in 
5  mM  NH3/NH4+. 

Table  I  summarizes  the  directions  of  the  expected  pH 
changes  for  the  three  possible  conditions:  pH  <  pHNuu, 
pH  =  pHNun,  and  pH  >  pH^,.1  These  pH  vs.  pHNun 
relationships  are  mirrored  by  corresponding  JswJJsiu*  vs- 
(7nh3A/nii,^)nuii  relationships.  Notice  that,  in  the  case  of 
inequalities,  the  directions  of  the  pH  changes  are  opposite 
on  the  side  from  which  the  NH3  and  NH4+  exit  vs  the  side 
which  NH3  and  NH4^  enter. 

Note  that  the  model  that  we  developed  here  only  allows 
us  to  conclude  that  a  NH3/NH4+-induced  fall  in  pHs  in  our 
experiments  is  associated  with  a  7nh,/7nh4-  that  is  >1/50. 
We  cannot  rule  out  an  influx  of  NH4+.  In  fact,  a  hypo¬ 
thetical  influx  of  NH4’f  could  be  49-fold  greater  than  the 
influx  of  NH3,  and  yet  pHs  still  would  fall.  On  the  other 
hand,  a  large  absolute  influx  of  NH4+  would  slow  the  fall 
in  pHs  and  reduce  the  magnitude  of  -pHs.  Such  kinetic 
issues  can  only  be  addressed  by  a  model  that  computes  the 
time  course  of  pHs. 

Model  B:  Is  the  Influx  of  NHj  Perfectly  Matched 
by  an  Influx  of  H*  (Fig.  7b)? 

Even  if  NH4^  did  not  enter  the  oocyte,  pH,  would  be  stable 
if  7m  1,  were  matched  by  a  comparable  7114_.  Each  entering 
NH3  would  undergo  the  reaction  NH3  -f  H  +  -►  NH4+,  and 
the  lost  cytosolic  H+  would  be  replenished  by  an  influx  of 
H  In  Fig.  7b,  we  assume  a  7n+  of  100  arbitrary  units  and 
a  7mi,  of  100,  and  we  also  assume  that  the  NH3  disap¬ 
pearing  from  the  extracellular  surface  of  the  cell  would  be 
replenished  by  an  NH3  diffusion  from  the  bulk  extracellular 
fluid  (ECF)  of  90  and  a  contribution  of  10  from  the 
extracellular  reaction  NH4+  — *  H+  4-  NH3.  Because  this 
model  results  in  a  net  consumption  of  H4  at  the  extracel¬ 
lular  surface  (100  -  10  =  90  in  this  example),  pHs  would 
rise  rather  than  fall  as  we  observe.  Moreover,  this  model 
would  predict  an  enormous  production  of  intracellular 
NH4+,  which  would  produce  cell  swelling;  we  observe  no 


1  In  our  analysis,  we  assume  that  NH3  and  Nhb"*"  move  in  the  same 
direction.  If  they  should  move  in  opposite  directions,  then  pH  would 
always  fall  on  the  side  of  the  membrane  toward  which  NH4+  moves, 
and  would  always  rise  on  the  opposite  side. 
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Table  1  Predicied  pH  changes 

Condition  Effect  on  compartment  from  Effect  on  compartment 

which  NH3  and  NH3/NH4  leave  lhai  NH3  and  NH3/NH4+  enler 


pH  <  pHNut,  ., 

..  lhai  is,  ./nh,/7nh4  * 

>  (7nh,/7nh40nuh 

1  pH 

t  pH 

pH  =  pHNu]]  . 

. .  thal  is,  7nHi/7jsjh4’ 

=  O/nHiA/nhu  ‘  )nuI1 

£pH 

A  PH 

pH  >  pHNui| 

lhai  is,  4 

<  (W7niv)nuii 

T  pH 

l  pH 

Note.  |pH,  decrease  in  pH;  |pH,  increase  in  pH;  ^  pH,  no  pH  change  in  this  comparlmcni 


tendency  for  the  oocyte  to  swell.  Finally,  the  massive 
accumulation  of  intracellular  NH44  would,  upon  removal 
of  extracellular  NH3/NH44,  lead  to  a  substantial  decline  in 
pH,  and  overshoot  of  pHs  (i.e.,  increase  beyond  the  pHBuik 
value  of  7.50),  neither  of  which  we  observe.  In  this  anal¬ 
ysis,  we  have  assumed  no  flux  of  NH4  + .  To  the  extent  that 
NH44  entered  the  oocyte,  we  would  require  less  H4  influx 
to  produce  no  pHj  change.  However,  any  substantial  H^ 
influx  would  be  inconsistent  with  our  pHs  data,  and  we 
have  already  seen  in  our  analysis  of  model  A  that  NH4+ 
influx  cannot  account  for  our  data.  It  is  theoretically  pos¬ 
sible  that  the  oocyte  could  export  NH4T  thereby  avoiding 
swelling.  However,  the  electrochemical  gradient  for  NH4+ 
is  inward,  and  in  any  case  pHs  would  rise  (rather  than  fall, 
as  we  observe).  Therefore,  model  B  is  incorrect. 

Model  C:  Is  the  Influx  of  NH3  Perfectly  Matched 
by  the  Metabolic  Production  of  /f4  (Fig.  7c)? 

This  analysis  is  similar  to  that  for  model  B  except  that  we 
can  ignore  the  incorrect  predictions  that  stem  from  IT 
influx  across  the  plasma  membrane.  Nevertheless,  we  are 
left  with  the  massive  accumulation  of  intracellular  NH44, 
which  would  lead  to  cell  swelling  and— upon  NH3/NH4  + 
removal — a  large  pH*  decline  and  a  large  pHs  overshoot, 
neither  of  which  we  observe. 

Model  D:  Is  Virtually  all  Entering  NH3  Metabolized 
to  a  Neutral  Product(s)?  (Fig.  7d)? 

If  an  amido-transferase  reaction  (e.g.,  the  conversion  of 
glutamate  to  glutamine)  consumed  the  entering  NH3,  pH; 
would  not  change.  Moreover,  NH44  would  not  accumulate 
inside  the  cell  during  the  NH3/NH4+  exposure.  Thus  fol¬ 
lowing  the  withdrawal  of  extracellular  NH3/NH44",  pH, 
would  not  fall  and  pHs  would  not  overshoot  pHBulk.  To  test 
model  D  we  attempted  to  detect  Gin  by  both  !H-[13C] 
(Fig.  5)  and  'H-^N-HSQC  NMR  spectroscopy,  but 
observed  little  Gin  or  other  likely  NH3  metabolite.  These 
data,  combined  with  our  demonstration  of  substantial 
NH3/NH4+  accumulation  (Fig.  6) — presumably  in  acidic 
intracellular  vesicles — make  it  unlikely  that  the  oocyte  per¬ 
forms  substantial  conversion  of  NH3  by  metabolic  processes. 


Model  E:  Is  Virtually  all  Entering  NHj  Sequestered 
in  an  Acidic  Subcotnparttnent  as  Nil 4  T  (Fig.  7e )  ? 

In  this  model,  a  small  fraction  of  entering  NH3  equilibrates 
with  H+  to  produce  NH44  in  a  cytoplasmic  microenvi¬ 
ronment  immediately  below  the  surface  of  the  plasma 
membrane.  Both  the  NH3  and  the  NH/  would  diffuse 
slightly  deeper  into  the  oocyte,  where  NH3  enters  an  acidic 
vesicle  and  becomes  trapped  as  NH44.  To  test  this 
hypothesis,  we  used  HSQC  NMR  spectroscopy  to 

measure  total  intercellular  NH3/NH44,  finding  values 
(Fig.  6)  that  were  far  too  high  to  represent  cytosolic  NH3/ 
NH44.  Thus,  we  conclude  that  the  vast  majority  of  intra¬ 
cellular  NH3/NH4+  must  be  trapped  as  NH4+  in  vesicles. 
Indeed  oocytes  contain  copious  yolk  granules  or  platelets 
(50%  of  oocyte  volume),  containing  yolk  proteins  (80%  of 
total  cell  proteins),  and  having  a  pH  of  ~5.6  (Fagotto  and 
Maxfield  1994a,  1994b;  Fagotto  1995).  By  microscopy,  a 
layer  of  vesicles  begins  within  ~  10  pm  of  the  plasma 
membrane  (Lu  et  al.  2006). 

We  predict  that — at  least  during  brief  experiments  such 
as  ours — NH3/NH4^  levels  would  rise  only  modestly  in 
the  small  subcompartment  between  the  membrane  and  the 
vesicles,  and  relatively  little  of  the  entering  NH3  could 
escape  the  aforementioned  vesicles  to  penetrate  deeper  into 
the  oocyte.  The  influx  of  NH3  into  the  vesicle  would  lead  to 
a  fall  in  [NH3]  at  the  vesicle  surface,  favoring  the  reaction 
NH44  -►  NH3  4-  H  ,  which  would  minimize  pH,  changes. 
It  is  possible  that  a  high  intravesicular  buffering  power 
(provided  by  abundant  yolk  proteins)  could,  by  itself, 
sufficiently  stabilize  intravesicular  pH  over  the  course  of 
our  experiments.  In  addition,  vesicular  lT  pumps  could 
replenish  intravesicular  H4,  with  cytosolic  metabolism 
providing  the  necessary  H4,  again  tending  to  stabilize 
cytosolic  pH 

The  NH^-trapping  hypothesis  would  account  for  all 
key,  paradoxical  findings  dealing  with  NHV  (1)  The 
extracellular  addition  of  0.5  nixM  NH3/NH44  produces  a 
pHj  change  of  virtually  nil  (the  conversion  of  incoming 
NH3  to  NH4+  occurs  in  a  subcompartment,  not  in  the 
cytosol).  (2)  The  fall  in  pHs  produced  by  the  extracellular 
addition  of  NH3/NH44  relaxes  very  slowly  (NH3  continues 
to  enter  the  subcompartment,  and  perhaps  deeper  into  the 
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oocyte,  for  many  tens  of  minutes).  (3)  The  extracellular 
removal  of  NH3/NH4+  causes  virtually  no  fall  in  pHj 
(virtually  no  NH4+  is  present  in  the  cytosol).  And  (4)  the 
extracellular  removal  of  NHj/NH**  produces  little  or  no 
ov  ershoot  of  pHs  (because  little  NH3/NH44  is  present  free 
in  the  cytosol  near  the  plasma  membrane,  the  efflux  of  NH3 
is  very  low). 

Models  of  NH3/NH4+  Handling  by  Xenopus  Oocytes: 

‘  High”  [NH,/NH4+]0 

Compared  to  the  analysis  of  data  from  experiments  with 
0.5  mM  NH3/NH4+,  that  of  data  from  experiments  with 
5  111M  NH1/NH4+  is  complicated  by  the  paradoxical  fall  in 
pH,  (Fig.  2),  which  creates  the  added  paradox  that  pH  falls 
on  both  sides  of  the  membrane. 

Model  A:  Does  the  Influx  of  NH4  4  Cause  the  Paradoxical 
Fall  in  pH,? 

Previous  investigators  favored  the  hypothesis  that  the 
paradoxical  fall  in  pH,  is  due  to  the  influx  of  NH44 — via 
either  an  Na/K/CI  cotransporter  (Keicher  and  Meech  1994) 
or  a  nonselective  cation  channel  (Burckhardt  and  Fromter 
1992) — followed  by  the  cytosolic  reaction  NH4+  — > 
NH3  +  H+  The  channel  hypothesis,  in  particular,  requires 
that  the  oocyte  dispose  of  the  newly  formed  NH3 — by 
either  NH^  efflux  or  intracellular  NH^  metabolism  or 
sequestration.  Our  pHs  data  (see  Fig.  2)  demonstrate  a  net 
movement  of  NH3  into  the  cell,  ruling  out  the  first  NH3- 
disposal  option.  In  the  introductory  section,  we  noted 
limitations  of  the  Na/K/Cl-cotransporter  hypothesis. 
Moreover,  although  we  agree  that  a  NH44  conductance  is  a 
reasonable  explanation  for  the  depolarization  triggered  by 
the  exposure  to  NH3/NH44,  two  arguments  will  lead  us  to 
conclude  here  that  the  hypothesized  NH4  +  entry  via  a 
nonselective  cation  channel  is  not  a  viable  explanation  for 
the  paradoxical  fall  in  pHj.  When  we  apply  5  mM  NH*/ 
NH44,  the  speed  of  the  pHj  decrease  (an  index  of  the 
hypothesized  NH4+  influx)  is  maximal  early  on  and  then 
gradually  wanes  as  pH,  stabilizes  (see  Fig.  2).  On  the  other 
hand,  the  positive  shift  in  Vm  is  small  early  on  and  slowly 
reaches  a  maximal  value  over  nearly  15  min  Thus,  if  the 
depolarization  is  an  index  of  permeability  to  NH4+,  we 
conclude  that  the  NH4+  conductance  develops  far  too 
slowly  to  account  for  the  decrease  in  pH,. 

Finally,  using  the  same  logic  as  we  did  for  model  A  in 
the  previous  section  (see  Fig.  7a),  we  can  conclude  that  it  is 
impossible  for  the  NH3/NH4+-influx  ratio  to  be — at  the 

We  cannot  rule  out  the  possibility  that  the  NH4+  conductance  is 
immediately  high,  but  that  other  conductances — also  initially  high — 
gradually  decline  to  allow  Vm  to  approach  Enu,-  . 


same  time — low  enough  to  cause  pH,  to  fall  and  yet  high 
enough  to  cause  pHs  to  fall.  For  example,  if  pH,  is  7.2  and 
pKa  is  9.2,  then  UunflJt<nA~)>tuii  would  be  1/100.  Thus,  for 
the  influx  of  NH44  to  lower  pH,-,  Jun  would  have  to 

be  <1/100.  However,  as  we  saw  earlier,  for  the  influxes  of 
NH3/NH4  +  to  lower  pHs  (7nii,A/nh4-)nu11  would  have  to 
be  >2/100.  Because  (J^ h,/^nh4*)nuII  cannot  simultaneously 
be  <1/100  and  >2/100,  the  influx  of  NH4+ — regardless  of 
mechanism — cannot  account  for  the  paradoxical  fall  in 
pH,. 

We  might  note  that,  at  least  in  theory,  it  would  be 
possible  for  the  parallel  influxes  of  NH3  and  NH44  to  cause 
both  pH,  and  pHs  to  fall.  As  noted  in  Table  I,  in  the 
compartment  that  NH^  and  NH44  enter,  the  pH  would  fall 
if  pH  >  pHNuU.  As  we  have  just  seen,  a  fall  in  pHs  from  an 
initial  value  of  7.5  requires  that  7nHj/^mu*  be  >2/100;  let 
us  assume  a  of  10/100.  For  this  ratio,  pHNul| 

would  be  8.2.  Thus,  if  the  initial  pH,  were  >8.2,  a  J mi  / 
7mh4-  of  10/100  would  cause  pH,  to  fall  along  with  pHs. 
(Of  course,  because  the  actual  initial  pHj  is  far  less  than 
8.2,  this  explanation  is  not  valid  for  our  data.)  Using  sim¬ 
ilar  logic,  we  could  predict  conditions  in  which  the  parallel 
influxes  of  NH3  and  NH4+  would  raise  both  pH,  and  pHs. 

Model  B:  Does  the  Influx  of  H+  Causes  the  Paradoxical 
Fall  in  pH,? 

As  noted  in  the  previous  section’s  model  B,  in  our  analysis 
of  data  for  0.5  mM  NH3/NH4+,  even  an  H  inllux  suffi¬ 
cient  to  stabilize  pH,  in  the  face  of  an  NH3  influx  (see 
Fig.  7b)  would  lead  to  a  rise — rather  than  the  observed 
fall — in  pHs.  Because  an  even  greater  H1  influx  would  be 
needed  to  produce  a  net  fall  in  pH, — and  such  a  greater  H" 
influx  would  cause  an  even  greater  increase  in  pHs — we 
can  rule  out  the  H+-influx  model. 

Model  C:  Does  the  Intracellular  Release  or  Production 
of  H+  Cause  the  Paradoxical  Fall  in  pH,? 

As  noted  in  the  previous  section’s  model  C  (see  Fig.  7c), 
the  intracellular  generation  of  HT — by  itself — during  the 
NH3/NH4+  exposure  would  lead  to  NH44  accumulation  in 
the  cytosol  and  thus  cell  swelling  (not  observed).  During 
NH3/NH4+  withdrawal,  the  accumulated  NH4+  would 
dissociate  to  produce  H+  (we  observed  no  fall  in  pH,)  and 
NH3,  which  would  exit  the  cell  (we  observed  no  substantial 
pHs  overshoot).  However,  if  H  '  generation  occurred  in 
parallel  with  trapping — in  acidic  intracellular  vesicles — of 
nearly  all  incoming  NH3,  then  cytosolic  H4  production 
would  promote  accumulation  of  NH44  in  vesicles  rather 
than  in  the  cytosol.  Thus,  this  hybrid  H ^-generation/ 
vesicular  NH44-trapping  model  (see  Fig.  7e)  would 
account  for  our  data  with  5  mM  NH3/NH4+. 
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Model  D:  Does  the  Closing  of  NH -Permeable  Channels 
Cause  pHl  to  Rise  Following  Withdrawal  of  High  Levels 
ofNHj/NHf? 

The  only  explanation  offered  by  previous  investigators  for 
the  observed  rise  in  pH,  with  NH3/NH4  +  removal  is  that 
the  sudden  repolarization  of  the  oocyte  membrane  would 
close  nonselective  cation  channels  and  thus  allow  a  slow 
efflux  of  accumulated  cytosolic  NH4+,  leading  to  a  slug¬ 
gish  pH,  recovery  (Burckhardt  and  Fromter  1992). 
Presumably,  the  release  of  previously  sequestered  NH3 
would  lead  to  the  cytosolic  reaction  NH3  +  HT  NH4+, 
which  would  lead  to  a  slow  rise  in  pH;.  One  argument 
against  the  channel  model  is  that,  as  already  noted,  our 
data'  indicate  that  [NH44]  in  the  cytosol  near  the  inner 
surface  of  the  plasma  membrane  must  be  very  low,  and 
thus  NH4f  efflux  could  not  produce  a  sustained  pH, 
increase.  Consistent  with  this  idea,  following  the  removal  of 
NH3/NH4\  the  Vm  undershoot  is  small  and  short-lived 
(Fig.  2).4  A  much  stronger  argument  flows  from  our  new 
analysis  of  (J nh/^nh4-  )nuii  values.  For  a  hypothetical  efflux 
of  NH4+  to  cause  a  rise  in  pHj — starting  from  an  initial  pH, 
of,  say,  6.9 — the  absolute  value  of  (./nh,A/nii4-)nuii  would 
have  to  be  <[NH1  ]/[ NH4^],  =  10<6'9  92)  =  0.5/100  (see 
Eq  3  and  Table  1).  On  the  other  hand,  even  though  pHs 
does  not  exhibit  a  substantial  overshoot  of  its  initial  value, 
the  small  pHs  overshoot  that  we  do  observe  indicates  that 
would  have  to  be  >10(7  5  9,2)  =  2/100  (see 
Table  1).  Because  ^  cannot  simultaneously 

be  <0.5/ 1 00  and  >2/100,  we  can  conclude  that  an  NH44 
efflux,  regardless  of  mechanism,  cannot  account  for  the 
slow  pH,  increase. 

Model  E:  Might  Endogenous  Na-H  Exchange  Cause  pHt  to 
Rise  Following  Withdrawal  of  High  Levels  of  NH  fNllC? 

After  an  intracellular  acid  load  induced  by  exposure  to  CO2 
or  butyric  acid,  Xenopus  oocytes  not  heterologously 
expressing  acid-base  transporters  have  very  low  rates  of 
acid  extrusion  (Romero  et  al.  1997;  Grichtchenko  et  al. 
2001;  Piermarini  et  al.  2007).  Thus,  the  pH,  recovery 
occurs  only  after  removal  of  NH^/NH^4.  Furthermore,  by 
analogy  to  the  argument  made  in  the  previous  section’s 

'  The  removal  of  NH3/NH4't’  causes  neither  a  fall  in  pH,  nor  a  rise  in 
pHs  lhat  subslaniially  overshoots  lhe  initial  value. 

4  Assuming  that  lhe  NH4+-eondueiive  pathway  remained  activate 
and  that  substantial  NH4'1'  were  present  in  lhe  cytosol,  lhe  removal  of 
extracellular  NH4+  would  create  a  diffusion  potential  that  would  drive 
Vm  10  well  below'  the  prc-NH3/NH4+  value.  Instead,  we  observed  a 
Vm  minimal  undershoot  that  decayed  rapidly,  presumably  reflecting 
either  NH4+  efflux  per  se,  or  NH^  efflux  followed  by  the  cytosolic 
reaction  NH4+  NIT  +  HT  We  have  no  data  that  bear  on  lhe 
decay  of  the  presumed  NH4^  conductance. 


point  B  (see  Fig.  7b),  the  extrusion  of  H4  would  have 
produced  a  fall  in  pHs,  rather  than  the  small  overshoot  that 
we  observed  (Fig.  2).  Thus,  acid  extrusion  by  an  endoge¬ 
nous  transporter  also  cannot  account  for  the  pH,  recovery. 

Model  F:  Might  the  Activation  (or  Inactivation)  of 
Cytosolic  7/ 4  Production  Cause  the  Fall  (or  rise) 
in  pH i  Caused  by  the  Application  (or  Removal)  of 
Extracellular  NH3/NH4+  ? 

Although  we  can  rule  out  NH4+  or  H4  uptake  as  an 
explanation  for  the  paradoxical  pHj  decrease  caused  by  the 
application  of  5  111M  NH3/NH4+  (and  presumably  higher 
concentrations  as  well) — and  NH44  or  H4  efflux  as  an 
explanation  for  the  paradoxical  pHj  increase  caused  by 
removal  of  high  levels  of  NH3/NH44 — we  have  no  delin- 
ltive  explanation  for  either  pH,  change.  Others  have 
proposed  that  members  of  the  Rh  family  function  as 
NH3/NH44  sensors  in  determining  the  choice  of  slug  vs. 
culmination  in  Dictyostelium  discoideum  (Kirsten  et  al 
2005,  2008;  Singleton  et  al.  2006).  We  propose  that 
oocytes  have  a  low-affinity  “sensor”  that  responds,  directly 
or  indirectly,  to  extracellular  NH3/NH4T  This  oocyte 
sensor  could  be  either  at  the  outer  surface  of  the  plasma 
membrane  or  somewhere  reasonably  close  to  the  inner 
surface,  but  ultimately  must  act  in  two  ways,  (a)  l  he  sensor 
triggers  the  production  of  cytosolic  H4 — perhaps  by  a 
metabolic  pathway.  The  influx  of  NH3  would  temper  the 
fall  in  pHj  and  lead  to  the  formation  of  some  NH4+,  which 
also  would  temper  the  depolarization  by  reducing  the 
inwardly  directed  NH44  diffusion  potential.  Indeed,  AmtB 
reduces  the  depolarization  (compare  AmtB  vs.  H30  bars  in 
Fig.  4d  and  h).  However,  ultimately,  the  overwhelming 
majority  of  the  incoming  NH3  must  be  trapped  as  NH4+  in 
a  presumably  acidic  intracellular  compartment  (b)  The 
oocyte  sensor  triggers  the  slow  activation  of  a  channel 
permeable  to  NH4  +  .  We  have  no  data  to  address  the  issue 
of  whether  this  activation  requires  the  attendant  fall  in  pH,, 
or  whether  the  putative  increase  in  NH44  permeability  and 
the  observed  fall  in  pHj  are  totally  independent.  Note  that 
NH44  need  not  enter  in  order  to  depolarize  the  cell.  We 
suggest  that  removal  of  extracellular  NH3/NH4  reverses 
this  production,  leading  to  consumption  of  H4  and  thus  the 
recovery  of  pHj.  Because  the  pHj  recovery  begins  so  soon 
after  NH3/NH44  removal,  we  suggest  that  the  NH3/NH44 
sensor  faces  the  extracellular  fluid. 

Vm  Changes 

As  noted  earlier,  the  slowly  developing  depolarization  that 
develops  in  the  presence  of  NH3/NH4+  could  reflect  per¬ 
meability  to  NH4+,  as  suggested  by  others  (Burckhardt  and 
Fromter  1992).  It  is  interesting  to  note  that,  at  5  both  mM 
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and  0.5  mM  NH^/NH4+,  expression  of  AmtB  substantially 
reduced  the  NHyNH4f-induced  depolarization  (Fig,  4d 
and  h)  without  significantly  affecting  pH,  (Fig,  4b  and  f). 
We  hypothesize  that  the  additional  influx  of  NH^  through 
AmtB  leads  to  the  formation  of  modest  amounts  of  NH4+ 
immediately  beneath  the  plasma  membrane,  reducing  the 
diffusion  gradient  for  NH44. 

Possible  Benefits  of  the  Oocyte’s  Unusual  Handling 
of  NHi/NH4+ 

An  intriguing  question  that  remains  is  why  the  oocyte 
should  handle  NH^  in  such  an  unusual  manner.  One  pos¬ 
sibility  is  that  the  oocyte’s  responses  to  extracellular 
NH-*/NH4+  are  an  adaptation  that  protects  the  cell — and 
perhaps,  more  importantly,  its  developmental  program — 
from  the  appearance  of  NH^  in  pond  water  that  contains 
decaying  organic  matter  and  thus  NHj/NH4+.  Levels  up  to 
at  least  0.5  mM  NH4/NH4+  cause  no  discernible  changes 
in  pH,,  and  even  much  higher  levels  lead,  at  most,  to 
modest,  slow,  and  fully  reversible  pH,  changes.  In  a  more 
typical  response,  an  exposure  to  NHi/NH44  might  mimic 
the  rise  in  pH,  caused  by  the  fertilization  of  a  Xenopus 
oocyte  (Webb  and  Nuccitelli  1981;  Nuccitelli  et  al  1981). 
The  Xenopus  oocyte  seems  to  be  particularly  adept  at 
avoiding  NH3/NH4 4 -induced  pH,  increases. 
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Abstract 

Aquaporins  arc  a  family  of  membrane  proteins  specialized  in  rapid  water  conduction  across  biological  membranes.  Whether  these 
channels  also  conduct  gas  molecules  and  the  physiological  significance  of  this  potential  function  have  not  been  well  understood.  Here 
we  report  140  ns  of  moleeular  dynamics  simulations  of  membrane-embedded  AQP1  and  of  a  pure  POPE  bilaycr  addressing  those  ques¬ 
tions.  The  permeability  of  AQPI  to  two  types  of  gas  molecules.  O2  and  CCT,  was  investigated  using  two  complementary  methods,  name¬ 
ly,  explicit  gas  dilTusion  simulation  and  implicit  ligand  sampling.  The  simulations  show  that  the  eentral  {tetra meric)  pore  of  \QP1  can  be 
readily  used  by  either  gas  molecule  to  permeate  the  channel  The  two  approaches  produced  similar  free  energy  profiles  associated  with 
gas  permeation  through  the  central  pore;  a  -0  4  to  —1.7  keal/mol  energy  well  in  the  middle,  and  a  3.6  4.6  keal/mol  energy  barrier  in  the 
periplasmic  vestibule.  The  barrier  appears  to  be  mainly  due  to  a  dense  duster  of  water  molecules  anchored  in  the  pcriplasmic  mouth  of 
the  central  pore  by  four  aspartate  residues.  Water  pores  show  a  very  low  permeability  to  O2,  but  may  contribute  to  the  overall  perme¬ 
ation  of  C02  due  to  its  more  hydrophilic  nature.  Although  the  central  pore  of  AQPI  is  found  to  be  gas  permeable,  the  pure  POPE  bilaycr 
provides  a  much  larger  cross-sectional  area,  thus  exhibiting  a  much  lower  free  energy  barrier  for  CO:  and  CL  permeation.  As  such,  gas 
conduction  through  AQPI  may  only  be  physiologically  relevant  cither  in  membranes  of  low  gas  permeability,  or  in  cells  where  a  major 
fraction  of  the  cellular  membrane  is  occupied  by  AQPs 
©  2006  Elsevier  Inc.  All  rights  reserved. 

Keywords:  Aquaporin;  AQPI,  Gas  permeability;  Oy  CO:;  Krce  energy  profile.  Gas  channels 


1.  Introduction 

Aquaporins  (AQPs)  arc  a  family  of  membrane  channels 
(Preston  et  al.,  1992,  Heymann  and  Engel,  1999;  de  Groot 
and  Grubmiiller,  2005;  Agre  et  al.,  1998;  Agre,  2006;  Borg- 
nia  et  al.,  1999;  Fujiyoshi  et  al.,  2002;  Hohmann  et  al., 
2001)  that,  through  modulating  the  water  permeability  of 
biological  membranes,  play  a  significant  role  in  water 
homeostasis  in  living  cells  (Agre  and  Kozono,  2003;  Deen 
and  van  Os,  1998;  Li  and  Verkman,  2001,  King  and  Agre, 
2004;  Agre  et  al.,  2004).  AQPs  organize  as  tetramers  in  the 


Corresponding  aulhor. 

E-mail  address:  emadfadife. uiuc.cdu  ( F .  Tajkhorshid). 

1 047-8477/S  -  see  front  maticr  ©  2006  hlscvicr  Inc.  All  rights  reserved 
doi:  10. 1016/j  jsb  2006.1 1.008 


membrane,  with  each  monomer  forming  a  functionally 
independent  water  pore  (Tajkhorshid  et  al.,  2002;  de  Groot 
and  Grubmiiller,  2001).  The  function  of  a  fifth  pore,  the 
'‘central  pore",  which  is  located  at  the  center  of  the  four 
monomers,  remains  unclear.  Due  to  their  structural  sim¬ 
plicity  and  the  approachable  timescale  of  their  basic  func¬ 
tion,  e.g.,  water  conduction,  AQPs  have  been  investigated 
extensively  with  molecular  dynamics  (MD)  simulations 
(Tornroth-Horseficld  ct  al.,  2006;  Tajkhorshid  et  al., 
2002,  2003,  2005a, b,  Jensen  et  al.,  2001,  2003;  Zhu  et  al., 
2001,  2004;  de  Groot  and  Grubmiiller,  2001,  2005;  Gray¬ 
son  et  al.,  2003).  These  studies  have  significantly  contribut¬ 
ed  to  our  understanding  of  the  molecular  basis  of  their 
function  and  selectivity.  A  particularly  intriguing  property 
of  AQPs  is  their  ability  to  block  protons  (Pohl  ct  al.,  2001; 
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Saparov  et  al.,  2005)  while  allowing  water  to  pass,  a  prob¬ 
lem  that  has  attracted  much  attention  from  theoretical 
studies  (Burykin  and  VVarshel,  2004;  Chakrabarti  et  al  , 
2004;  Tajkhorshid  et  al.,  2002;  de  Groot  et  al.,  2003;  Han 
et  al.,  2004;  de  Groot  and  Grubmuller,  2005;  Jensen 
et  al  ,  2003;  Burykin  and  Warshel,  2003). 

Water  permeation  is  by  far  the  best  characterized  physio¬ 
logical  role  of  AQPs,  however,  they  also  participate  in  other 
diverse  cellular  functions.  Permeabilities  of  other  small  sub¬ 
strates  in  AQPs  have  been  reported  for  glycerol  ( Borgnia  and 
Agre,  2001;  Grayson  et  al.,  2003;  Heller  et  al.,  1980),  urea 
(Borgnia  et  al  ,  1999),  nitrate  (Ikeda  et  al.,  2002),  arsenite 
(Liu  et  al.,  2002)  and  even  ions  (Yool  and  Weinstein, 
2002).  Furthermore,  some  AQPs  have  been  suspected  to  par¬ 
ticipate  in  gas  conduction  across  cellular  membranes  (Hanba 
et  al  ,2004;  Nakhouletal  ,  1998;  Herrera  et  al.,  2006;  Cooper 
and  Boron,  1998;  Endeward  et  al.,  2006;  Cooper  et  ah,  2002; 
Blank  and  Ehmke,  2003,  Terashima  and  Ono,  2002;  Uehlein 
et  ah,  2003,  Prasad  et  ah,  1998).  For  instance,  a  recent  study 
on  tobacco  aquaporin  NtAQPl  has  suggested  that  NtAQPl 
not  only  facilitates  C02  permeation  when  expressed  in  Xeno- 
pus  oocytes,  but  also  has  a  significant  function  in  photosyn¬ 
thesis  and  in  stomatal  opening  (Uehlein  et  ah,  2003). 
Another  study  on  Xenopus  oocytes  has  shown  that  the  C02 
permeability  is  increased  by  about  40%  when  AQP1  is 
expressed  (Nakhoul  et  ah,  1998;  Cooper  and  Boron,  1998). 
AQP1  reconstituted  in  proteoliposomes  has  been  found  to 
increase  both  water  and  C02  permeabilities  markedly  (Pra¬ 
sad  et  ah,  1998).  However,  as  lipid  membranes  are  often 
gas-permeable,  characterizing  the  gas  permeability  of  AQPs 
is  a  challenging  experimental  problem  (Verkman,  2002; 
Fang  et  ah,  2002;  Yang  et  ah,  2000).  Moreover,  the  pathway 
through  which  gas  molecules  might  permeate  AQPs  cannot 
be  identified  directly  through  presently  available  experimen¬ 
tal  methods.  Therefore,  it  is  of  great  interest  as  well  as  impor¬ 
tance  to  further  explore  the  gas  permeability  of  these 
membrane  channels  using  computational  methodologies. 
A  recent  simulation  study  has  employed  umbrella  sampling 
to  calculate  the  energetics  associated  with  C02  permeation 
through  AQP1  (Hub  and  de  Groot,  2006). 

To  shed  light  on  the  gas  permeability  of  AQP1  and  on 
its  physiological  relevance,  we  performed  a  total  of 
140  ns  MD  simulations  on  a  membrane-embedded  AQPI 


tetramer,  as  well  as  on  pure  POPE  bilayers.  Using  two  dif¬ 
ferent  methodologies,  namely  explicit  gas  diffusion  simula¬ 
tion  and  implicit  ligand  sampling  (Cohen  et  al.,  2006),  we 
have  analyzed  potential  pathways  for  permeation  of  02 
and  C02  through  AQPI,  and  calculated  the  associated 
energetics.  In  the  gas  dillusion  simulations,  100  C02  or 
02  molecules  were  initially  placed  in  water  on  both  sides 
of  the  membrane  and  allowed  to  move  freely.  The  implicit 
ligand  sampling,  on  the  other  hand,  predicts  the  PMF 
purely  based  on  simulations  of  membrane-embedded 
AQPI  in  the  absence  of  explicit  gas  molecules.  Although 
the  two  methods  work  through  disparate  mechanisms,  they 
have  yielded  similar  results.  The  PMFs  associated  with  gas 
permeation  revealed  a  —0  4  to  —  1 .7  keal/mol  energy  well  in 
the  middle  of  the  AQPI  central  pore,  and  an  energy  barrier 
of  3.6-L6  keal/mol  in  the  periplasmie  entrance.  Our  results 
suggest  that  the  hydrophobic  central  pore  of  AQPI  is 
indeed  permeable  to  both  C02  and  02.  Compared  with 
the  well-characterized  water  pores,  the  central  pore  of 
AQPs  is  poorly  understood,  and  it  remains  a  long-standing 
question  whether  it  has  any  physiological  importance.  Our 
simulations  suggest  that  the  AQPI  central  pore  provides  a 
potential  pathway  for  gas  permeation,  and  that  this  path¬ 
way  can  play  a  physiological  role  either  in  membranes  with 
low  intrinsic  gas  permeabilities,  or  when  a  major  fraction 
of  the  membrane  is  occupied  by  AQPs. 

2.  Methods 

2.1 .  Explicit  gas  diffusion  simulations 

Explicit  gas  dillusion  simulations  were  performed  on 
an  AQPI  tetramer  embedded  in  a  POPF’  bilayer,  with 
100  C02  or  02  molecules  initially  placed  in  the  water 
solution.  Control  simulations  were  performed  on  pure 
POPE  bilayers  without  AQPI.  Altogether,  four  systems, 
AQP1-C02,  AQPl-02,  POPE-C02  and  PQPE-02,  were 
investigated  using  two  types  of  simulations,  equilibrium 
simulations  and  induced  pressure  simulations,  as 
described  below  and  summarized  in  Table  1.  The  first 
two  systems  consisted  of  a  membrane-embedded  AQPI 
tetramer,  while  the  last  two  were  composed  of  pure 
POPE  bilayers 


Table  I 


The  lour  systems  investigated  using  explicit  gas  diffusion  simulations 


System 

/eq(  US) 

ACq 

E iw 

D  { 10  '  cmz/s) 

/p  (11s) 

POPE-C()2 

21 

25 

4.0  ±0.1 

1.80  ±0  01 

10 

203 

POPL-Oi 

15 

16 

4.7  ±0.2 

2  14  ±001 

10 

258 

AQPI -CO: 

30 

1 

11  I'd 

1  43  ±0  09 

10 

80(6.4) 

AQPI-O2 

26 

3 

11 /a 

2.39  ±0  06 

10 

83(1.6) 

n/a.  not  applicable. 

The  equilibration  time  (/eq)  and  the  mini  her  of  permeation  events  during  equilibrations  (  AW  are  listed.  Part il ion  coefficients  between  lipids  and  water  (A'lw) 
are  shown  for  POPL-CCN  and  POPH-02  Gas  diffusion  coefficients  {£))  in  water  are  listed  for  POPE-CCF  and  POPE-02;  gas  diffusion  coefficients  in  the 
central  pore  of  AQPI  (D)  arc  shown  for  AQPl-CCN  and  AQP1»02  The  simulation  time  for  induced  pressure  simulations  (rp)  and  ihc  number  of 
permeation  events  (Arp)  are  also  listed  Numbers  in  parenthesis  are  permeation  events  through  the  monomeric  water  pores  and  ihe  central  pore  of  AQPI, 
respectively  A.lvv  was  not  calculated  for  the  AQP  1/bilayer  systems  due  to  the  inhomogeneity  of  the  membrane  composition. 
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2. 1  1.  Modeling  o 

A  100  x  100  A2  POPE  bilaycr  patch  was  initially  built 
using  VMD’s  membrane  plugin  (Humphrey  et  al.,  1996). 
Two  20-A  layers  of  water  were  then  added  using  the  solvate 
plugin  of  VMD  (Humphrey  et  al.,  1996)  to  fully  hydrate 
the  system;  15  sodium  and  15  chloride  ions  were  then 
added  using  VMD's  autoionize  plugin  (Humphrey  et  al  , 
1996),  generating  a  50  mM  ionic  concentration.  After 
5000  steps  of  minimization,  the  lipid  tails  were  “melted" 
over  the  course  of  a  500  ps  of  constant  temperature 
(310  K)  simulation,  in  which  only  the  lipid  tails  were 
allowed  to  move  freely  while  the  rest  of  the  system  was 
fixed.  Another  500  ps  simulation  was  then  performed  under 
constant  temperature  (310  K)  and  constant  pressure 
(1  atm)  conditions  with  all  the  atoms  free  to  move.  Using 
the  last  100  ps  of  this  simulation,  we  calculated  the  free  vol¬ 
ume  inside  lipids  along  the  membrane  normal,  with  VMD's 
volmap  plugin  (Humphrey  et  al.,  1996).  The  equilibrated 
POPE  bilayer  was  then  used  to  build  the  POPE-CO2  and 
POPE-O2  systems,  where  100  C02  or  02  molecules  were 
initially  placed  in  the  bulk  region  among  water  molecules 
on  both  sides  of  the  membrane.  The  AQPl-COi  and 
AQPl-02  systems  were  built  in  a  similar  fashion,  with 
the  POPE-AQP1  structure  taken  from  a  450  ps  equilibra¬ 
tion  performed  by  Zhu  et  al.  (2001,  2004).  During  the 


a 

20  A 


wpe"Bla.: 

«  ** 


sm 


cytoplasm 


equilibration  of  POPE-AQP1  (Zhu  et  al.,  2001;  Zhu 
et  al.,  2004),  a  constraint  on  the  H-bond  between  the  guan- 
idinium  group  and  the  backbone  oxygen  of  Argl97  was 
used  to  prevent  the  flipping  of  the  side  chain  which  would 
block  the  channel  (Wang  et  al  ,  2005;  Jiang  et  al.,  2006). 
This  constraint  was  kept  in  our  explicit  gas  diffusion  simu¬ 
lations.  The  final  systems:  POPE-C02  (82,989  atoms), 
POPE-O,  (82,889  atoms),  AQPl-CCF  (81,065  atoms)  and 
AQPI-O2  (81,052  atoms)  share  a  similar  cross-sectional 
area,  with  the  membrane  normal  along  the  r-axis. 


2. 1  2.  Equilibrium  simulations 

All  the  four  systems  were  first  minimized  for  5000  steps. 
As  shown  in  Table  1,  POPE-C02  and  POPE-02  were  then 
equilibrated  at  constant  temperature  for  23  11s  and  15  ns, 
respectively;  AQP1-C02  and  AQPl-02  were  equilibrated 
under  constant  temperature  and  constant  pressure  condi¬ 
tions  for  30  ns  and  26  ns,  respectively.  Fig.  1  shows  snap¬ 
shots  of  the  AQP1-C02  and  POPE-CO?  systems  after  500 
ps  of  equilibration. 


2.1.3.  Induced  pressure  simulations 

After  the  first  10  ns  of  equilibrium  simulation  in  each 
system,  we  applied  a  force  of  20  pN  on  the  center  of  mass 
of  every  gas  molecule  along  the  — r  direction  in  a  10  11s  4  in- 
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Pig  1.  Side  (left  panels)  and  lop  (right  panels)  views  of  the  systems  AQP1-CCF  (a,b)  and  POPE-CCF  (e.d)  after  500  ps  of  equilibration  The  two  systems 
have  a  similar  cross-sectional  area,  with  the  membrane  normal  along  the  r-axis.  C02  molecules,  initially  placed  in  the  bulk  region,  are  drawn  in  vdW 
representation 
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dueed  pressure”  simulation  (Table  1).  Previous  simulations 
(Zhu  et  al.,  2002,  2004)  had  shown  that  such  an  external 
force  can  generate  an  effect  equivalent  to  a  concentration 
gradient  across  the  membrane,  favoring  the  difiusion  of 
gas  molecules  in  one  direction.  Using  the  force  applied  in 
our  simulations,  a  pressure  difference  of  ~20  MPa  was  gen¬ 
erated  across  the  membrane. 


Time  (ns) 


2.1  4.  Free  energy  profile  calculation 

The  free  energy  profile  associated  with  gas  permeation 
was  calculated  from  equilibrium  simulations  using  the  rela¬ 
tion  (Marrink  and  Berendsen,  1994): 

AC(.-)  = -firing  (1) 

L  bulk 

where  R  is  the  gas  constant,  Tis  the  temperature,  Cbuik  is  the 
concentration  of  gas  molecules  in  bulk  water,  and  C\z)  is  the 
concentration  of  gas  molecules  in  each  1  A  thick  layer  paral¬ 
lel  to  the  x)  plane  along  the  r-axis.  For  AQP1-C02  and 
AQPI-O2,  C(r)  is  the  gas  concentration  inside  the  central 
pore,  except  in  the  bulk  water  region  (r  ^  29  A  or 
-26  A),  where  gas  concentration  in  the  entire  layer 
was  calculated  and  then  normalized  according  to  the  area 
of  an  AQP1  tetramer  (50  nnr).  This  normalization  ensures 
that  our  calculated  PMF  corresponds  to  a  definite  density 
of  AQP1  ( 1  AQP1  per  50  nnr  of  bilayer).  Asa  result,  by  com¬ 
paring  the  PMFs  for  AQP1  and  POPE,  we  can  directly  deter¬ 
mine  whether  replacing  a  patch  of  lipids  by  AQP1  will 
increase  or  decrease  the  transmembrane  barrier  to  gas 
permeation. 

2.1  5.  Diffusion  coefficient  calculation 

The  gas  difiusion  coefficient  D  was  calculated  from  the 
slope  of  the  mean  square  displacement  against  time  using 
the  relationship 

<(/•(/  +  Al)  -  r(l)f)  =  2ilDAt  (2) 

where  r  is  the  coordinates  of  gas  molecules,  cl  is  the  dimen¬ 
sion,  and  At  is  an  arbitrary  time  interval.  As  gas  molecules 
may  travel  through  regions  with  different  diffusion  coeffi- 
cients  in  the  chosen  time  interval,  we  are  limited  to  small  val¬ 
ues  for  At  (here  At  =  5  ps),  so  that  D  best  represents  the  local 
gas  difiusion  coefficient  in  the  region  it  is  calculated  for 
(Marrink  and  Berendsen,  1994;  Shinoda  et  al.,  2004)  In 
determining  D  from  the  slope  of  the  mean  square  displace¬ 
ment  against  time,  the  first  2  ps  was  disregarded,  in  order 
to  ensure  that  the  calculated  diffusion  coefficient  is  not  con¬ 
taminated  b>  the  velocity  correlation  of  gas  molecules.  The 
time  origin  t  was  shifted  to  make  use  of  the  entire  trajectory 
and  improve  statistics.  A  resolution  of  2  A  was  used  to  com¬ 
pute  gas  diffusion  coefficients  in  lipid  bilayers. 

Both  the  free  energy  profile  and  diffusion  coefficients 
must  be  calculated  after  the  system  reaches  equilibrium. 
The  gas  density  of  the  lipid  bilayer  becomes  a  natural 
indicator  for  the  state  of  the  system,  since  we  put  all 
100  gas  molecules  in  bulk  water  at  the  beginning  of  sim¬ 
ulations,  the  system  is  approximately  at  equilibrium 


z  (A) 

Fig.  2.  Partition  of  gas  molecules  inlo  lhe  membrane  during  the  explicit 
gas  diffusion  simulations,  (a)  Number  of  gas  molecules  inside  the  POPE 
bilayer  over  lime  (POPE-COi:  black,  POPE-Oy.  green).  The  lipid  water 
interface  is  defined  as  the  average  position  of  phosphorus  aloms  of  lhe 
POPE  monolayers,  (b)  Lateral  (solid  lines)  and  traversal  (dashed  lines) 
diffusion  cocthcients  of  gas  molecules  (C02:  black.  02:  green)  in  POPE. 
The  percentage  of  free  volume  along  the  membrane  normal  in  the  POPE 
bilaycr  is  shown  in  blue.  Statistical  errors  of  the  diffusion  coefficients  are 
smaller  than  0.2  x  10  ■  cm  !/s  al  the  lipid  headgroup  region,  and  smaller 
than  0.1  x  10  '  cm‘/s  elsewhere,  (c)  Free  energy  profiles  calculated  using 
explicit  gas  diffusion  simulations  (exp)  and  implied  ligand  sampling  (imp). 
Slalislical  errors  of  lhe  PMFs  (exp)  are  smaller  lhan  0.25  keal/mol  al  the 
lipid  headgroup  region,  and  smaller  lhan  0. 18  keal/mol  elsewhere;  lhe 
statistical  error  for  PMF  (imp)  is  smaller  than  0  25  keal/mol  al  all  poinis. 
In  both  (b)  and  (e),  r  0  is  defined  as  lhe  mid-point  of  lhe  phosphorus 
planes  of  the  two  POPE  monolayers 

when  the  net  flux  of  gas  molecules  across  the  lipid-water 
interfaces  is  zero,  and  the  number  of  gas  molecules  inside 
the  lipid  bilayer  remains  constant.  Using  this  indicator, 
both  the  POPE-CO2  and  POPE-02  simulations  are  in 
equilibrium  during  their  last  5  ns  of  equilibration 
(Fig.  2a).  For  the  AQP1 -embedded  membranes,  we 
allowed  for  a  longer  equilibration  time  and  discarded 
the  first  10  ns  of  the  simulation,  We  then  investigated 
the  convergence  of  the  PMFs  calculated  using  the 
remaining  data.  The  PMFs  calculated  using  the  last  20, 
15  and  10  ns  of  the  simulations  were  rather  similar  for 
the  system  AQP1-C02  (data  not  shown).  Therefore,  we 
chose  to  use  the  last  20  ns  for  the  analysis,  in  order  to 
maximize  our  sampling.  As  for  AQPl-CK  we  resorted 
to  a  better  equilibrated  system  by  using  the  last  10  ns 
simulation,  except  for  the  periplasmic  part  of  the  central 
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pore  ( 16  A  ^  r  ^  29  A),  where  we  used  the  last  16  ns  due 
to  a  limited  sampling  in  this  region 


2.2.  Implicit  ligand  sampling 


Implicit  ligand  sampling  is  a  method  to  infer  the  3D 
PMF  of  a  weakly  interacting  ligand,  based  purely  on  an 
equilibrium  simulation  of  the  ligand-free  system.  This 
method  overcomes  the  problem  of  low  sampling  caused 
by  the  slow'  diffusion  of  explicit  ligands  exploring  a  protein 
and  has  proven  useful  in  finding  gas  migration  pathways  in 
proteins  (Cohen  et  al.,  2005,  2006).  Multiple  conformations 
of  an  imaginary  ligand  molecule  are  placed  at  every  point 
of  a  grid,  on  top  of  an  equilibrium  trajectory  of  a  simulated 
protein  system.  At  each  time  step,  and  for  each  ligand  posi¬ 
tion  and  orientation,  the  protein  ligand  interactions  are 
computed.  These  interaction  energies  are  then  used  to 
reconstruct  the  PMF  of  placing  the  ligand  at  any  point 
inside  the  protein  using  a  perturbative  approach,  as 
described  in  detail  in  Cohen  et  al  (2006). 

In  order  to  compute  the  PMFs,  equilibrium  simulations 
of  pure  POPE  and  membrane-embedded  AQP1  in  the 
absence  of  gas  molecules  were  performed  for  1  and  5  ns, 
respectively.  The  simulation  protocols  are  identical  with 
those  used  in  explicit  gas  diffusion  simulations,  except  that 
the  temperature  was  set  to  300  K,  and  the  H-bond  con¬ 
straint  on  Argl97  was  not  applied.  For  each  of  the  two  sys¬ 
tems,  a  3D  PMF  map  at  1  A  resolution  was  computed, 
using  the  volmap  feature  of  VMD  (Humphrey  et  al , 
1996)  with  a  chargeless  model  for  02  (Table  SI).  For  the 
lipid  bilayer  and  bulk  water  regions,  the  projected  PMFs 
along  the  r-axis.  Cfr),  were  then  calculated  from  the  3D 
PMF  map,  using  the  formula: 


F(=) 


-RT\nJ2 
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C-F(x.y;)/RT 


LXLy 


(3) 


where  F(x,)\z)  is  the  local  3D  PMF  computed  by  implicit 
ligand  sampling,  and  Lv  and  Lv  are  the  dimensions  of  the 
PMF  map  in  units  of  A.  In  the  case  of  AQP1,  the  calcula¬ 
tion  is  the  same  as  in  Eq  (3),  except  that  the  summations 
were  performed  over  restricted  cross-sectional  areas,  in  or¬ 
der  to  isolate  the  targeted  gas  channels.  The  resulting  PMF 
was  then  shifted  by  —  RT\w(LxLj A{)),  where  LxLy  is  the 
area  of  the  summation,  and  Aq  is  roughly  the  area  occupied 
by  an  AQP1  tetramer  (Aq  —  50  nm2).  This  is  done  so  that 
our  absolute  values  of  the  PMF  can  be  compared  with 
those  from  explicit  gas  diffusion  simulations.  To  compare 
our  PMFs  with  experiments  or  other  simulations,  one 
would  simply  shift  our  barriers  by  the  term  —  RT\ii(Aj 
//i),  where  A  \  is  the  cross-sectional  area  of  the  system  being 
used  for  comparison. 


2.3.  Error  analysis 

For  the  explicit  gas  diffusion  simulations,  the  statistical 
errors  in  the  calculation  of  a  variable  X  were  estimated 


by  dividing  the  trajectory  into  two  sub-trajectories  of  the 
same  length  X  was  then  calculated  using  each  sub-trajecto¬ 
ry.  We  used  the  quantity  std(3f)/\/2  as  a  measure  of  the 
statistical  error  on  A' in  our  calculation  using  the  entire  tra¬ 
jectory.  For  the  implicit  ligand  sampling,  the  estimated 
error  depends  on  both  the  measured  PMF  and  the  number 
of  independent  simulation  snapshots  used  to  compute  it, 
and  was  calculated  as  described  in  Cohen  et  al.  (2006). 

2.4.  MD  simulation  protocols 

For  all  simulations,  the  program  NAMD  2.6  (Phillips 
et  al.,  2005)  and  the  CHARMM27  parameter  set  (Schlenk- 
rich  et  al.,  1996;  MacKerell  et  al.,  1998)  were  used.  Parame¬ 
ters  for  C02  or  02  were  taken  from  the  literature  ( Harris  and 
Yung,  1995;  Schlenkrich  cl  al.,  1996;  MacKerell  et  al  ,  1998) 
and  are  listed  in  Supplemental  material  (Table  SI).  Assum¬ 
ing  periodic  boundary  conditions,  the  Particle  Mesh  Ewald 
(PME)  method  (Darden  et  ah,  1993)  with  a  grid  density  of 
at  least  1/A3  was  employed  for  computation  of  long-range 
electrostatic  forces.  All  simulations  were  performed  with 
time  steps  of  1,  2,  and  4  fs  for  bonded,  noil-bonded  and 
PME  calculation,  respectively.  A  Langevin  piston  (Feller 
et  ah,  1995)  was  employed  to  maintain  the  pressure  at 
1  atm  in  constant  pressure  simulations.  Langevin  dynamics 
was  used  to  keep  the  temperature  constant  at  310  K  (damp¬ 
ing  coefficient  5  ps  !),  except  in  the  implicit  ligand  sampling, 
where  the  temperature  was  set  to  300  K.  During  induced 
pressure  simulations,  gas  molecules  were  not  coupled  to 
the  temperature  bath,  and  counter  forces  were  applied  to 
Cy  atoms  of  AQP1  and/or  heavy  atoms  of  the  POPE  bilayer 
to  prevent  net  translocation  of  the  system  along  the  direction 
of  external  force.  The  total  counter  force  on  the  membrane 
was  equal  to  the  total  external  force  applied  on  gas  mole¬ 
cules.  In  the  AQP 1  -C02  and  AQP 1  -02  systems,  the  counter 
force  was  distributed  between  the  protein  and  the  lipids 
according  to  the  ratio  of  their  areas  in  the  membrane. 

3.  Results  and  discussion 

3. 1 .  Gas  di) fusion  in  POPE  hi  layer 

Free  energy  profiles  associated  with  gas  permeation 
through  POPE  are  shown  in  Fig.  2c.  In  general,  the  bilayer 
has  a  rather  “fiat"  free  energy  profile  for  both  C02  and  02, 
with  the  largest  energy  barrier  being  1 .4  kcal/mol  and 
1.9  kcal/mol,  respectively.  Such  a  profile  allows  gas  mole¬ 
cules  to  readily  diffuse  through  the  lipid  bilayer,  a  behavior 
observed  in  our  simulations  and  reflected  in  the  number  of 
permeation  events  during  the  equilibrium  simulations;  on 
the  average  one  permeation  event  per  nanosecond  was 
observed  for  both  C02  and  02  (Tabic  I).  As  shown  in 
Fig.  2c,  the  PMFs  calculated  using  explicit  gas  diffusion  sim¬ 
ulations  and  implicit  ligand  sampling  are  in  very  good  agree¬ 
ment.  An  energy  well  of  -  L0  kcal/mol  for  C02  and  -  1  3  to 
—  1.5  kcal/mol  for  02  is  located  at  the  bilayer  center,  while 
a  small  barrier  of  0. 1-0.4  kcal/mol  for  C02  and  0.1- 
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0.6  kcal/mol  lor  02  has  to  be  overcome  at  the  headgroup 
region  where  gas  molecules  cross  the  lipid-water  interface. 
Our  free  energy  profiles  calculated  from  explicit  gas  diffusion 
simulations  arc  mostly  symmetrical  about  the  center  of  the 
bilayer;  the  slight  asymmetry  at  the  lipid  headgroup  region, 
where  a  0.3  keal/mol  ( for  C02)  or  0.5  keal/mol  ( for  02)  ener¬ 
gy  difference  between  the  two  halves  of  the  PMFs  is 
observed,  is  most  likely  due  to  limited  sampling. 

It  is  noteworthy  that  the  headgroup  region  is  also  the 
region  in  which  gas  molecules  exhibit  the  least  mobility. 
As  shown  in  Fig.  2b,  inside  the  bilayer,  traversal  (along 
the  r-axis)  gas  diffusion  is  almost  always  faster  than  the 
corresponding  lateral  diffusion;  however,  both  the  traversal 
and  lateral  diffusions  are  significantly  slower  in  the  head- 
group  region.  Fig.  2b  also  shows  the  available  free  volume 
calculated  from  a  100  ps  equilibration  of  the  POPE  bilayer. 
It  can  be  expected  that  gas  diffusion  across  a  lipid  bilayer  is 
proportional  to  the  available  free  spaee  in  the  bilayer  In 
our  simulations,  the  diffusion  coefficients  correlate  elosely 
with  the  free  volume  inside  the  bilayer  in  all  regions. 

In  our  explicit  gas  diffusion  simulations,  100  gas  mole¬ 
cules  were  used  in  a  ~8  x  10"  A3  simulation  box,  generating 
an  overall  concentration  of  ~200  niM.  Although  this  con¬ 
centration  is  higher  than  physiological  concentrations,  giv¬ 
en  the  weak  interaction  potentials  and  the  large  average 
distance  between  the  gas  molecules  only  marginal  interac¬ 
tions  between  them  are  expected  As  such,  we  do  not  expeet 
the  results  of  the  explicit  gas  diffusion  simulations  to  be  sig¬ 
nificantly  a  fleeted  by  the  high  concentration  of  gas  mole¬ 
cules.  The  calculated  diffusion  coefficients  for  C02  and 
02  in  water  are  1 .80  x  10  5  enr/s  and  2.14x  10  5  enr/s, 
respectively  (Table  I),  which  are  smaller  than  the  experi¬ 
mental  values  (2.93  x  10  5  cnr/s  lor  C02  at  313  K  (Tam* 
imi  et  ah,  1994),  and  3.0  x  10“ 5  cm2/s  for  02  at  310  K 
(Subczynski  et  ah,  1992)). 

I  he  02  partition  coefficient  calculated  from  our  simula¬ 
tions  is  4.7,  elose  to  the  experimentally  determined  value 
(3.9)  for  02  in  liposomes  (Moller  et  ah,  2005).  However, 
our  calculated  C02  partition  coefficient  (4.0)  is  bigger  than 
the  experimental  value  for  C02  in  olive  oil  ( 1 .7)  (Simon  and 
Gutknecht,  1980),  which  is  very  likely  caused  by  the  statis¬ 
tical  error  in  our  simulation  and  the  force  field  parameters 
for  C02.  Therefore,  it  is  expected  that  the  depth  of  the  free 
energy  well  of  C02  inside  the  POPE  bilayer  (Fig.  2e)  is  also 
overestimated. 

3.2.  CO 2  and  02  permeability  across  AQPI 

For  C02  and  02  permeability  across  AQPI,  we  again 
resort  to  a  combined  approach,  simulating  explicit  diffusion 
of  C02  (30  ns)  and  02  (26  ns)  in  AQPI,  and  supporting  our 
results  with  an  independent  calculation  of  the  02  PMF  using 
implicit  ligand  sampling.  As  was  the  ease  for  the  lipid  bilay¬ 
ers,  PMFs  calculated  using  both  approaches  are  in  a  good 
agreement,  exeept  that  the  implicit  ligand  sampling  deteets 
additional  favorable  areas  not  explored  in  the  explicit  gas 
diffusion  simulations,  which  are  far  from  the  water  solution 
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Fig.  3.  (a)  Side  view  of  the  simulated  AQPl-02  system,  showing  areas  of 
high  02  occupancy  predicted  hy  implicit  ligand  sampling  (0.6  keal/mol 
energy  isosurfaec  shown  in  black  wireframe)  compared  with  02  locations 
revealed  by  explicit  gas  diffusion  simulations  (in  red,  frames  taken  every 
10  psol  the  26ms  simulation).  The  02  Irajectories  were. symmetrized  over  ihe 
AQPI  let ramer  such  that  each  monomer  shows  the  locations  explored  hy  02 
in  all  four  monomers,  (b)  Top  view  of  ihe  AQPI  tetramcr.  along  with  i he 
regions  used  to  compute  projected  PM  Fs  for  each  of  AQP  l  \s  various  pores: 
(A)  ihe  central  pore.  { B)  ihe  water  pores.  (C)  the  side  pores  TheO  kcal/mol 
energy  isosurface  of  02  PMF  in  these  pores  is  shown  in  yellow. 

and  the  bilayer.  This  agreement  is  highlighted  in  Fig.  3a, 
where  favorable  locations  for  02  revealed  by  the  two 
approaches  are  shown.  In  our  assessment  of  the  gas  perme¬ 
ability  of  AQPI,  we  have  divided  the  AQPI  tet ramer  into 
various  regions  of  interest  which  we  describe  separately: 
the  central  pore,  the  monomerie  water  pores,  a  side  pore  that 
we  propose  as  an  additional  possible  gas  ehannel,  and  the 
protein-lipid  interface.  In  all  of  our  analyses,  we  define  these 
regions  as  shown  in  Fig.  3b. 

3.2.1  The  central  pore 

Two  quadruplets  of  hydrophobie  residues,  Val52  and 
Leu  172,  form  the  outermost  doorways  of  the  AQPI  eentral 
pore.  Water  molecules  can  hardly  pass  through  these 
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residues  and  were  kept  away  from  the  inside  of  the  pore  dur¬ 
ing  our  simulations,  a  behavior  also  revealed  by  other  MD 
studies  (Zhu  et  al.,  2001 ;  Yu  et  al.,  2006).  For  gas  molecules, 
however,  the  purely  hydrophobic,  empty  central  pore  pro¬ 
vides  an  ideal  reservoir.  At  the  end  of  the  explicit  gas  diffu¬ 
sion  simulations,  8  C02  and  5  02  molecules,  respectively, 
were  found  inside  the  central  pore,  as  seen  in  Fig.  4a  and  b. 

As  shown  in  Fig.  5a,  PMFs  calculated  from  both  the 
explicit  gas  diffusion  simulations  and  the  implicit  ligand 
sampling  suggest  that  the  central  pore  is  the  most  favorable 
location  for  both  C02  and  02.  However,  a  more  favorable 
central  pore  for  02  with  a  —  1.7kcal/mol  energy  well  was 
predicted  by  implicit  ligand  sampling,  compared  to  a 
-0.4  kcal/mol  energy  well  revealed  by  the  explicit  gas  diffu¬ 
sion  simulations.  This  difference  is  probably  caused  by  the 
limited  timescale  of  our  explicit  gas  diffusion  simulations: 
while  the  gas  distribution  between  the  water  and  lipid  phas¬ 
es  is  well  equilibrated,  the  gas  molecules  need  longer  time 
to  reach  an  equilibrated  population  inside  the  protein.  This 
is  especially  true  given  that  the  central  pore  is  capped  by 
areas  of  relatively  high  energy  barriers,  and  that  the  flow 
of  gas  molecules  across  the  central  pore  was  low  compared 
to  the  gas  population  inside  the  central  pore  during  the 
course  of  the  explicit  gas  diffusion  simulations. 

According  to  our  PMFs,  the  major  barrier  of  the  central 
pore  to  gas  permeation  is  about  3.6-4. 6  kcal/mol,  located 
at  the  periplasmic  side  of  the  central  pore 


Fig  4  (;i  and  b)  Accumulation  of  gas  molecules  inside  the  cenlral  pore  of 
\QP1  \t  ihe  end  of  the  equilibrium  Simula  lions,  eight  CO:  (a)  and  five 
02  (b)  were  found  in  the  central  pore,  respectively.  Only  Ihree  AQP1 
monomers  are  shown,  with  the  from  monomer  being  removed  for  clariiy. 
Residues  Val52  and  Leu  172.  which  form  the  outermost  doorways  of  lhe 
hydrophobic  central  pore,  are  highlighied  in  licorice  represeniaiion.  (c) 
Snapshot  showing  C02  in  one  of  lhe  monomeric  water  pores  during 
equilibrium  simulations.  Water  and  C02  molecules  are  shown  in  vdW 
representation.  Argl97  at  the  SF  (see  text)  and  Asnl94  and  Asn78  from 
the  NPA  motifs  are  shown  in  Licorice  representation. 


Fig.  5  (a)  PMFs  associated  with  gas  permeation  through  lhe  AQP1 
central  pore,  calcuiaied  from  the  explicit  gas'  diffusion  simulations  (exp), 
and  implicit  ligand  sampling  (imp).  The  doited  lines  correspond  lo  lhe 
PMFs  of  gas  molecules  in  water  of  the  same  area  as  the  central  pore  (b) 
CL  PMF  calculated  by  implicit  ligand  sampling  for  regions  shown  in 
Fig  3.  along  with  the  PMF  in  lhe  POPE  bilayer.  The  PMFs  were 
computed  relative  to  a  value  of  0  kcal/mol  for  CL  in  vacuum.  In  both  (a) 
and  (b)  r  =  0  is  defined  as  the  center  of  Asn78  and  Asnl94  of  the  NPA 
motifs.  All  PMFs  were  calcuiaied  assuming  a  density  of  1  AQP1  per 
50  nm"  of  bilayer,  as  described  in  Section  2  Statistical  errors  for  CCL 
PMF  (exp)  and  02  PMF  (exp)  are  less  lhan  0.4  kcal/mol  except  at  the 
energy  barriers,  where  due  lo  limited  sampling,  the  siaiisiieal  errors  cannot 
be  calculated  using  ihe  same  formula  mentioned  in  Section  2.  The  upper 
bound  error  of  the  PMF  (imp)  is  +0.25  kcal/mol.  lhe  lower  bound  errors 
are  0.25,  -0.6,  1.7  kcal/mol,  respectively,  for  PMF  values  below  4,  6, 

and  8  kcal/mol  (Cohen  ct  al.,  2006) 


( 13  A  ^  r  ^  19  A).  This  barrier,  consistently  found  by  both 
our  approaches,  as  well  as  by  another  study  (Hub  and  dc 
Groot,  2006),  surprisingly,  does  not  correspond  to  a  region 
that  is  sterically  blocked  directly  by  the  protein.  As  shown 
in  Fig  6a,  this  barrier  is  located  above  the  region  of  max¬ 
imum  protein  contraction  where  the  four  hydrophobic  res¬ 
idues  Val52  reside:  rather,  it  corresponds  to  a  region  that  is 
populated  solely  by  water  We  have  created  a  volumetric 
map  of  the  local  occupancy  of  water,  as  shown  in 
Fig.  6b  It  is  clear  that  the  barrier  corresponds  to  a  dense 
layer  of  water  molecules  surrounded  by  four  aspartate 
(Asp50)  residues  (Fig.  6c).  With  a  higher  density  than  the 
bulk  water,  this  water  layer  reduces  the  chance  of  gas  mol¬ 
ecules  to  access  the  central  pore.  If  these  aspartate  residues 
will  be  mutated  to  neutral  residues,  e.g  ,  alanines  or  aspar¬ 
agines.  the  strong  electrostatic  effects  of  the  quadruplets 
may  be  eliminated  and  a  less  dense  water  structure  could 
be  expected,  which  might  result  in  a  better  gas-conductivc 
central  pore  of  AQP1  Simulations  of  these  mutants  arc 
currently  underway. 
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Fig.  6.  (a)  Close-up  of  lhc  energy  barrier  (uidicaled  by  l he  black  arrow)  al  lhe  periplasmic  moulh  of  lhe  AQP1  central  pore  The  0  kcal/mol  (yellow)  and 
3  kcal/mol  (black)  energy  isosurfaces  calculated  from  implicit  ligand  sampling  are  shown,  along  with  02  positions  obtained  from  explicit  gas  diffusion 
simulations  (red),  (b)  VValcr  occupancy  per  A\  showing  in  blue  (he  volumes  occupied  by  water  during  at  least  75%  of  lhe  lime  during  the  26-ns  AQPl-CT 
simulation.  A  layer  of  water,  denser  than  lhe  bulk,  is  found  around  the  Asp50  residues,  corresponding  to  lhe  location  of  the  energy  barrier  shown  in  (a). 
Water  occupancy  was  calculated  using  lhe  volmap  plugin  of  VMD  (Humphrey  el  al  1996).  The  protile  of  lhc  central  pore  is  shown  in  black  lines,  (e) 
Snapshot  of  water  molecules  around  Asp50  during  equilibrium  simulations  of  AQPI-CT.  The  water  molecules  are  shown  in  both  surface  and  CPK 
representations. 


3.2.2.  The  w  etter  pores 

Compared  with  the  central  pore,  the  four  monomeric 
water  pores  in  AQP1  were  found  to  be  less  permeable  to 
gas  molecules.  In  our  gas  diffusion  simulations,  although 
both  C02  and  02  could  enter  the  two  vestibules  of  water 
pores  and  reach  the  NPA  motifs  (see  Fig.  4c),  no  gas  mol¬ 
ecule  passed  through  the  selectivity  filter  (SF),  />.,  the  nar¬ 
rowest  part  of  an  AQP  channel.  The  single  file  of  water 
molecules  at  the  SF  was  almost  always  kept  intact  during 
the  simulations.  This  may  be  explained  by  the  steric  barrier 
imposed  by  the  protein,  as  well  as  the  strong  tendency  of 
water  molecules  to  keep  their  hydrogen-bonded  structure; 
as  the  SF  has  a  diameter  of  only  —2  A,  gas  molecules  could 
hardly  pass  this  region  without  perturbing  the  water  file, 
which  is  held  strongly  by  hydrogen  bonds  between  water 
molecules  and  with  the  protein  Using  the  implicit  ligand 
sampling  approach,  we  measured  the  PMF  for  02  perme¬ 
ation  across  the  four  water  pores  counted  as  one.  As  shown 
in  Fig.  5b,  the  highest  energy  barrier  is  located  at  the  SF  and 
is  found  to  be  5.7  kcal/mol  with  respect  to  the  water  solution 
This  is  comparable  with  a  recent  study  which  measures  a 
—5.5  kcal/mol  energy  barrier  for  C02  permeation  through 
the  water  pore  of  AQP1  (Hub  and  de  Groot,  2006). 

The  central  pore  also  appears  to  be  more  gas-perme¬ 
able  than  the  water  pores  in  our  induced  pressure  simula¬ 
tions.  As  shown  in  Fig.  SI,  the  single  file  of  water  was 
largely  disturbed  when  C02  penetrated  the  water  pores 
during  induced  pressure  simulations.  In  contrast,  no  obvi¬ 
ous  disturbance  to  the  protein  was  observed  when  C02 
permeated  the  central  pore.  During  the  induced  pressure 
simulations,  four  C02  and  six  02  permeated  the  central 
pore,  while  six  C02  and  one  02  passed  through  the  water 


pores  (Table  1).  Given  that  there  are  four  water  pores,  the 
ratio  of  gas  permeation  between  the  central  pore  and  a 
water  pore  is  —3:1  for  C02  and  24  1  for  02.  The  differ¬ 
ence  between  C02  and  02  suggests  that  the  water  pores 
of  AQP1  are  less  favorable  to  the  more  hydrophobic  02 
molecules.  Indeed,  in  our  equilibrium  simulations,  the 
probability  of  finding  02  in  a  water  pore  of  AQP1  is 
much  lower  than  that  of  C02  (data  not  shown).  There¬ 
fore,  while  C02  might  take  advantage  of  both  the  central 
pore  and  water  pores  to  permeate  AQP1,  02  is  primarily 
conducted  through  the  central  pore. 

3.2.3.  The  side  pores 

An  interesting  and  somewhat  unexpected  result  of  our 
implicit  ligand  sampling  calculations  is  the  identification 
of  a  favorable  gas  channel  spanning  AQP1,  located 
between  each  pair  of  AQP1  monomers  near  the  protein 
lipid  interface.  These  “side  pores”,  which  to  our  knowledge 
have  not  been  considered  before,  have  periplasmic  and 
cytoplasmic  exits  which  are  completely  enclosed  inside 
the  AQP1  tetramer.  The  shape  of  the  side  pores  can  be  seen 
on  either  sides  of  Fig.  3a  The  02  PMF  for  the  four  side 
pores  (Fig.  5c)  suggests  that  it  offers  low  entry  and  exit  bar¬ 
riers  to  02  in  AQP1  along  with  a  flat  energy  profile  inside 
the  pore,  making  it  a  strong  candidate  for  a  favorable  gas 
pathway  across  AQP1 .  Furthermore,  because  the  middle  of 
these  four  pores  are  directly  connected  to  the  center  of  the 
lipid  bilayer,  the  side  pores  can  increase  the  rates  of  02 
entry  into  and  exit  from  the  bilayer,  further  facilitating 
permeation  across  the  membrane.  Although  no  complete 
permeation  events  through  the  side  pores  were 
recorded  in  our  equilibrium  simulations,  we  observed 
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several  half-permeation  events  where  gas  molecules  entered 
or  exited  the  bilayer  center  via  these  side  pores. 

3.2.4.  Protein  lipid  interface 

According  to  both  explicit  gas  diffusion  simulations  and 
implicit  ligand  sampling,  the  interfacial  region  between  the 
AQPf  tetramer  and  the  lipid  bilayer  has  an  increased  affinity 
for  02  as  compared  to  the  rest  of  the  bilayer.  Because  of 
AQPl’s  uneven  hydrophobic  surface  inside  the  bilayer,  the 
packing  with  the  lipids  is  not  optimal,  leaving  plenty  of  room 
for  small  hydrophobic  molecules.  However,  our  3D  PMF 
map  for  02  indicates  that  the  gas  density  in  this  region  is 
localized  in  many  disconnected  pockets  around  the  protein, 
and  there  is  no  single  continuous  pathway  connecting  the 
sides  of  the  membrane.  For  this  reason,  while  the  protein - 
lipid  interface  may  possibly  store  02,  it  does  not  necessarily 
conduct  it  well.  As  the  lipids  do  not  migrate  or  diffuse  at  the 
timescales  probed  by  either  of  our  analyses,  the  PMFs  we 
observed  may  not  be  accurate  and  the  possibility  of  gas  per¬ 
meation  at  the  protein-lipid  interface  cannot  be  completely 
discounted.  However,  we  consider  high-throughput  02 
migration  in  this  region  to  be  unlikely  on  the  basis  of  our 
computations. 

3.2.5.  AQPl  compared  with  POPE 

Our  free  energy  profiles  revealed  a  higher  energy  barrier 
for  gas  permeation  across  AQPl  than  a  patch  of  POPE 
bilayer  with  the  similar  area.  Besides  the  PMFs,  one  can 
compare  the  gas  permeabilities  of  AQPl  and  POPE  in  a 
rather  straightforward  way.  As  shown  in  Table  I,  on  aver¬ 
age  one  permeation  per  nanosecond  is  observed  for  POPE- 
C02  and  P0PE-02  during  equilibrium  simulations,  while 
few  permeation  events  were  recorded  for  AQPI-C02  and 
AQPl-02.  The  number  of  gas  permeation  events  naturally 
reflects  the  gas  permeability  of  the  system.  However,  better 
statistics  are  required  to  meaningfully  compare  the  AQP1- 
embedded  membrane  with  the  pure  POPE  bilayer.  This 
was  achieved  by  induced  pressure  simulations,  in  which 
gas  molecules  were  forced  to  cross  the  membrane  by  an 
external  force,  and  the  number  of  gas  permeation  events 
was  increased  by  approximately  an  order  of  magnitude, 
allowing  us  to  compare  the  gas-permeability  of  the  tw'o 
types  of  membranes  directly.  As  shown  in  Table  1,  203 
C02  and  258  02  permeation  events  were  observed  for  pure 
POPE  bilayers,  while  for  the  AQPl -embedded  membranes, 
the  number  is  80  and  83,  respectively.  Since  all  the  four  sys¬ 
tems  share  a  similar  x v  area,  the  insertion  of  AQPl  into  the 
POPE  bilayer  has  clearly  reduced  the  total  gas  permeability 
of  the  membrane. 

4.  Conclusions 

We  have  studied  the  gas  permeability  of  AQPl  by  per¬ 
forming  MD  simulations  on  a  membrane-embedded  AQPl 
tetramer  as  well  as  on  pure  lipid  bilayers.  Two  different 
approaches  were  used  to  calculate  the  free  energy  profiles 
associated  with  gas  permeation  through  AQPl  and  yielded 


consistent  results.  In  general,  AQPl  appears  to  be  less  gas- 
permeable  than  a  patch  of  POPE  bilayer  with  a  similar  area. 
However,  the  central  pore  of  AQPl  is  found  to  be  permeable 
to  both  C02  and  02  with  a  favorable  energy  well  of  -0.4  to 
-1.7  kcal/mol  at  the  center  and  a  3. 6^4. 6  kcal/mol  barrier 
on  the  periplasmic  side  It  is  noteworthy  that  this  barrier, 
revealed  by  both  of  our  approaches,  is  not  directly  caused 
by  the  protein,  but  rather  by  a  dense  layer  of  water  molecules 
‘‘anchored"  at  the  periplasmic  mouth  of  the  central  pore  by 
four  aspartate  residues. 

Both  the  gas  diffusion  simulations  and  the  implicit 
ligand  sampling  method  suggest  that  the  monomeric  water 
pore  of  AQPl  is  less  gas-permeable  than  the  central  pore, 
likely  due  to  the  strong  hydrogen  bonds  between  the  pro¬ 
tein  and  water  molecules  inside  the  channel.  A  somewhat 
unexpected  discovery  is  the  identification  of  “side  pores" 
located  in  between  AQPl  monomers,  which  may  facilitate 
gas  conduction  across  the  membrane.  However,  due  to  the 
fluctuation  of  the  protein,  these  side  channels  might  not  be 
permanently  open.  It  is  therefore  of  interest  and  impor¬ 
tance  to  further  investigate  the  role  of  these  side  pores  in 
gas  conduction  mediated  by  AQPs. 

Our  simulations  suggest  that  aside  from  its  conventional 
role  as  a  water  channel,  AQPl  might  provide  a  pathway  for 
small,  neutral  gas  molecules  across  the  cellular  membrane. 
However,  compared  with  a  patch  of  pure  POPE  bilayer 
with  a  similar  area,  AQPl  is  less  permeable  to  C02  and 
02.  Therefore,  gas  conduction  through  AQPl  may  be  of 
physiological  importance  only  in  membranes  with  low 
intrinsic  gas  permeability  or  where  a  major  fraction  of 
the  surface  area  of  the  membrane  is  occupied  by  AQPs. 
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The  water  channel  aquaporin  1  (AQP1)  and  certain  Rh-family 
members  are  permeable  to  CO2  and  NH3.  Here,  we  use  changes  in 
surface  pH  (pHs)  to  assess  relative  CO2  vs.  NH3  permeability  of 
Xenopus  oocytes  expressing  members  of  the  AQP  or  Rh  family. 
Exposed  to  CO2  or  NH3,  AQP1  oocytes  exhibit  a  greater  maximal 
magnitude  of  pHs  change  (ApHs)  compared  with  day-matched 
controls  injected  with  H20  or  with  RNA  encoding  SGLT1,  NKCC2,  or 
PepTI.  With  CO2,  AQP1  oocytes  also  have  faster  time  constants  for 
pHs  relaxation  (tpHs).  Thus,  AQP1,  but  not  the  other  proteins, 
conduct  C02  and  NH3.  Oocytes  expressing  rat  AQP4,  rat  AQP5, 
human  RhAG,  or  the  bacterial  Rh  homolog  AmtB  also  exhibit 
greater  ApHs(C02)  and  faster  tpHs  compared  with  controls  Oocytes 
expressing  AmtB  and  RhAG,  but  not  AQP4  or  AQP5,  exhibit  greater 
ApHs(NH3)  values.  Only  AQPs  exhibited  significant  osmotic  water 
permeability  (Pf).  We  computed  channel-dependent  (*)  ApHs  or  Pf 
by  subtracting  values  for  H20  oocytes  from  those  of  channel¬ 
expressing  oocytes.  For  the  ratio  ApHs(C02)*/Pf\  the  sequence  was 
AQP5  >  AQP1  3  AQP4.  For  ApHs(C02)*/ApHs(NH3)*,  the  sequence 
was  AQP4  AQP5  >  AQP1  >  AmtB  >  RhAG.  Thus,  each  channel 
exhibits  a  characteristic  ratio  for  indices  of  CO2  vs.  NH3  permeabil¬ 
ity,  demonstrating  that,  like  ion  channels,  gas  channels  can  exhibit 
selectivity. 

gas  channel  |  oocyte  |  permeability  |  signal  peptide  | 
surface  pH  measurement 

Gas  transport  through  membranes  is  of  fundamental  impor¬ 
tance  for  nutritive  transport,  photosynthesis,  oxidative  me¬ 
tabolism,  and  signaling.  For  most  of  the  past  century,  we 
assumed  that  gas  molecules  cross  biological  membranes  merely 
by  diffusing  through  the  lipid  phase.  This  dogma  was  challenged 
by  2  observations:  (t)  Apical  membranes  of  gastric-gland  cells 
have  no  demonstrable  permeability  to  CO:  or  NH3  (1).  (ii) 
Heterologous  expression  of  the  water  channel  aquaporin  I 
(AQPI)  increases  the  CO:  permeability  of  Xenopus  oocytes  (2). 
Cooper  and  Boron  (3)  and  Prasad  et  al.  (4)  confirmed  and 
extended  this  observation.  Uehlcin  (5)  showed  that  an  AQP 
plays  a  physiological  role  by  enhancing  CO:  uptake  by  plants. 
Endcward  et  al  (6)  demonstrated  that  AQP1  accounts  for  ^60% 
of  the  CO2  permeahility  of  human  red  blood  cells  (RBCs). 
Molecular  dynamics  simulations  suggest  that  CO2  can  pass 
through  the  4  aquaporcs  of  an  AQPi  tetramer  (7)  and  especially 
through  the  central  pore  between  the  4  monomers  (7).  Addi¬ 
tional  data  indicate  that  AQPI  is  permeable  to  nitric  oxide  (8), 
and  that — when  expressed  in  Xenopus  oocytes  (9,  10)  or  when 
reconstituted  into  planar  lipid  bilaycrs  (11) — AQPI  AQP3, 
AQPS,  AQP9,  and  the  plant  aquaporin  TFP2;  1  arc  all  permeable 
to  NH.v 

The  AmtB/MEP/Rh  proteins  represent  a  second  family  of  gas 
channels  (12-15).  Early  work  showed  that  AmtB  and  MEP 
transport  NH3  or  NH4,  thereby  playing  a  nutritive  role  in 
archaea,  bacteria,  and  fungi  ( 16, 1 7).  The  crystal  structures  of  the 
bacterial  AmtB  (18-20)  and  Rh50  (21)  and  the  fungal  Amt-1 
(22)  are  consistent  with  the  idea  that  NH^  passes  through  a  pore 
in  each  monomer  of  the  homotrimcr.  Indeed,  reconstituted 
AmtB  conducts  NH3  (14),  and  RhAG  is  necessary  for  NH3 


transport  in  mammalian  RBCs  (23).  Soupene  et  al.  found  that 
RhI  deficiency  impairs  the  grow  th  of  the  green  alga  C.  reinhardtii 
(24)  and  suggested  that  Rhl  plays  a  role  in  CO:  transport.  In 
RBCs,  RhAG  accounts  for  **50%  of  CO:  transport  (25). 

In  2006,  we  introduced  an  approach  (6)  to  assess  CO2  trans¬ 
port  by  pushing  a  blunt  microclectrodc  against  the  surface  of  an 
oocyte,  while  monitoring  surface  pH  (pHs)-  Introducing  extra¬ 
cellular  CO2  causes  a  transient  pHs  increase,  the  maximum 
magnitude  of  which  (ApHs)  is  an  index  of  maximal  CO2  influx. 
Earlier,  De  Hcmptinnc  and  Hugucnin  (26)  had  observed  such  a 
CCE-induccd  transient  while  monitoring  extracellular  pH  (pH()) 
of  rat  soleus  muscle.  Moreover,  Chcslcr  (27)  had  found  that 
exposing  lamprey  neurons  to  NH3  causes  a  transient  decreases  in 
pH0.  Here,  we  exploit  CO:-  and  NIH-induccd  pHs  transients  to 
study  the  CO2  vs.  NH3  permeability  of  4  channels  abundantly 
expressed  in  cells  that  mediate  high  rates  of  gas  transport:  human 
AQPI  (RBCs’  ref.  28),  the  M23  variant  of  rat  AQP4  (astrocytic 
endfeet  at  the  blood-brain  barrier,  ref.  29),  rat  AQP5  (alveolar 
type  I  pneumocytcs;  ref.  30),  and  human  RhAG  (RBCs,  ref.  31). 
We  also  studied  bacterial  AmtB.  Our  results  show  that  all  5 
channels  arc  permeable  to  CO2,  and  all  but  AQP4  and  AQPS  are 
permeable  to  NHv  A  relative  index  of  CO^NFH  permeability 
varied  widely:  AQP4  =  AQPS  >  >  AQPI  >  AmtB  >  RhAG. 
Thus,  as  is  true  for  ion  channels,  gas  channels  exhibit  substantial 
solute  selectivity,  which  could  play  an  important  physiological 
role  in  controlling  gas  fluxes. 

Results 

pHs  Transients  Caused  by  Applying  C02  vs.  NH3.  Fig.  \A  illustrates 
schematically  how-  the  influx  of  CO:  leads  to  a  fall  in  [CO:]  near 
the  extracellular  surface  of  the  membrane  ([CO:]s),  which  in 
turn  leads  to  a  rise  in  pHs-  Fig.  I B  shows  how  the  influx  of  NH3 
leads  to  a  fall  in  pH$.  As  described  in  ref.  6,  exposing  an 
AQPl-cxprcssing  oocyte  to  a  solution  containing  5%  CO:/33 
mM  HC03  at  a  constant  pH  of  7.50  causes  a  transient  rise  in 
pHs,  followed  by  an  exponential  decay  (Fig.  1C  Left ,  green 
record).  After  the  washout  of  CO:  (sec  SI  Text  and  Fig.  S3), 
exposing  the  same  oocyte  to  0.5  mM  NH3/NIII  causes  a 
transient  fall  in  pHs  (Fig.  1C  Right,  green  record),  as  noted 
elsewhere  (32).  Additional  data  are  consistent  with  the  hypoth¬ 
esis  that  Xenopus  oocytes  handle  NTH  in  an  unusual  way, 
sequestering  most  incoming  NH^  in  an  intracellular  compart¬ 
ment  as  NH4  (32). 

The  maximal  pHs  transients  are  much  smaller  in  day-matched 
oocytes  injected  with  H:0  (orange)  or  cRNA  encoding  the 
Na/glucosc  cotransporter  SGLTI  (black),  and  arc  totally  lacking 
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Fig.  1.  Basis  of  surface  pH  changes.  (4)  CO2  influx.  At  the  outer  surface  of  the  membrane,  the  CO2  influx  creates  a  CO2  deficit  The  reaction  HC03  +  H+  — *C02 
+  H20,  in  part,  replenishes  the  C02,  raising  pHs.  (8)  NH3  influx  The  reaction  NH4  — *•  NH3  +  H\  in  part,  replenishes  NH3,  lowering  pHs  (C -E)  Representative  pHs 
transients  from  oocytes  injected  with  H2O  or  expressing  AQP1  (record  repeated  in  the  3  images),  SGLT1,  NKCC2,  or  PepTI  All  data  in  C-f  were  obtained  on  the 
same  day,  from  the  same  batch  of  oocytes,  exposed  first  to  C02/HC03  and  then  (after  C02  removal)  to  NH^NHT  Also  shown  in  C  are  records  with  no  oocyte 
present.  Before  and  after  solution  changes,  we  retracted  the  pH  electrode  to  the  bulk  extracellular  solution  (pH  7.50)  for  calibration.  ( F  and  G)  Summary  of 
extreme  excursions  of  pHs  (ApHs)  for  C02  and  NH3  data.  (H  and  /)  Representative  pHs  transients  from  day-matched  H20  or  AQP1  oocytes  exposed  to  C02/HC03 
or  30  mM  butyrate.  ( J  and  K)  Summary  of  ApHs  for  experiments  like  those  in  H  or  /  Values  are  means  i  SE,  with  numbers  of  oocytes  in  parentheses.  For  F  and 
G,  statistical  comparison  between  H2O-  injected  controls  and  other  oocytes  (separately  for  C02  and  NH3  data)  were  made  using  a  1  -way  ANOVA  for  5  groups, 
followed  by  Dunnett's  multiple  comparison.  ApH5  values  for  H20,  NKCC2,  and  PepTI  do  not  differ  from  one  other.  For  J  and  K,  statistical  comparisons  were  made 
using  unpaired  2-tailed  t  tests 


m  the  absence  of  an  oocyte  (gray).  Moreover,  oocytes  expressing 
the  Na/K/Cl  cotransportcr  NKCC2  or  the  H/oligopeptide  co¬ 
transporter  PepTI  have  ApH$  values  similar  to  those  of  H2O- 
injected  oocytes  (Fig.  1  D  and  E). 

Fig.  1  F  and  G  summarize  ApHs  data  for  a  larger  number  of 
experiments  like  those  in  Fig.  1  C-E  and  show  that  ApHs  for 
AQP1  is  significantly  greater  than  for  all  other  oocyte  groups.  As 
expected,  the  time  constant  (rpns)  for  the  decay  of  pH$  from  its 
peak,  an  index  of  the  time  required  for  CO2  to  equilibrate  across 
the  membrane,  has  a  pattern  that  is  the  inverse  of  that  for  ApHs 
(Fig.  SI  I  and  /?).  Because  of  the  oocyte’s  unusual  NFF  handling, 
the  pFIs  relaxation  during  NFF  exposures  is  prolonged,  preclud¬ 
ing  the  calculation  of  a  tphs  for  NH3.  In  the  C02  protocol,  the 
smaller  ApHs  for  SGLTl  vs.  H20  oocytes  could  reflect  a 
decrease  in  the  expression  of  other  proteins  or  an  increase  in 
membrane-protein/lipid  ratio. 

Given  our  observations  with  C02,  one  might  ask  whether 
AQP1  also  would  enhance  the  flux  of  a  permeant  weak  acid  like 
butyric  acid  (33).  While  confirming  that  AQP1  increases  the 
ApHs  in  oocytes  exposed  to  C02  (Fig.  1  H  and  /),  we  found  that 
the  channel  has  no  effect  either  on  ApHs  (Fig.  1  /  and  K)  or  tpMs 
(Fig.  SI  C  and  /))  in  oocytes  exposed  to  butyric  acid. 

In  principle,  AQP1  could  enhance  the  C02-induced  pHs  spike 
in  Fig  1  Cor/7  not  hecause  AQP1  is  a  CO2  channel,  hut  hecause 
it  has  unanticipated  carbonic  anhydrase  (CA)  activity  and  thus 
catalyzes  the  extracellular  reaction  HCOf  +  H^  — >C02  +  H20 
(see  Fig.  \A).  To  test  the  CA  hypothesis,  we  injected  oocytes 
with  cRNA  encoding  CA  IV,  coupled  via  a  GPI  linkage  to  the 


outer  surface  of  the  memhrane.  Increased  CA-1V  expression 
causes  a  graded  increase  in  ApHs  (Fig.  S3),  and  we  determined 
the  dose  of  CA-1V  cRNA  that  produces  the  same  ApHs  as  a 
typical  AQP1  oocyte.  Fig.  2  shows  that  membrane  preparations 
of  oocytes  injected  with  this  dose  of  CA-1V  cRNA,  compared 
with  those  from  IFO  oocytes,  require  a  much  shorter  time  to 
achieve  a  pH  endpoint  in  a  CA  assay.  However,  memhrane 
preparations  of  AQP1  oocytes  are  indistinguishable  from  those 


100  pMMTZ 

Fig.  2.  Carbonic  anhydrase  activities  of  Xenopus  oocytes  CA  activity  was 
determined  from  membrane  preparations  of  100  oocytes  injected  with  H2O, 
0  25  ng  of  cRNA  encoding  hCA  IV,  or  25  ng  of  cRNA  encoding  hAQPl  We 
divided  each  membrane  preparation  into  aliquots  containing  20  ^g  of  total 
protein  and  performed  a  colorimetric  CA  assay  on  each  aliquot  (see  5/  Text). 
The  sample  mixtures  containing  CA-IV  were  run  ±  100  ;iM  methazolamide 
(MTZ),  a  CA  inhibitor  Each  N  refers  to  a  CA  assay  on  1  aliquot. 
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Fig.  3.  Surface  pH  changes  caused  by  CO2  and  NH3  influx  in  oocytes  expressing  gas-channel  proteins  (A-E)  Typical  pHs  transients  from  oocytes  injected  with 
H20  or  expressing  AQPI,  AQP4,  AQP5,  WT  AmtB  or  its  inactive  D160A  mutant,  or  WT  RhAG  or  its  inactive  D167A  mutant.  The  protocol  was  the  same  as  in  Fig. 
1.  (F-J)  Summary  of  extreme  excursions  of  pHs  (ApHs)  caused  by  C02  influx.  Each  image  represents  mean  values  for  day-matched  oocytes.  ( K-0 )  Summary  of  ApHs 
caused  by  NH3  influx  Each  image  ( F-0 )  represents  mean  values  for  day-matched  oocytes.  Some  H20  oocytes  served  as  controls  in  more  than  1  panel  (total  number 
of  H20  oocytes:  S4  for  C02,  61  for  NH3)  Values  are  means  t  SE,  with  numbers  of  oocytes  in  parentheses.  For  F-H  and  K-M,  statistical  comparisons  were  made 
using  unpaired  2  tailed  t  tests  For  /and  J  and  A/ and  O,  statistical  comparisons  were  made  using  1-way  ANOVAsfor  3  groups,  followed  by  Student-Newman-Keuls 
analyses. 


of  H2O,  allowing  us  to  rule  out  the  CA  hypothesis.  Thus,  on  the 
basis  of  ApHs  and  rplis  measurements,  AQP1,  but  not  SGLT1, 
NKCC1,  and  PepTI,  acts  as  a  channel  for  CO2  and  NHv 

Comparison  of  AQPs  with  AmtB  and  RhAG.  Using  the  same  protocol 
shown  in  Fig.  I  C-E ,  we  systematically  examined  the  effects  of 
sequential  exposures  to  CO2/HCO3  and  then  NH3/NH4  on 
oocytes  expressing  AQP1,  AQP4,  AQP5,  AmtB,  or  RhAG.  Fig. 
3 A  shows  again  that  CO2  and  NH3  elicit  larger  pHs  spikes  in 
AQP1  (green)  than  in  day-matched  H2O  (orange)  oocytes.  Both 
AQP4  (Fig.  3 B)  and  AQP5  (Fig.  3C)  enhance  the  pHs  spike  with 
CO2  but  not  with  NH3.  Both  AmtB  (Fig.  3D)  and  RhAG  (Fig. 
3 E)  enhance  the  pH$  spike  with  COj  but  arc  especially  effective 
with  NH3.  Asp1611  is  vital  for  AmtB  activity  (19,  34),  and  the 
homologous  Asp167  is  required  by  RhAG  (35).  We  found  that  the 
inactive  D 1 60A  mutant  of  AmtB  (34,36)  and  the  inactive  D167N 
mutant  of  RhAG  (35)  are  inactive  as  either  CO2  or  NH3 
channels,  presumably  beeause  the  mutations  cause  major  struc¬ 
tural  changes  (35).  Fugs.  S3  and  S4  A  and  /)  show  full-length 
experiments  similar  to  those  in  Fig.  3  A-E. 

In  the  above  experiments,  all  AmtB  and  RhAG  constructs 
were  C-terminally  tagged  with  EGFP  (enhanced  GFP),  and 
fluorescence  measurements  confirmed  trafficking  to  near  the 
oocyte  surface.  The  tagged  and  untagged  constructs  yielded 
identical  results  in  pHs  assays  (Fig.  S4). 

Fig.  3  E-J  arc  analogous  to  Fig.  IF,  except  that  in  Fig.  3  F-J, 
we  pair  each  oocyte  expressing  a  WT  or  mutant  channel  with  its 
day-matched  I^O-injectcd  control.  Each  WT  protein  yields  a 
ApHs  that  is  significantly  greater  than  the  H20  control  or  (as 
applicable)  the  mutant  protein.  Moreover,  the  ApHs  values  of 
the  mutants  are  not  different  from  the  corresponding  H^O 
oocytes.  Fig.  3  K-O  is  a  summary  of  the  NH3  data.  The  results 
are  comparable  to  the  CO2  data,  except  that  the  magnitudes  of 
ApHs  for  AQP4  and  AQP5  in  the  NH3  protocol  are  not  different 
from  those  of  their  corresponding  H20-injected  controls.  Fig.  S5 
shows  that  the  relationship  for  the  TpHs  values  is  the  inverse  of 
that  for  the  ApHs  values  in  Fig.  3.  Thus,  each  of  the  proteins, 
AQPI,  AQP4,  AQP5,  AmtB,  and  RhAG,  is  permeable  to  CO2. 
Moreover,  AQPI,  AmtB,  and  RhAG,  but  not  AQP4  or  AQP5, 
are  permeable  to  NH3. 


Cleavage  of  AmtB  Signal  Peptide  by  Oocytes.  Our  AmtB  cDNA 
encodes  a  signal  peptide  that  Escherichia  coli  naturally  cleaves 
(37),  so  that  the  new  N  terminus  is  extracellular.  The  A22K  point 
mutation  in  AmtB  prevents  the  cleavage  in  E.  coli,  although  the 
mutant  AmtB  still  forms  trimers  and  is  active  (37).  To  verify  that 
Xenopus  oocytes  also  cleave  the  signal  peptide,  we  added  a 
C-terminal  His  tag  to  WT  AmtB,  A22K-AmtB,  and  an  AmtB 
variant  with  the  signal  peptide  truncated.  Fig.  4  shows  Western 
blots  of  plasma-membrane  preparations  from  oocytes  expressing 
the  3  constructs.  The  molecular  mass  of  the  major  band  of  WT 
AmtB  is  appropriately  less  than  that  of  A22K-ArntB,  but  the 
same  as  truncated  AmtB.  Thus,  oocytes  do  indeed  cleave  the 
signal  peptide  of  WT  AmtB.  Densitometry  indicates  that  >90% 
of  the  AmtB  in  oocytes  is  appropriately  cleaved.  Additional  data 
reveal  that  AmtB  His  is  active  as  both  a  CO2  and  an  NH? 
channel. 

Pf  in  Oocytes  Expressing  Different  Channels.  So  that  we  could  relate 
our  CO2  and  NIT  data  to  the  wealth  of  information  on  the 
osmotic  water  permeability  (Pf)  of  AQP-expressing  oocytes,  we 
determined  Pf  for  each  AQP  oocyte  and  its  day-matched  control 
from  the  dataset  in  Fig.  3.  As  summarized  by  the  3  pairs  of  bars 
on  the  left  side  of  Fig.  5,  the  mean  Pf  value  for  each  AQP  was 
significantly  and  substantially  greater  than  the  matched  controls. 
We  separately  assessed  Pf  for  matched  oocytes  expressing  AQPI , 
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Fig.  4.  Western  blots  testing  cleavage  of  the  AmtB  signal  peptide  in  Xenopus 
oocytes.  (A)  Wild-type  AmtB  vs.  uncleavable  A22K  mutant  (B)  Wild  type  AmtB 
vs  AmtB  with  truncated  signal  peptide  Data  are  representative  of  4  similar 
experiments  All  constructs  were  His  tagged  at  the  C  terminus  and  detected 
with  an  anti-His  antibody. 
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Fig.  S.  Osmotic  water  permeabilities  of  Xenopus  oocytes.  Pf  (cm/s)  of  oocytes 
injected  with  H2O  or  cRNA  encoding  AQPI,  AQP4,  AQP5,  AmtB,  or  RhAG 
Values  are  means  t  5E,  with  numbers  of  oocytes  In  parentheses  Statistical 
comparison  between  H20  vs  AQP  oocytes  were  made  using  unpaired  2  tai led 
t  tests.  Statistical  comparisons  among  the  4  groups  were  made  using  a  1-way 
ANOVA,  followed  by  Student-Newman-Keuls  analyses. 


AmtB,  or  RhAG  vs.  day-matchcd  H2O  oocytes.  The  right  side  of 
Fig.  5  shows  that  only  for  AQP1  oocytes  was  the  mean  Pf  value 
significantly  greater  than  that  for  FFO  oocytes;  the  Pf  values  for 
ITO,  AmtB,  and  RhAG  oocytes  were  not  significantly  different 
from  one  another.  Thus,  despite  the  hypothesized  presence  of 
FBO  in  the  NHj  pore  of  AmtB  (20,  22),  the  2  Rh  proteins  do  not 
function  as  water  channels. 

Discussion 

Channel-Dependent  ApHs  and  Pf  Values.  Given  our  experimental 
design,  the  magnitude  of  ApHs  is  a  semiquantitative  index  of 
hoth  the  flux  of,  and  membrane  pcrmcahility  to,  CO2  or  NH3. 
For  CO2,  the  same  is  true  of  tpiiv  Note  that  the  quantitative 
relationship  between  Apl  ls  on  the  one  hand  and  absolute  flux  or 
permeability  on  the  other  is  likely  to  be  different  for  CO2  vs. 
NH.v  The  portion  of  the  CC^-induccd  ApHs  signal  that  we  can 
ascrihe  to  a  particular  channel  is  the  difference  hetween  the 
ApHs  of  each  channel-expressing  oocyte  (e.g.,  green  record  in 
Fig.  3 A)  and  the  ApHs  of  its  day-matched  tBO-mjected  control 
(e.g.,  orange  record  in  Fig.  3/1).  Fig.  6A  summarizes  these 
differences,  the  channel-dependent  signal  (ApHs*)co2»  for  the 
CO2  data,  computed  oocyte  hy  oocyte.  Similarly,  Fig.  6/?  sum¬ 
marizes  the  analogous  differences,  the  channel-specific  signal 
(Apl  Is*)nh3»  f°r  the  NIB  data.  Note  that  the  mean  (ApHs*)iMm 
values  for  AQP4  and  AQP5  arc  not  significantly  different  from 
zero. 

Ratios  of  Indices  of  Permeability.  Because  we  do  not  know  the 
number  of  AQP  molecules  at  the  plasma  membrane  in  each 
oocyte,  it  is  impossible  to  normalize  our  (ApHs*)co2  or 
(ApHs* )niu  data  to  protein  abundance.  However,  for  each  AQP 
oocyte,  wc  also  have  a  channel-dependent  Pf  (Pf),  summarized 
in  Fig.  6C.  For  each  oocyte,  wc  divided  (ApHs*)co2  or 
(Apl  Is*)niii  hy  Pt* .  Fig.  6 D  summarizes  these  mean  values,  which 
represent  semiquantitative  indices  of  the  CO2/H2O  or  NH3/H2O 
permeahility  ratios.  By  a  factor  of  2,  AQP5  has  the  highest 
(ApHs*)co2/Pf>  and  the  values  for  AQP1  and  AQP4  arc  indis¬ 
tinguishable. 

Because  we  have  both  (ApHs*)co2  and  (ApHs*)Nin  for  each 
of  a  large  numher  of  oocytes,  it  is  also  possihlc  to  compute  the 
ratio  (ApHs*)co2/(ApHs*)NH3,  a  relative  index  of  the  CO2/NH3 
permeahility  ratio,  for  AQP1,  AmtB,  and  RhAG.  Fig.  6 E 
summarizes  these  values.  Because  (ApHs*)NH3  for  AQP4  and 
AQP5  do  not  differ  from  zero,  the  ratios  for  these  proteins  arc 
theoretically  infinite.  Thus,  among  the  channels  tested,  AQP4 
and  AQP5  have  the  highest  CO2/NH3  permeahility  ratios,  fol¬ 
lowed  hy  AQP1,  AmtB.  and  RhAG.  Conversely,  RhAG  has  the 
highest  NH3/CO2  permeability  ratio  (sec  Fig.  $6). 
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Fig.  6.  Comparison  of  channel-dependent  properties.  (A)  Indices  of  channel- 
dependent  C02  permeability  For  each  ApHs  from  a  channel  expressing  00 
cyte,  we  subtracted  the  mean,  day-matched  ApHs  for  H20  oocytes.  Bars 
represent  mean  subtracted  values,  the  channel-dependent  ApH5  for  CO2,  or 
(ApHS*)co2  Note  Oocytes  in  A  are  the  same  as  those  in  B  and  C.  (B)  Indices  of 
channel  dependent  NH?  permeability  We  computed  (ApHs*)nh3  using  the 
same  approach  as  in  A.  (Q  Channel-dependent  water  permeabilities.  For  each 
Pf  from  a  channel-expressing  oocyte,  we  subtracted  the  mean,  day  matched  Pf 
for  H20  oocytes  Bars  represent  mean  subtracted  values,  the  channel- 
dependent  Pf,  or  Pj.  (D)  Indices  of  channel-dependent  C02  and  NH3  perme¬ 
ability,  normalized  to  Pf.  For  each  oocyte,  we  divided  (ApHs*)C02  and 
(ApHs*)nh3  by  its  Pf.  (£)  Indices  of  gas  selectivity.  For  each  oocyte  expressing 
AQP1,  AmtB,  or  RhAG,  we  divided  (ApHs*)Co2  by  ~(ApHs*)NH3.  Because 
(ApHs*)NH3  was  not  significantly  different  from  zero  for  oocytes  expressing 
AQP4  or  AQP5,  we  represent  these  ratios  as  "infinity."  Statistical  comparisons 
were  made  using  1  -way  ANOVAs  for  3  groups,  followed  by  Student-Newman- 
Keuls  analyses 


Thus,  compared  with  the  AQPs  tested,  the  Rh-like  proteins 
tested  arc  relatively  more  selective  for  NH3,  whereas  the  AQPs 
arc  relatively  more  selective  for  CCK 

Significance.  Our  data  demonstrate  that  channel  proteins  can 
exhibit  gas  selectivity  hy  channel  proteins.  The  basis  of  the 
selectivity  is  probably  not  the  size  of  the  transiting  molecules — 
H20,  CO2,  and  NH3  have  similar  minimum  diameters — hut 
rather  their  chemistries  and  the  chemistries  of  the  monomeric 
pores  vs.  the  pore  at  the  center  of  the  multimcrs.  The  electronic 
configuration  of  NH3  is  identical  to  that  of  FTO,  which  moves 
exclusively  through  the  monomeric  aquaporcs  of  AQPI.  Thus, 
the  hydrophilic  NH3  probably  also  moves  exclusively  through  the 
monomeric  aquaporcs  of  AQPI  and  through  the  monomeric 
ammonia  pores  of  AmtB  and  RhAG.  The  less  hydrophilic  CO2, 
however,  could  move  through  the  hydrophohic  central  pores  of 
all  5  channels.  NO  is  also  known  to  move  through  AQPI  (8),  and 
indirect  evidence  is  consistent  with  the  idea  that  O2  moves 
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through  AQPl  (38).  We  suggest  that  the  hydrophobic  NO  and 
O:  move  through  the  central  pores.  Crystallographic  data  show 
that  xenon  can  enter  the  central  pore  of  the  bacterial  Rh50  (21 ). 
Moreover,  in  the  case  of  AQP1,  molecular-dynamics  simulations 
show  that  CO2  could  penetrate  the  4  aquapores  or,  with  greater 
ease,  the  central  pore  (7).  We  hypothesize  that  the  CO2/NH3 
selectivities  that  we  observe  reflect  the  relative  permeabilities  of 
the  2  gases  through  the  monomeric  vs.  the  central  pores  of  the 
5  channels  wc  studied 

In  a  cell  like  the  human  RBC,  whose  plasma-membrane  lipid 
has  an  intrinsically  low  gas  permeability  (6),  the  gas  selectivity  of 
AQP1  and  the  Rh  complex  would  provide  control  over  dissolved 
gases  crossing  the  membrane.  The  NH3  permeability  of  AQP1 
and  RhAG  could  enhance  the  ability  of  RBCs  to  pick  up  NHi  in 
various  tissues  (where  the  NH3  gradient  would  favor  NIH 
uptake)  and  then  to  off-load  it  in  the  liver  (where  the  gradient 
would  favor  NHi  efflux  from  RBCs  and  uptake  by  hcpatocytcs) 
for  NH3  detoxification.  In  the  hypertonic  renal  medulla,  this 
NH3  permeability  could  reduce  the  reflection  coefficient 
for  N1H  and  thereby  reduce  cell-volume  changes.  However,  the 
low  NH3  permeability  of  AQP4  could  protect  the  brain  from 
rising  blood  levels  of  NH3,  while  still  allowing  CO2,  and  perhaps 
NO  and  O2,  to  pass. 

Materials  and  Methods 

Molecular  Biology,  AQPs.  Human  AQPl  cDNA  (GenBank  accession  no. 
NMJ 98098),  cDNA  encoding  the  rat  AQP4/M23  splice  variant  (GenBank  ac¬ 
cession  no.  NM  012825),  and  human  AQPS  cDNA  (GenBank  accession  no 
NM  012779)  were  gifts  of  Peter  Agre  (Johns  Hopkins  University,  Baltimore) 
AmtB.  We  cloned  E  co li  AmtB  cDNA  (GenBank  accession  no.  ECU40429)  by  PCR 
from  genomic  DNAand  subcloned  the  — 1  3-kb  PCR  product  into  the  Xenopus 
expression  vector  pGH19  (39).  Using  PCR,  we  created  an  additional  construct 
in  which  we  replaced  the  nucleotides  encoding  the  signal  sequence  (i.e.,  first 
22  residues)  with  ATG.  At  the  3  end  of  some  constructs,  we  added  in-frame 
cDNA  encoding  either  EGFP  (Clontech,  ref.  40)  or  a  His  tag. 

RhAG  Human  RhAG  cDNA  in  pT7T5  (GenBank  accession  no  NM  000324,  a  gift 
of  Baya  Ch£rif-Zahar,  Universite  Rene  Descartes,  IN5ERM,  Paris)  (41)  was 
subdoned  into  pGH19.  We  tagged  RhAG  at  its  3  end  with  EGFP 
Other  tDNAs.  Rabbit  NKCC2  (42)  was  a  gift  of  Biff  Forbush  (Yale  University,  New 
Haven,  CT)  SGLTl  (43)  and  PepTl  (44)  were  gifts  of  Matthias  Hediger  (Brigham 
and  Women’s  Hospital  and  Harvard  Medical  School,  Boston). 

Site-Directed  Mutagenesis.  We  used  the  QuikChange  Site  Directed  Mutagen 
esis  Kit  (Stratagene),  following  the  manufacturer's  instructions 

cRNA  Preparation.  We  generated  cRNA  using  the  Message  Machine  kit  (Am- 
bion)  and,  unless  otherwise  stated,  injected  oocytes  with  50  nL  of  0  5  ng/nL  of 
cRNA  or  25  nL  of  1  ng/nL  of  cRNA. 

Western  Blot  Analysis.  Plasma-membrane  proteins  were  prepared  from  oo¬ 
cytes  (45),  separated  on  a  13%  SDS  polyacrylamide  gel,  blotted  on  a  PVDF 
membrane,  probed  with  a  monoclonal  anti-His  antibody  (Catalog  no. 
70796-3,  Novagen),  and  detected  using  ECL  plus  Western  Blotting  Detection 
Reagents  (Amersham  Biosciences) 

Solutions  for  Physiological  Assays.  The  ND96  solution  contained  96  mM  NaCl,  2 
mM  KCI,  1  mM  MgCU,  1.8  mM  CaCI^  and  5  mM  Hepes,  pH  7.50,  osmolality  195 
mOsm.  For  Pf  assays,  we  used  a  hypotonic  ND96  variant  (100  mOsm)  that  con¬ 
tained  only  43  mM  NaCl.  The  CO2/HCO3  solution  was  identical  to  ND96  except 
that  33  mM  NaHC03  replaced  33  mM  NaCl,  and  thesolution  was  bubbled  with  5% 
COybalance  02-  The  NH*/NH4+  solution  was  a  variant  of  ND96  in  which  we 


replaced  0.S  mM  NaCl  with  0.5  mM  NH4CI.  The  butyrate  solution  was  a  variant  of 
ND96  in  which  we  replace  30  mM  Na-butyrate  with  30  mM  NaCl 

Carbonic-Anhydrase  Assay.  Carbonic-anhydrase  activity  was  assessed  in  20  pg 
of  membrane  preparation  of  CA-IV  or  AQPl  oocytes  using  a  colorimetric 
technique  (46)  The  assay  measures  the  rate  at  which  the  pH  of  a  weakly 
buffered  alkaline  solution  (imidazole-Tris,  50%  C02,  with  p-nitrophenol  as 
indicator  at  0  °C),  falls  in  the  presence  or  absence  of  CA,  noted  by  a  color 
change  from  yellow  to  clear,  due  to  the  reaction  C02  +  H20  —  HC03  +  H\ 

Measurement  of  Oocyte  Water  Permeability.  We  used  a  volumetric  assay  (47, 
48)  to  measure  osmotic  water  permeability  (Pf).  Briefly,  afterdropping  oocytes 
into  a  Petri  dish  containing  the  hypotonic  solution,  we  acquired  video  images 
every  1-2  s,  obtaining  the  time  course  of  the  projection  area  of  the  oocyte. 
Assuming  the  oocyte  to  be  a  sphere,  and  the  true  surface  area  (S)  to  be  8  fold 
greater  than  the  idealized  area  (49),  we  computed  Pf  as: 

V  d(V/V^ 

0  dt 

Pf  "  S  •  AOsm  •  V* 

where  V0  is  initial  oocyte  volume,  d(V/V0)/dt  is  the  maximal  fractional  rate  of 
volume  increase,  AOsm  is  the  osmotic  gradient  across  the  membrane,  and  Vw 
is  the  molar  volume  of  water. 

Measurement  of  Surface  pH  We  used  microelectrodes  to  measure  pHs  (6,  50). 
Briefly,  the  pH  electrode  had  a  tip  diameter  of  15  pm,  was  filled  at  its  tip  with  H  + 
ionophore  mixture  B  (Catalog  no  95293,  Fluka),  and  was  connected  to  a  FD223 
electrometer  (World  Precision  Instruments).  The  extracellular  reference  electrode 
was  a  glass  micropipette  filled  with  3  M  KCI  and  connected  via  a  calomel  half  cell 
to  a  7S0  electrometer  (World  Precision  Instruments).  The  extracellular  solution 
flowed  at  3  mL/min,  and  the  sampling  rate  was  1  per  500  ms  Using  an  MPC-200 
system  micromanipulator  (Sutter  Instrument),  we  positioned  the  pHs  electrode 
tip  either  in  the  bulk  extracellular  fluid  or  dimpling  —40  fxm  onto  the  oocyte 
surface,  in  the  "shadow"  of  the  oocyte.  Although  not  displayed  in  the  figures, 
membrane  potential  (Vm)  and  intracellular  pH  were  also  monitored  (see  5/  Text). 
All  oocytes  had  initial  V™  values  at  least  as  negative  as  -40  mV 

We  verified  delivery  of  EGFP-tagged  proteins  to  a  region  near  the  plasma 
membrane  by  using  a  96-well  plate  reader  (BMG  Labtechnologies)  to  assess 
fluorescence  (40). 

Data  Analysis.  Before  applying  CO2/HCO3  or  NH3/NH4,  we  computed  pH$from 
the  preceding  calibration,  withthe  electrodetip  inthe  bulk  phaseof  the  ND96 
solution  (pH  7.50).  After  applying  CCVHC03  or  NH3/NH4 ,  we  computed  pHs 
from  a  second  calibration  in  the  bulk  phase  of  the  new  solution  (also  pH  7  50) 
The  maximum  pH$  excursion  (ApH$)  was  taken  as  the  maximum  pH$  after  the 
application  of  C02/HC03  (or  the  minimum  pHs  after  application  of  NH3/NH4 ) 
minus  the  pHs  prevailing  just  before  the  solution  change  from  ND96, 

Statistics.  Data  are  presented  as  mean  1.  5EM.  To  compare  the  difference 
between  2  means,  we  performed  Student's  f  tests  (two  tails).  To  compare  more 
than  2  means,  we  performed  a  1-way  ANOVA  followed  by  a  Dunnett's  or  a 
Student-Newman-Keuls  posthoc  analysis,  using  KaleidaGraph  (Version  4, 
Synergy  Software)  P  <  0  05  was  considered  significant 
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